Compute the Expected F_{ST}

Description

Compute the expected F_ST in a finite island model.

Usage

1
	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

1
2
3
4
## 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)