Abstract
Objective: Exact tests for Hardy-Weinberg equilibrium are widely used in genetic association studies. We evaluate the mid p-value, unknown in the genetics literature, as an alternative for the standard p-value in the exact test.
Method: The type 1 error rate and the power of the exact test are calculated for different sample sizes, sigificance levels, minor allele counts and degrees of deviation from equilibrium. Three different p-value are considered: the standard two-sided p-value, the doubled one-sided p-value and the mid p-value. Practical implications of using the mid p-value are discussed with HapMap datasets and a data set on colon cancer.
Results: The mid p-value is shown to have a type 1 error rate that is always closer to the nominal level, and to have better power. Differences between the standard p-value and the mid p-value can be large for insignificant results, and are smaller for significant results. The analysis of empirical databases shows that the mid p-value uncovers more significant markers, and that the equilibrium null distribution is not tenable for both databases.
Conclusion: The standard exact p-value is overly conservative, in particular for small minor allele frequencies. The mid p-value ameliorates this problem by bringing the rejection rate closer to the nominal level, at the price of ocasionally exceeding the nominal level.
This study was supported by grants ECO2011-28875 and CODA-RSS MTM2009-13272 of the Spanish Ministry of Education and Science. The colorectal cancer study was supported by the Instituto de Salud Carlos III grants FIS PI08/1635, FIS PI08/1359 and CIBERESP CB07/02/2005. The authors thank the referees for their comments which have helped to improve the paper.
References
Agresti, A. (2002): Categorical data analysis. New York: John Wiley & Sons, second edition.10.1002/0471249688Search in Google Scholar
Attia, J., A. Thakkinstian, P. McElduff, E. Milne, S. Dawson, R. J. Scott, N. de Klerk, B. Armstrong and J. Thompson (2010): “Detecting genotyping error using measures of degree of Hardy-Weinberg disequilibrium,” Stat. Appl. Genet. Mol. Biol., 9, doi: 10.2202/1544-6115.1463.10.2202/1544-6115.1463Search in Google Scholar PubMed
Ayres, K. L. and D. J. Balding (1998): “Measuring departures from Hardy-Weinberg: a markov chain monte carlo method for estimating the inbreeding coefficient,” Heredity, 80, 769–777.10.1046/j.1365-2540.1998.00360.xSearch in Google Scholar PubMed
Barnard, G. A. (1989): “On alleged gains in power from lower p-values,” Stat. Med., 8, 1469–1477.Search in Google Scholar
Berry, G. and P. Armitage (1995): “Mid-p confidence intervals: a brief review,” J. Roy. Stat. Soc. D, 44, 417–423.Search in Google Scholar
Chakraborty, R. and Y. Zhong (1994): “Statistical power of an exact test of Hardy-Weinberg proportions of genotypic data at a multiallelic locus,” Hum. Hered., 44, 1–9.Search in Google Scholar
Elston, R. C. and R. Forthofer (1977): “Testing for Hardy-Weinberg equilibrium in small samples,” Biometrics, 33, 536–542.10.2307/2529370Search in Google Scholar
Emigh, T. H. (1980): “A comparison of tests for Hardy-Weinberg equilibrium,” Biometrics, 36, 627–642.10.2307/2556115Search in Google Scholar
Graffelman, J. (2010): “The number of markers in the hapmap project: some notes on chi-square and exact tests for Hardy-Weinberg equilibrium,” Am. J. Hum. Genet., 86, 813–818.10.1016/j.ajhg.2009.11.019Search in Google Scholar PubMed PubMed Central
Graffelman, J. and J. Morales-Camarena (2008): “Graphical tests for Hardy-Weinberg equilibrium based on the ternary plot,” Hum. Hered., 65, 77–84.Search in Google Scholar
Guo, W. S. and E. A. Thompson (1992): “Performing the exact test of Hardy-Weinberg proportion for multiple alleles,” Biometrics, 48, 361–372.10.2307/2532296Search in Google Scholar
Haldane, J. B. S. (1954): “An exact test for randomness of mating,” J. Genet., 52, 631–635.Search in Google Scholar
Hirji, K. F. (1991): “A comparison of exact, mid-p, and score tests for matched case-control studies,” Biometrics, 47, 487–496.10.2307/2532140Search in Google Scholar
Hosking, L., S. Lumsden, K. Lewis, A. Yeo, L. McCarthy, A. Bansal, J. Riley, I. Purvis, and C. Xu (2004): “Detection of genotyping errors by Hardy-Weinberg equilibrium testing,” Eur. J. Hum. Genet., 12, 395–399.Search in Google Scholar
Lancaster, H. O. (1961): “Significance tests in discrete distributions,” J. Am. Stat. Assoc., 56, 223–234.Search in Google Scholar
Landi, S., F. Gemignania, V. Moreno, L. Gioia-Patricola, A. Chabrier, E. Guino, M. Navarro, J. de Oca, F. Capella and F. Canzian (2005): “A comprehensive analysis of phase i and phase ii metabolism gene polymorphisms and risk of colorectal cancer,” Pharmacogenetics and Genomics, 15, 535–546.10.1097/01.fpc.0000165904.48994.3dSearch in Google Scholar PubMed
Lee, W. C. (2003): “Searching for disease-susceptibility loci by testing for Hardy-Weinberg disequilibrium in a gene bank of affected individuals,” Am. J. Epidemiol., 158, 397–400.Search in Google Scholar
Levene, H. (1949): “On a matching problem arising in genetics,” Ann. Math. Stat., 20, 91–94.Search in Google Scholar
Lindley, D. V. (1988): Statistical inference concerning Hardy-Weinberg equilibrium. In Bernardo, J. M., M. H. DeGroot, D. V. Lindley and A. F. M. Smith, (Eds.), Bayesian Statistics, 3, Oxford: Oxford University Press, pp. 307–326.Search in Google Scholar
Okamoto, M. and G. Ishii (1961): “Test of independence in intraclass 2×2 tables,” Biometrika, 48, 181–190.10.1093/biomet/48.1-2.181Search in Google Scholar
Rohlfs, R. V. and B. S. Weir (2008): “Distributions of Hardy-Weinberg equilibrium test statistics,” Genetics, 180, 1609–1616.10.1534/genetics.108.088005Search in Google Scholar PubMed PubMed Central
Salanti, G., G. Amountza, E. E. Ntzani and J. P. A. Ioannidis (2005): “Hardy-Weinberg equilibrium in genetic association studies: an empirical evaluation of reporting, deviations, and power,” Eur. J. Hum. Genet., 13, 840–848.Search in Google Scholar
Shoemaker, J., I. Painter and B. S. Weir (1998): “A bayesian characterization of Hardy-Weinberg disequilibrium,” Genetics, 149, 2079–2088.10.1093/genetics/149.4.2079Search in Google Scholar PubMed PubMed Central
Smith, C. A. B. (1986): “Chi-squared tests with small numbers,” Ann. Hum. Genet., 50, 163–167.Search in Google Scholar
The International HapMap Consortium (2003): “The international hapmap project,” Nature, 426, 789–796.10.1038/nature02168Search in Google Scholar PubMed
The International HapMap Consortium (2005): “A haplotype map of the human genome,” Nature, 437, 1299–1320.10.1038/nature04226Search in Google Scholar PubMed PubMed Central
The International HapMap Consortium (2007): “A second generation human haplotype map of over 3.1 million snps,” Nature, 449, 851–861.Search in Google Scholar
Wakefield, J. (2010): “Bayesian methods for examining Hardy-Weinberg equilibrium,” Biometrics, 66, 257–265.10.1111/j.1541-0420.2009.01267.xSearch in Google Scholar PubMed PubMed Central
Weir, B. S. (1996): Genetic Data Analysis II. Massachusetts: Sinauer Associates.Search in Google Scholar
Wigginton, J. E., D. J. Cutler and G. R. Abecasis (2005): “A note on exact tests of Hardy-Weinberg equilibrium,” Am. J. Hum. Genet., 76, 887–893.Search in Google Scholar
Yates, F. (1984): “Tests of significance for 2×2 contingency tables,” J. Roy. Stat. Soc A., 147, 426–463.Search in Google Scholar
Yu, C., S. Zhang, C. Zhou, and S. Sile (2009): “A likelihood ratio test of population Hardy-Weinberg equilibrium of case-control studies,” Genet. Epidemiol., 33, 275–280.Search in Google Scholar
©2013 by Walter de Gruyter Berlin Boston
