DirichNormSat {HWEBayes}R Documentation

Evaluates the normalizing constant (as used in the denominator of a Bayes factor) for a conjugate prior

Description

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

Usage

DirichNormSat(nvec, bvec)

Arguments

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

Value

The normalizing constant.

Author(s)

Jon Wakefield (jonno@u.washington.edu)

References

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

See Also

DirichSampSat, DirichSampHWE, DirichNormSat, DirichNormHWE, HWEDirichBF2

Examples

data(DiabRecess)
DirichNormSat(nvec=DiabRecess,bvec=rep(1,10))

[Package HWEBayes version 1.0 Index]