DirichNormHWE: Evaluates the normalizing constant under the HWE model, with...

In HWEBayes: Bayesian investigation of Hardy-Weinberg Equilibrium via estimation and testing.

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 ([e-mail]).

References

Wakefield, J. (2010). Bayesian methods for examining Hardy-Weinberg equilibrium. Biometrics; Vol 66:257-65

See Also

DirichNormSat

Examples

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

HWEBayes documentation built on May 2, 2019, 11:07 a.m.