expected.fst: Compute the Expected F_{ST}

In infeR: A pedagogical tool to teach the basics of inference in population genetics

Description

Compute the expected F_ST in a finite island model.

Usage

	expected.fst(total.number.of.demes,effective.size,migration.rate,mutation.rate)

Arguments

total.number.of.demes	total number of demes in the island model
effective.size	effective population size of each deme
migration.rate	rate of migration among demes
mutation.rate	mutation rate. The model assumes a 2-alleles model of mutations

Details

Computes the expected F_ST in a finite island model of population structure, with equal sized demes and isotropic migration across demes. The mutation model is a 2-alleles model. This function is based upon eq.(6) in Rousset (1997), adapted for haploids.

Value

Expected F_ST in a finite island model of population structure, with equal sized demes and isotropic migration across demes

References

Rousset, F. (1996) Equilibrium values of measures of population subdivision for stepwise mutation processes, Genetics 142: 1357-1362.

Examples

## This is to get the expected level of population differentiation in an island model with 100 demes,
## each of size 1000, that exchange migrants at a rate of 0.001. The mutation rate is assumed to be 0.000001.

expected.fst(total.number.of.demes = 100,effective.size = 1000,migration.rate = 0.001,mutation.rate = 0.000001)

infeR documentation built on May 2, 2019, 4:45 p.m.