HWE.exact {genetics}

Exact Test of Hardy-Weinberg Equilibrium for 2-Allele Markers

Description

Exact test of Hardy-Weinberg Equilibrium for 2 Allele Markers.

Usage

HWE.exact(x)

Arguments

x

Genotype object

Value

Object of class 'htest'.

Note

This function only works for genotypes with exactly 2 alleles.

Author(s)

David Duffy [e-mail] with modifications by Gregory R. Warnes [e-mail]

References

Emigh TH. (1980) "Comparison of tests for Hardy-Weinberg Equilibrium", Biometrics, 36, 627-642.

See Also

HWE.chisq, HWE.test, diseq, diseq.ci

Examples

example.data   <- c("D/D","D/I","D/D","I/I","D/D",
                    "D/D","D/D","D/D","I/I","")
g1  <- genotype(example.data)
g1

HWE.exact(g1)
# compare with
HWE.chisq(g1)




g2 <- genotype(sample( c("A","C"), 100, p=c(100,10), rep=TRUE),
               sample( c("A","C"), 100, p=c(100,10), rep=TRUE) )
HWE.exact(g2)