HWETriangBF2 {HWEBayes}R Documentation

Evaluates the Bayes factor in the k=2 allele case with a "triangular" prior under the null

Description

Function to evaluate the Bayes factor Pr(n|mbox{ HWE })/Pr(n|mbox{ saturated model }) in the k=2 allele case and with a conjugate (Dirichlet) priors under the saturated model and a "triangular" distribution under the null. the latter is the marginal prior distribution under the (1,1,1) Dirichlet prior under the saturated model.

Usage

HWETriangBF2(nvec)

Arguments

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

Value

Bayes factor is returned.

Author(s)

Jon Wakefield (jonno@u.washington.edu)

References

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

See Also

TriangNormHWE

Examples

HWETriangBF2(nvec=c(88,10,2))

[Package HWEBayes version 1.0 Index]