Compute the *F_ST* for a sample of genes.

1 |

`sample` |
an object generated by the |

Once the `sim.inference.model`

or by the `sim.coalescent`

command lines have been executed, `fst`

can be used to compute the Weir and Cockerham's (1984) estimate of *F_ST* for that sample. The multilocus estimator of *F_ST* is computed as in Genepop (Rousset 2008).

Weir and Cockerham's (1984) estimate of *F_ST* for the sample of genes contained in the object ‘sample’

Rousset, F. (2008) genepop'007: a complete re-implementation of the genepop software for Windows and Linux, *Molecular Ecology Notes* **8**: 103-106.

Weir, B. S. and Cockerham, C. C. (1984) Estimating *F*-statistics for the analysis of population structure, *Evolution* **38**: 1358-1370.

1 2 3 4 5 6 7 8 9 10 | ```
## This is to simulate a sample of genes (at a single locus), using the inference model, with
## 50 genes collected in each of 10 sampled demes. In this example, the product of
## twice the effective population size and migration rate is 2,
## and the frequency of allele A in the migrant pool is 0.5
sample <- sim.inference.model(number.of.sampled.demes = 10,sample.sizes = 50,M = 2,pi = 0.5)
## This is to compute Weir and Cockerham's (1984) estimate for that sample
fst(sample)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.