Exact test of Hardy–Weinberg.

Share:

Description

Computes a p-value for the exact test of Hardy–Weinberg proportions, for genotype counts for a biallelic locus, as described by Wiggington et al. (2005)

Usage

1
snphwe(g)

Arguments

g

Vector of genotypes.

Details

This function is a wrapper for snphweCounts, which uses code by Wiggington et al. to compute a P-value for an exact test of Hardy–Weinberg proportions.

This function requires that when g is sorted, the order is hom/het/hom. For example g can be additively coded genotypes (0, 1, 2), or text coded genotypes like AA, AC, CC. For text coded genotypes the heterozygote MUST ALWAYS be coded the same way, e.g. a mixture of AC and CA are not allowed.

Value

The p-value.

Author(s)

Toby Johnson Toby.x.Johnson@gsk.com

Examples

1
snphwe(rbinom(100, 2, 0.2))

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.