DirichNormHWE {HWEBayes}R Documentation

Evaluates the normalizing constant under the HWE model, with a conjugate prior

Description

Function to evaluate the normalizing constant given a conjugate Dirichlet prior and the HWE model.

Usage

DirichNormHWE(nvec, bvec0)

Arguments

nvec vector of genotype frequencies in the order n_{11}, n_{12},..., n_{1k},n_{22} ..., n_{2k},..., n_{kk}.
bvec0 vector of length k Dirichlet prior parameters, where k is the number of alleles.

Value

The normalizing constant.

Author(s)

Jon Wakefield (jonno@u.washington).

References

Wakefield, J. (2009). Bayesian methods for examining Hardy-Weinberg equilibrium.

See Also

DirichNormSat

Examples

data(DiabRecess)
DirichNormHWE(nvec=DiabRecess,bvec0=rep(1,4))

[Package HWEBayes version 1.0 Index]