dgraffelmanweir.bi: Calculate Graffelman-Weir exact density for bi-allelic...

View source: R/dgraffelmanweir.bi.R

dgraffelmanweir.biR Documentation

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 jan.graffelman@upc.edu

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. \Sexpr[results=rd]{tools:::Rd_expr_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)

HardyWeinberg documentation built on May 29, 2024, 6:17 a.m.