Learn R Programming

polyfreqs (version 1.0.2)

point_Hexp: Estimation of expected heterozygosity

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.

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)

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.