HWETriangBF2 {HWEBayes} | R Documentation |
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.
HWETriangBF2(nvec)
nvec |
vector of genotype frequencies in the order n_{11}, n_{21}, n_{22}. |
Bayes factor is returned.
Jon Wakefield (jonno@u.washington.edu)
Wakefield, J. (2009). Bayesian methods for examining Hardy-Weinberg equilibrium. Biometrics.
TriangNormHWE
HWETriangBF2(nvec=c(88,10,2))