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.
Institute of Statistics and Decision Sciences, Duke University, Durham, North Carolina 27708-0251, USA. mhuber@math.duke.edu
Mark
M
Huber
Yuguo
Y
Chen
Ian
I
Dinwoodie
Adrian
A
Dobra
Mike
M
Nicholas
D000465:Algorithms; D001699:Biometry; D005828:Genetics, Population / Q000379:methods*; D009010:Monte Carlo Method; D013997:Time Factors