Allows the estimation of hierarchical F-statistics from haploid or diploid genetic data
with any numbers of levels in the hierarchy, following the algorithm of Yang (Evolution, 1998, 52(4):950-956;
. Functions are also given to test via randomisations the significance of each F and variance components,
using the likelihood-ratio statistics G.