dgraffelmanweir.bi {HardyWeinberg}

Calculate Graffelman-Weir exact density for bi-allelic X-chromosomal variant

Description

Function dgraffelmanweir calculate the probability P(NAB=nab and MA=ma|NA=na) for a bi-allelic X-chromosomal variant.

Usage

dgraffelmanweir.bi(x, y)

Arguments

x

vector with male genotype counts (A,B)

y

vector with female genotype counts (AA,AB,BB)

Value

a single real number

Author(s)

Jan Graffelman [e-mail]

References

Graffelman, J. and Weir, B.S. (2016) Testing for Hardy-Weinberg equilibrium at bi-allelic genetic markers on the X chromosome. Heredity 116(6) pp. 558–568. doi:10.1038/hdy.2016.20

See Also

HWExact, HWExactStats

Examples

males <- c(A=392, B=212)
females <- c(AA=275, AB=296, BB=80)
prob <- dgraffelmanweir.bi(males,females)
print(prob)