Description Usage Arguments Value Author(s) References See Also Examples
View source: R/pop_specificFST.R
This function estimates genome-wide poppulation-specific FST based on Weir and Goudet's Method (Weir&Goudet 2017) from a GENEPOP data object (Rousset 2008). Missing genotype values in the GENEPOP file ("0000" or "000000") are simply ignored.
1 | pop_specificFST(popdata, cov = FALSE)
|
popdata |
Population data object created by |
cov |
A Logical argument indicating whether variance-covariance matrix of estimated FST should be calculated. |
fst |
Estimated genome-wide population-specific FST and standard error. |
cov |
Variance-covariance matrix of estimated FST |
Reiichiro Nakamichi, Shuichi Kitada, Hirohisa Kishino
Rousset F (2008) Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol. Ecol. Resources, 8, 103-106.
Weir BS, Goudet J (2017) A Unified Characterization of Population Structure and Relatedness. Genetics, 206, 2085-2103.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | # Example of GENEPOP file
data(jsmackerel)
jsm.ms.genepop.file <- tempfile()
jsm.popname.file <- tempfile()
cat(jsmackerel$MS.genepop, file=jsm.ms.genepop.file, sep="\n")
cat(jsmackerel$popname, file=jsm.popname.file, sep=" ")
# Data load
# Prepare your GENEPOP file and population name file in the working directory.
# Replace "jsm.ms.genepop.file" and "jsm.popname.file" by your file names.
popdata <- read.GENEPOP(genepop=jsm.ms.genepop.file, popname=jsm.popname.file)
# FST estimation
result.popspFST <- pop_specificFST(popdata)
print(result.popspFST$fst)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.