DirichNormHWE {HWEBayes} | R Documentation |
Function to evaluate the normalizing constant given a conjugate Dirichlet prior and the HWE model.
DirichNormHWE(nvec, bvec0)
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. |
The normalizing constant.
Jon Wakefield (jonno@u.washington).
Wakefield, J. (2009). Bayesian methods for examining Hardy-Weinberg equilibrium.
DirichNormSat
data(DiabRecess) DirichNormHWE(nvec=DiabRecess,bvec0=rep(1,4))