Testing Hardy‐Weinberg Equilibri

2021-01-19 本文已影响0人 MeUpUp

Testing Hardy‐Weinberg Equilibrium

reference：Allele Frequencies and Hardy‐Weinberg Equilibrium，Summer Institute in Statistical Genetics 2013，Module 8，Topic 2，66

Why should we test HWE locus?

When a locus is not in HWE, then this suggests one or more of the Hardy-Weinberg assumptions is false.
Departure from HWE has been used to infer the existence of natural selection, argue for the existence of assortative(non‐random) mating, and infer genotyping errors.

Two popular way to test HWE

Chi‐Square test

Exact test

Chi-square test

Testing procedures:

Compares observed genotype counts with the values expected under Hardy‐Weinberg for a locus with two alleles, in which $p=p(A)=(n_{Aa}+2n_{AA})/2n$

the expected is calculated from:

	'A'(p)	'a'(q)	TOTAL
'A'(p)	AA( $p^{2}$ )	Aa(pq)
'a'(q)	Aa(pq)	aa( $q^{2}$ )
			n

test statistic:

$X^{2}=\sum_{genotypes}\frac{(observed count-expected count)^2}{expected count}=n(\frac{4n_{AA}n_{aa}-n^{2}_{Aa}}{2n_{AA}+n_{Aa}(2n_{aa}+n_{Aa})})^{2}\sim \chi^{2}_{1}$

$df=(2-1)\times (2-1)=1$

When to use?

the expected count should be at least 5 in every cell.
If allele frequencies are low, and/or sample size is small, and/or there are many alleles at a locus, this may be a problem.

Exact test

The Hardy‐Weinberg exact test is based on calculating probabilities P(genotype counts | allele counts) under HWE.

$P(n_{Aa}|n_{A},n_{a},HWE)=\frac{n!}{n_{AA}!n_{Aa}!n_{aa}!}\frac{2^{n_{Aa}}n_{A}!n_{a}!}{(2n)!}$

To get the p‐value, sum the probabilities of all configurations with probability equal to or less that the observed configuration.

Compare two tests

The exact test is always conservative; the chi-square test can be either conservative or anticonservative.
Exact Test should be preferred: smaller sample sizes and/or multiallelic loci, since the chi-square test is prima facie not valid in these cases
(rule of thumb: must expect at least 5 in each cell)
The coarseness of Exact Test means it is conservative.