calcFStats {relSim}

Calculate locus-wise and population \theta = F_{ST}, F = F_{IT}, and f = F_{IS} values

Description

This procedure uses the method of Weir and Cockerham to estimate \theta = F_{ST}, F = F_{IT}, and f = F_{IS} for a population with known substructure

Usage

calcFStats(Pop, subPopIdx = NULL)

Arguments

Pop

An object type 'population'

subPopIdx

If this vector is not null, then it must consist of N elements with values from 1 to n_s representing which subpopulation each member of Pop$profiles belongs to. If it is null then it is assumed that the population consists of n_s subpopulations of equal size N_s so that n_s\times N_s = N

Value

A vector of length n_{loci}+1 with locus-wise \theta values and an overall \theta value for the population

Author(s)

James M. Curran

References

Weir, B.S., Genetic Data Analysis II, (1996) p.173–179, Sinauer, Sunderland, MA.

See Also

breedFst

Examples

data(USCaucs)
set.seed(123)
p = breedFst(USCaucs)
fstats = calcFStats(p)
fstats