HardyWeinberg (version 1.6.1)

dgraffelmanweir.bi: 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

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.

See Also

HWExact, HWExactStats

Examples

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