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

Share:

Description

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

Usage

1
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. (2010). Bayesian methods for examining Hardy-Weinberg equilibrium. Biometrics; Vol 66:257-65

See Also

DirichSampSat, DirichSampHWE, DirichNormSat, DirichNormHWE, HWEDirichBF2

Examples

1
2