TriangNormHWE {HWEBayes}R Documentation

Evaluates the normalizing constant under the HWE model, for the "triangular" prior distribution

Description

Function to evaluate the normalizing constant given a "triangular" prior and the HWE model, in the k=2 allele case. This prior results from marginalizing the conjugate Dirichlet prior with parameters (1,1,1) on the genotype frequencies under the alternative.

Usage

TriangNormHWE(nvec)

Arguments

nvec vector of genotype frequencies in the order n_{11}, n_{21}, n_{22}.

Value

Normalizing constant is returned.

Author(s)

Jon Wakefield (jonno@u.washington.edu)

References

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

See Also

DirichNormHWE, DirichNormSat, DirichSampHWE, DirichSampSat, HWEDirichBF2, TriangNormHWE

Examples

nvec <- c(88,10,2)
TriangNormHWE(nvec)

[Package HWEBayes version 1.0 Index]