fst: Generic function to calculate Fst

In MarcoAndrello/MetaPopGen_0.0.8: Simulation of metapopulation genetics

Usage

fst(N)

Arguments

N	A 2-D array. It gives the number of individuals of each genotype (first dimension) in each group (second dimension)

Value

This function calculates the genetic differentiation index Fst among the groups contained in N as ratio of heterozygosities: (Ht - Hs) / Ht. Ht is the expected heterozygosity over the population formed by pooling the groups together. Hs is the average expected heterozygosity over the groups. The average is weighted by the number of individuals in each group.

Examples

# Data: three groups with two alleles
# Genotypic composition of first, second and third groups
N1 <- c(18,299,277)
N2 <- c(284,14,145)
N3 <- c(115,228,254)
N <- array(c(N1,N2,N3),dim=c(3,3),dimnames=list(genotype=c("A1A1","A1A2","A2A2"),group=c(1:3)))
fst(N)

MarcoAndrello/MetaPopGen_0.0.8 documentation built on May 25, 2019, 12:23 p.m.