Monte Carlo algorithms for Hardy-Weinberg proportions.
The Hardy-Weinberg law is among the most important principles in the study of biological systems. Given its importance, many tests have been devised to determine whether a finite population follows Hardy-Weinberg proportions. Because asymptotic tests can fail, Guo and Thompson developed an exact test; unfortunately, the Monte Carlo method they proposed to evaluate their test has a running time that grows linearly in the size of the population N. Here, we propose a new algorithm whose expected running time is linear in the size of the table produced, and completely independent of N. In practice, this new algorithm can be considerably faster than the original method.
Huber
Mark
M
Institute of Statistics and Decision Sciences, Duke University, Durham, North Carolina 27708-0251, USA. mhuber@math.duke.edu
Chen
Yuguo
Y
Dinwoodie
Ian
I
Dobra
Adrian
A
Nicholas
Mike
M
D000465:Algorithms; D001699:Biometry; D005828:Genetics, Population / Q000379:methods*; D009010:Monte Carlo Method; D013997:Time Factors