Estimate variance components and hierarchical F-statistics over all loci
Return multilocus estimators of variance components and F-statistics
varcomp.glob(levels=levels,loci=loci,diploid=TRUE)
levels |
a data frame containing the different levels (factors) from the outermost (e.g. region) to the innermost before the individual |
loci |
a data frame containing the different loci |
diploid |
Specify whether the data are coming from diploid or haploid organisms (diploid is the default) |
loc |
The variance components for each locus |
overall |
The variance components summed over all loci |
F |
a matrix of hierarchical F-statistics type-coefficients with the first line corresponding to F_{(n-1)/n},F_{(n-2)/n}...F_{i/n} and the diagonal corresponding to F_{(n-1)/n}, F_{(n-2)/(n-1)},F_{i/2} |
Jerome Goudet DEE, UNIL, CH-1015 Lausanne Switzerland
Weir, B.S. (1996) Genetic Data Analysis II. Sinauer Associates.
Yang, R.C. (1998). Estimating hierarchical F-statistics. Evolution 52(4):950-956
Goudet J. (2005). Hierfstat, a package for R to compute and test variance components and F-statistics. Molecular Ecology Notes. 5:184-186
#load data set data(gtrunchier) attach(gtrunchier) varcomp.glob(data.frame(Locality,Patch),gtrunchier[,-c(1,2)])
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.