Fst {snpStats} | R Documentation |
This function calculates the fixation index Fst for each SNP, together with its weight in the overall estimate (as used by the Internation HapMap Consortium).
Fst(snps, group, pairwise=FALSE)
snps |
an object of class |
group |
a factor (or object than can be coerced into a factor),
of length equal to the number of rows of |
pairwise |
if |
See vignette.
A list:
Fst values for each SNP
The weights for combining these into a single index
Uncertain genotypes are treated as missing
David Clayton dc208@cam.ac.uk
## Analysis of some HapMap data data(for.exercise) f <- Fst(snps.10, subject.support$stratum) weighted.mean(f$Fst, f$weight)