calcFStats: Calculate locus-wise and population theta = F_ST, F = F_IT,...

In relSim: Relative Simulator

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

relSim documentation built on May 29, 2024, 5:33 a.m.