EBDJ: Empirical Bayes estimator of Jost's D

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/EBDJ.R

Description

This function estimates pairwise D (Jost 2008) among subpopulations using empirical Bayes method (Kitada et al. 2007). This function accepts two types of data object, GENEPOP data (Rousset 2008) and allele (haplotype) frequency data (Kitada et al. 2007). Missing genotype values in the GENEPOP file ("0000" or "000000") are simply ignored.

Usage

1
EBDJ(popdata, num.iter=100)

Arguments

popdata

Genotype data object of populations created by read.genepop function from a GENEPOP file. Allele (haplotype) frequency data object created by read.frequency function from a frequency format file also is acceptable.

num.iter

A positive integer value specifying the number of iterations in empirical Bayes simulation.

Details

Frequency format file is a plain text file containing allele (haplotype) count data. This format is mainly for mitochondrial DNA (mtDNA) haplotype frequency data, however nuclear DNA (nDNA) data also is applicable. In the data object created by read.frequency function, "number of samples" means haplotype count. Therefore, it equals the number of individuals in mtDNA data, however it is the twice of the number of individuals in nDNA data. First part of the frequency format file is the number of subpopulations, second part is the number of loci, and latter parts are [population x allele] matrices of the observed allele (haplotype) counts at each locus. Two examples of frequency format files are attached in this package. See jsmackerel.

Value

Matrix of estimated pairwise Jost's D.

Author(s)

Reiichiro Nakamichi, Hirohisa Kishino, Shuichi Kitada

References

Jost L (2008) Gst and its relatives do not measure differentiation. Molecular Ecology, 17, 4015-4026.

Kitada S, Kitakado T, Kishino H (2007) Empirical Bayes inference of pairwise FST and its distribution in the genome. Genetics, 177, 861-873.

Rousset F (2008) Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol. Ecol. Resources, 8, 103-106.

See Also

read.genepop, read.frequency

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# Example of GENEPOP file
data(jsmackerel)
cat(jsmackerel$MS.genepop, file="JSM_MS_genepop.txt", sep="\n")
cat(jsmackerel$popname, file="JSM_popname.txt", sep=" ")

# Data load
# Prepare your GENEPOP file and population name file in the working directory
# (Here, these files were provided as "JSM_MS_genepop.txt" and "JSM_popname.txt".)
popdata <- read.genepop(genepop="JSM_MS_genepop.txt", popname="JSM_popname.txt")

# Jost's D estimation
result.EBDJ <- EBDJ(popdata)
write.csv(result.EBDJ, "result_EBDJ.csv", na="")
print(as.dist(result.EBDJ))

Example output

Computing global differenciation with ML method... done.
Computing pairwise differenciation with EB method... pop 1:21:31:41:51:61:71:82:32:42:52:62:72:83:43:53:63:73:84:54:64:74:85:65:75:86:76:87:8done.
!!! FINISH !!!
            OS        HM2         BS        HI1        HI2        HAa
HM2 0.02047266                                                       
BS  0.01901979 0.02111564                                            
HI1 0.01817407 0.02231051 0.02092934                                 
HI2 0.01926271 0.01962844 0.01710790 0.02025211                      
HAa 0.05472664 0.04683020 0.05122904 0.05640782 0.05075922           
HAb 0.12681136 0.14591252 0.12366942 0.13253804 0.14623220 0.15968290
HAr 0.03647948 0.03213550 0.03185911 0.03636448 0.03388405 0.05146561
           HAb
HM2           
BS            
HI1           
HI2           
HAa           
HAb           
HAr 0.12530047

FinePop documentation built on May 22, 2018, 5:07 p.m.