HWE.exact: Exact Test of Hardy-Weinberg Equilibrium for 2-Allele Markers

Description Usage Arguments Value Note Author(s) References See Also Examples

Description

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

Usage

1

Arguments

x

Genotype object

Value

Object of class 'htest'.

Note

This function only works for genotypes with exactly 2 alleles.

Author(s)

David Duffy davidD@qimr.edu.au with modifications by Gregory R. Warnes greg@warnes.net

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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)

kindlychung/genetics documentation built on May 20, 2019, 9:58 a.m.