Description Usage Arguments Value Author(s) References
Computes F statistics according to Weir and Cockerham's estimators. Missing values are allowed but as of today, the NA code is only treated as an extra allele which might bias the result. This function should not be used on haploid data.
1 | Fstat.output(coordinates,genotypes,ploidy,burnin,path.mcmc)
|
coordinates |
Matrix with one line per individual and two columns |
genotypes |
Genotypes of individuals. A matrix with one line per individual and 2 columns per locus |
burnin |
Integer: number of saved iterations to discard. |
ploidy |
Integer: 1 or 2 (default is 2). Do not use for haploid data. |
path.mcmc |
Path to output files directory |
A list with components
Pairwise.Fis |
A matrix of real numbers estimating the pairwise Fis |
Pairwise.Fst |
A matrix of real numbers estimating the pairwise Fst |
Arnaud Estoup for original code in Turbo Pascal.
Translation in Fortran and interface with R by Gilles Guillot
Weir, B.S. and C.C. Cockerham, Estimating F-statistics for the analysis of population structure, Evolution, 1984, vol. 38, 1358-1370.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.