FstBoot: Bootstrap sampler of Fst
In FinePop: Fine-Scale Population Analysis

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

This function provides bootstrapped estimators of Fst to evaluate the environmental effects on the genetic diversity. See Details.

1	FstBoot(popdata, fst.method = "EBFST", bsrep = 100, log.bs = F, locus = F)

`popdata`	Genotype data object of populations created by read.genepop function from a GENEPOP file.
`fst.method`	A character value specifying the Fst estimation method to be used. Currently, "EBFST", "EBGstH", "EBDJ", "GstN", "GstNC", "GstH", "DJ" and "thetaWC.pair" are available.
`bsrep`	A positive integer value specifying the trial times of bootstrapping.
`log.bs`	A logical value specifying whether the bootstrapped data of each trial should be saved. If TRUE, GENEPOP format files named "gtdata_bsXXX.txt" (XXX=trial number) are saved in the working directory.
`locus`	A Logical argument indicating whether locus-specific Fst values should be calculated.

FinePop provides a method for regression analyses of the pairwise Fst values against geographical distance and the differences to examine the effect of environmental variables on population differentiation (Kitada et. al 2017). First, FstBoot function resamples locations with replacement, and then, we also resample the member individuals with replacement from the sampled populations. It calculates pairwise Fst for each bootstrap sample. Second, FstEnv function estimates regression coefficients (lm function) of for the Fst values for each iteration. It then computes the standard deviation of the regression coefficients, Z-values and P-values of each regression coefficient. All possible model combinations for the environmental explanatory variables were examined, including their interactions. The best fit model with the minimum information criterion (TIC, Takeuchi 1976, Burnham & Anderson 2002) is selected. Performance for detecting environmental effects on population structuring is evaluated by the R2 value.

`bs.pop.list`	List of subpopulations in bootstrapped data
`bs.fst.list`	List of genome-wide pairwise Fst matrices for bootstrapped data.
`org.fst`	Genome-wide pairwise Fst matrix for original data.
`bs.fst.list.locus`	List of locus-specific pairwise Fst matrices for bootstrapped data. (If locus = TRUE)
`org.fst.locus`	Locus-specific pairwise Fst matrix for original data. (If locus = TRUE)

Reiichiro Nakamichi, Hirohisa Kishino, Shuichi Kitada

Burnham KP, Anderson DR (2002) Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. Springer, New York.

Kitada S, Nakamichi R, Kishino H (2017) The empirical Bayes estimators of fine-scale population structure in high gene flow species. Mol. Ecol. Resources, DOI: 10.1111/1755-0998.12663.

Takeuchi K (1976) Distribution of information statistics and criteria for adequacy of Models. Mathematical Science, 153, 12-18 (in Japanese).

read.genepop, FstEnv, EBFST, EBGstH,EBDJ, GstN, GstNC, GstH, DJ, thetaWC.pair, herring

# Example of genotypic and environmental dataset
data(herring)

# Data bootstrapping and Fst estimation
# fstbs <- FstBoot(herring$popdata)

# Effects of environmental factors on genetic differentiation
# fstenv <- FstEnv(fstbs, herring$environment, herring$distance)

# Since these calculations are too heavy, pre-caluculated results are included in this dataset.
fstbs <- herring$fst.bootstrap
fstenv <- herring$fst.env
summary(fstenv)