# HWEDirichBF2: Evaluates the Bayes factor in the k=2 allele case under... In HWEBayes: Bayesian investigation of Hardy-Weinberg Equilibrium via estimation and testing.

## Description

Function to evaluate the Bayes factor \Pr(n|\mbox{ HWE })/\Pr(n|\mbox{ saturated model }) in the k=2 allele case and with conjugate (Dirichlet) priors under HWE and saturated models.

## Usage

 `1` ```HWEDirichBF2(nvec, bvec0, bvec1) ```

## Arguments

 `nvec` vector of genotype frequencies in the order n_{11}, n_{21}, n_{22}. `bvec0` vector of length k=2 Dirichlet prior parameters for the prior under the null, where k is the number of alleles. `bvec1` vector of length k(k+1)/2=3 Dirichlet prior parameters for the prior under the saturated model, where k is the number of alleles.

## Value

Bayes factor is returned.

## Author(s)

Jon Wakefield ([email protected])

## References

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

`DirichNormHWE`, `DirichNormSat`, `DirichSampHWE`, `DirichSampSat`, `HWETriangBF2`, `TriangNormHWE`
 `1` ```HWEDirichBF2(nvec=c(88,10,2),bvec0=c(1,1),bvec1=c(1,1,1)) ```