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; <DOI:10.2307/2411227>. Functions are also given to test via randomisations the significance of each F and variance components, using the likelihood-ratio statistics G.
copied from cf-staging / r-hierfstat