DirichNormSat: Evaluates the normalizing constant (as used in the...

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

References

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

See Also

DirichSampSat, DirichSampHWE, DirichNormSat, DirichNormHWE, HWEDirichBF2

Examples

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

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