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

Share:

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

1

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

See Also

TriangNormHWE

Examples

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