Estimation of expected heterozygosity

Usage

point_Hexp(p_samp, genotypes, ploidy)

Arguments

p_samp
A posterior sample of allele frequencies from polyfreqs.
genotypes
Matrix of genotypes sampled during MCMC.
ploidy
The ploidy level of individuals in the population (must be >= 2).

Value

Returns the per locus estimates of expected heterozygosity (per_locus_Hexp)

Description

INTERNAL: Estimates a posterior distribution for the per locus expected heterozygosity using the unbiased estimator of Hardy (2015) and the poterior samples of allele frequencies calculated by polyfreqs.

Details

Posterior distributions for the per locus expected heterozygosity are automatically calculated and returned by the polyfreqs function.

References

Hardy, OJ. 2015. Population genetics of autopolyploids under a mixed mating model and the estimation of selfing rate. Molecular Ecology Resources, doi: 10.1111/1755-0998.12431.