Abstract:To the Editor: It is well known that Hardy-Weinberg equilibrium (HWE) is an important property in population genetics. Deviation from HWE among cases can provide evidence for a valid association.1Feder J.N. Gnirke A. Thomas W. Tsuchihashi Z. Ruddy D.A. Basava A. Dormishian F. Domingo Jr., R. Ellis M.C. Fullan A. et al.A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis.Nat. Genet. 1996; 13: 399-408Crossref PubMed Scopus (3278) Google Scholar, 2Nielsen D.M. Ehm M.G. Weir B.S. Detecting marker-disease association by testing for HW disequilibrium at a marker locus.Am. J. Hum. Genet. 1998; 63: 1531-1540Abstract Full Text Full Text PDF PubMed Scopus (317) Google Scholar, 3Czika W. Weir B.S. Properties of the multiallelic trend test.Biometrics. 2004; 60: 69-74Crossref PubMed Scopus (16) Google Scholar, 4Wittke-Thompson J.K. Pluzhnikov A. Cox N.J. Rational inferences about departures from HW equilibrium.Am. J. Hum. Genet. 2005; 76: 967-986Abstract Full Text Full Text PDF PubMed Scopus (324) Google Scholar Thus, it would be advisable to incorporate information from the HWE test for the improvement of power in detecting associated variants in genetic association studies. In the July 2008 issue of The Journal, Wang et al.5Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (31) Google Scholar described a test statistic, the tail-strength (TS) measure,6Taylor J. Tibshirani R. A tail strength measure for assessing the overall univariate significance in a dataset.Biostatistics. 2006; 7: 167-181Crossref PubMed Scopus (61) Google Scholar for evaluation of the global null hypothesis, that the SNP was not associated with disease, which is a function of two p values: one from a logistic-regression test in a genetic association study and one from a HWE test in cases. The authors further extended the mean-based TS measure to a median-based measure (TSM) by measuring the deviation of each p value from its median value instead of its expected value. On the basis of simulation studies and real disease studies, the authors stated that the adopted TS measure was more powerful than the traditional logistic-regression test and that the type I error was also well controlled. However, we have two main concerns about these conclusions. First, the two assumptions, which are required for deriving the exact distribution of TS and TSM statistics under the null hypothesis, hold only in certain scenarios. Violation of the two assumptions fails to obtain the exact distribution that the authors derived. The first assumption is that the two p values, of the HWE test and the logistic-regression test, are independent. This assumption may be violated, given that the two tests use the same case data. The authors mention that the adopted TS measure allows dependence between individual tests; yet, how to take the correlation into consideration was not discussed. The second assumption presumes the null distribution of HWE test p values to be uniform (0, 1), and this is breached when the exact test is applied for assessment of HWE.7Rohlfs R.V. Weir B.S. Distributions of Hardy-Weinberg Equilibrium Test Statistics.Genetics. 2008; 180: 1609-1616Crossref PubMed Scopus (78) Google Scholar To evaluate the validity of these two assumptions, we generated 10,000 data sets of cases and controls (500 and 500) at various minor-allele frequencies (MAFs) (0.05, 0.1, 0.2, 0.3, 0.4, and 0.5) under the null hypothesis of no association between the SNP and disease status. We first assessed the correlations between the two p values for the HWE exact test and the association test (likelihood-ratio test) under four different genetic models, including additive, genotypic, dominant, and recessive models, respectively (Table 1). “Genotypic model” means that the genotypes are coded as categorical variables and the major-allele homozygote is taken as a reference. We observed nonignorable correlations between the two p values when a genotypic effect was modeled, whereas the correlations were low when an additive effect was modeled at different MAFs. For dominant and recessive models, the correlations between the two p values were not stable but were, instead, dependent on the MAFs. We also found that the HWE exact test p values did not correspond to a uniform distribution when we used the same simulated data, as described above (Figure 1). At the MAF of 0.05, the p values of the HWE tests were skewed to the left and did not follow a uniform distribution, and the degree of skew gradually diminished as the MAF increased. It is well known that under the null hypothesis, the p value based on a continuous test statistic has a uniform distribution over the interval [0,1], regardless of the sample size of the experiment.8Hung H.M. O'Neill R.T. Bauer P. Köhne K. The behavior of the P-value when the alternative hypothesis is true.Biometrics. 1997; 53: 11-22Crossref PubMed Scopus (131) Google Scholar What we observed in Figure 1 is due to the discreteness of the test. In essence, the HWE exact test is based on a discrete hypergeometric distribution of the data under the null hypothesis of HWE. Given the sample size and the MAF, only a finite number of possible distinct genotype configurations exists, and therefore, a finite number of possible p values generates a coarse distribution. When the sample size or MAF is small, the number of possible distinct genotype configurations will be small, with specific observed probabilities that deviate the distribution of p values from uniform. In particular, a spike close to a p value of 1 will appear for the most frequent sample configuration. As the sample size or MAF increases, more distinct genotype configurations increasingly resemble a uniform distribution.Table 1The p Value Correlation Coefficient between the HWE Exact Test and the Likelihood-Ratio Association TestGenetic ModelMinor-Allele Frequency0.050.10.20.30.40.5Additive−0.029−0.00630.00140.0140.012−0.0004GenotypicaThe genotypic model means that the genotypes are coded as categorical variables and the major-allele homozygote is taken as a reference.0.140.300.270.270.280.27Dominant−0.00370.0240.0410.0630.100.12Recessive0.160.290.220.170.140.095The results are based on 10,000 simulated data sets under the null hypothesis.a The genotypic model means that the genotypes are coded as categorical variables and the major-allele homozygote is taken as a reference. Open table in a new tab The results are based on 10,000 simulated data sets under the null hypothesis. To further evaluate the influence of these two assumptions on the analytical distribution of the TS and TSM statistics that the authors derived in the paper, we compared the analytical distribution of TS and TSM with the empirical distribution under different genetic models and MAF combinations, using the same simulated data sets. Results are shown in Figure S1 (available online). Specifically, for the additive model, when the MAF equals 0.05 or 0.1, the analytical p value is larger than the empirical p value at the right-hand tail of the distribution, which is more evident at a MAF of 0.05. However, the analytical and empirical distributions match pretty well at MAFs ≥ 0.2. Notably, when a genotypic model is used for calculation of the association test p value, the corresponding analytical and empirical distributions of the TS and TSM statistics do not fit well, regardless of the MAF, which can be attributed mainly to the nonignorable correlation between the two tests at various MAFs (violation of the independence assumption of the two tests). A similar result exists for a dominant or recessive genetic model when either the correlation between the two individual tests is nonignorable or the MAF is low (an MAF = 0.05 or 0.1). As an example, Figure 2 shows the empirical and the analytical distribution of the TS and TSM statistics at a MAF of 0.2 under the null hypothesis. We found that the analytical and the empirical distributions fit well for the additive model. However, it is obvious that the empirical and analytical distributions mismatch for a genotypic-effect model. The exact p value is less than empirical p value at the tails of the distribution for both TS and TSM. Therefore, we suggest the use of empirical distributions rather than exact distributions for TS and TSM measures in the practice of genetic association studies. Second, we are concerned with the power of the TS and TSM statistics. To evaluate the empirical power, we performed another 10,000 simulations using the same parameters implemented in Wang et al.; i.e., simulation of data 1 in model 1 with β0=−2 and β1=0.3. We calculated the power of the two marginal tests. Under the 0.05 significance level, the power for the HWE test is 0.053, which is only slightly larger than the type I error, and the power of additive effects in logistic regression is 0.59. Surprisingly, the empirical power of additive effects for the TS and TSM statistics, which combine the p values of the HWE exact test and the association test, is only 0.25 (0.27), which is lower than what the authors reported in their Table 4. We are unable to explain this discrepancy. In summary, the TS statistic is useful for combining the information from HWE test and case-control association test to improve the power of detecting SNP effects in genetic association studies. However, one needs to be cautious when using this statistic. On the basis of simulation results, we found that the analytical distributions of the TS and TSM statistics are influenced both by the MAF and by genetic models used in association tests. We suggest using the empirical p value, rather than the exact p value, in real situations. A more generalized statistic that does not depend on HWE-test significance in cases should be developed for the incorporation of HWE information and improvement of the power of genetic association studies. This work was supported by the U.S. National Institutes of Health (career development award K01 DA024758 to B.Z.Y.). We are grateful for Joel Gelernter's valuable discussion and comments on this letter. Download .pdf (.35 MB) Help with pdf files Document S1. One Figure Tail Strength to Combine Two p Values: Their Correlation Cannot Be IgnoredZang et al.The American Journal of Human GeneticsFebruary 13, 2009In BriefTo the Editor: The population-based case-control study is a useful approach to evaluating genetic association with many common and complex diseases. In general, one first uses the generalized linear model to fit the data and then uses an asymptotic test to detect the true association. In addition to this regression-based analysis, when Hardy-Weinberg equilibrium (HWE) holds in the population, testing HWE in cases has been used for indicating the association. Because the regression-based analyses (including the trend test and the likelihood-ratio test) are generally more powerful than testing HWE in cases, they are often employed in case-control studies. Full-Text PDF Open ArchiveResponse to Zang et al. and Han et al.Wang et al.The American Journal of Human GeneticsFebruary 13, 2009In BriefTo the Editor: In July 2008, we proposed a powerful test for the study of genetic association that incorporates information about deviation from Hardy-Weinberg proportions (HWP) in cases.1 Two approaches were proposed: the mean-based tail-strength (TS) measure and the median-based tail-strength (TSM) measure. These measures combined p values from the likelihood ratio test (LRT) for association and the exact test for HWP. For both measures, we derived exact formulas to compute p values, and we also provided an approach for obtaining empirical p values with the use of a resampling procedure. Full-Text PDF Open Archive

Abstract:To the Editor: The population-based case-control study is a useful approach to evaluating genetic association with many common and complex diseases. In general, one first uses the generalized linear model to fit the data and then uses an asymptotic test to detect the true association. In addition to this regression-based analysis, when Hardy-Weinberg equilibrium (HWE) holds in the population, testing HWE in cases has been used for indicating the association. Because the regression-based analyses (including the trend test and the likelihood-ratio test) are generally more powerful than testing HWE in cases, they are often employed in case-control studies. Less attention is paid to testing HWE in cases. In the July 2008 issue of The American Journal of Human Genetics, Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar proposed a novel approach of using the tail strength to combine the p value of the likelihood-ratio test (LRT) for association and the p value of an exact test for the deviation from HWE in cases. Taylor and Tibshirani2Taylor J. Tibshirani R. A tail strength measure for assessing the overall univariate significance in a dataset.Biostatistics. 2006; 7: 167-181Crossref PubMed Scopus (62) Google Scholar originally proposed the tail strength as a measure of the overall strength of association for a large number of hypotheses in microarray analyses and genome-wide association studies (GWAS). Compared to Fisher's combination of p values3Elston R.C. On Fisher method of combining p-values.Biometrical J. 1991; 33: 339-345Crossref Scopus (49) Google Scholar, which weights each p value equally, the tail strength weights each ordered p value by its expectation under the null hypothesis. The tail strength can be used for combining independent and dependent p values and is not restricted to any special genetic model underlying the data. Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar combined the two p values by using the tail strength and extended the original tail strength by using the medians of the ordered p values as weights. They derived asymptotic null distributions for the tail strengths by applying the additive model and using the mean and median as weights, respectively. Their results showed significant improvement in terms of the power when the tail strengths were used. They also showed that the type I errors were under control, although we notice that almost all reported type I errors in their tables are less than the nominal levels. Normally, when the tail strength is used as a test statistic, its asymptotic null distribution is approximated by Monte-Carlo simulation procedures. Simulation-based approaches to determining the tail probabilities or p values of complex statistics have limitations for applications in GWAS.4Sladek R. Rocheleau G. Rung J. Dina C. Shen L. Serre D. Boutin P. Vincent D. Belisle A. Hadjadj S. et al.A genome-wide association study identifies novel risk loci for type 2 diabetes.Nature. 2007; 445: 881-885Crossref PubMed Scopus (2366) Google Scholar, 5Conneely K.N. Boehnke M. So many correlated tests, so little time! Rapid adjustment of p values for multiple correlated tests.Am. J. Hum. Genet. 2007; 81: 1158-1168Abstract Full Text Full Text PDF PubMed Scopus (339) Google Scholar In this situation, deriving their asymptotic distributions is important. Although Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar derived the asymptotic null distributions and critical values for their tail-strength statistics, they assumed in their derivations that the two p values were independent even though in the introduction section they mentioned that they would use the tail strength to combine two dependent p values. When the two p values are correlated, their asymptotic null distributions may be inappropriate. Using two test statistics different from those in Wang and Shete,1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar Zheng and Ng6Zheng G. Ng H.K.T. Genetic model selection in two-phase analysis for case-control association studies.Biostatistics. 2008; 9: 391-399Crossref PubMed Scopus (69) Google Scholar noticed that the correlation between the p values of the trend test and testing HWE between cases and controls (HWDTT7Song K. Elston R.C. A powerful method of combining measures of association and Hardy-Weinberg disequilibrium for fine-mapping in case-control studies.Stat. Med. 2006; 25: 105-126Crossref PubMed Scopus (84) Google Scholar) could also vary from the recessive (REC) model to the additive (ADD) model, the multiplicative (MUL) model, and the dominant (DOM) model. As we mentioned before, Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar considered the tail strengths based only on the ADD model. However, the performance of testing HWE in cases would also vary across the genetic models. For example, it is known that testing HWE cannot detect association under the MUL model even though testing HWE has been used for detecting association.8Nielsen D.M. Ehm M.G. Weir B.S. Detecting marker-disease association by testing for Hardy-Weinberg disequilibrium at a marker locus.Am. J. Hum. Genet. 1998; 63: 1531-1540Abstract Full Text Full Text PDF PubMed Scopus (321) Google Scholar, 9Wittke-Thompson J.K. Pluzhnikov A. Cox N.J. Rational inferences about departure from Hardy-Weinberg equilibrium.Am. J. Hum. Genet. 2005; 76: 967-986Abstract Full Text Full Text PDF PubMed Scopus (331) Google Scholar, 10Wang T. Zhu X. Elston R.C. Improving power in contrasting linkage-disequilibrium patterns between cases and controls.Am. J. Hum. Genet. 2007; 80: 911-920Abstract Full Text Full Text PDF PubMed Scopus (30) Google Scholar Therefore, the performance of the tail strength of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar can be potentially affected by two factors that were either ignored or not examined in their article. One is the correlation between the two p values of the LRT and the test for HWE in cases, and the other is the unknown underlying genetic models. In this letter, using Monte-Carlo simulation procedures, we study the correlations between the p values of the LRT and the exact test for Hardy-Weinberg proportion in cases under the four genetic models. If the two p values are indeed correlated, we examine the performance of the tail-strength statistics of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar under the null and alternative hypotheses. The analytical formula of the correlation, if any, between the LRT and the exact test for HWE used in Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar is difficult to obtain. Therefore, we consider the combination of the p values of the trend test and chi-square test for HWE between cases and controls (HWDTT), from which the asymptotic correlation between the two p values has been obtained.6Zheng G. Ng H.K.T. Genetic model selection in two-phase analysis for case-control association studies.Biostatistics. 2008; 9: 391-399Crossref PubMed Scopus (69) Google Scholar, 7Song K. Elston R.C. A powerful method of combining measures of association and Hardy-Weinberg disequilibrium for fine-mapping in case-control studies.Stat. Med. 2006; 25: 105-126Crossref PubMed Scopus (84) Google Scholar This new tail strength with the correlation is denoted by TSC. We further derive its asymptotic null distribution and critical value (see Appendix A). Comparison between our TSC and that of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar is obtained by Monte-Carlo simulations under the null and alternative hypotheses. We also denote the tail strengths based on the mean and median in Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar by TS and TSM, respectively. Here we report the main results from our simulation study. In the simulation, we assumed HWE holds in the population. In each replicate, 500 cases and 500 controls were generated under the null hypothesis with the baseline penetrance fixed at 0.02 (the probability of disease with a genotype of zero risk alleles), and minor-allele frequency (MAF) increases from 0.1 to 0.5 in increments of 0.1. We used a total of 10,000 replicates to estimate the null correlations between the two p values, the type I error rates, and power. The nominal levels 0.01 and 0.05 were used. For LRT statistics, we considered 1-degree-of-freedom tests. Therefore, for each genetic model under the alternative hypothesis (REC, ADD/MUL, and DOM), an optimal test is available for the LRT or trend test. In the simulation, we consider three LRTs and three trend tests, optimal for the three genetic models. Therefore, a total of nine tail strengths were considered in the simulation: TS, TSM, and TSC each have three model choices depending on the targeted genetic model. The results of the null correlations between the two p values in Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar and corresponding type I errors are reported in Table 1 for the nominal level 0.01 and in Table 2 for the nominal level 0.05.Table 1Simulated Null Correlations of the Two p Values of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar and the Asymptotic Type I Errors with Nominal Level 0.01MAFModelSimulated Null CorrelationsTSTSMTSC0.1REC0.27020.02620.02720.0067ADD−0.00490.00560.00570.0081DOM0.02380.00720.0070.0090.2REC0.23270.02480.02510.0123ADD0.00180.00920.00920.0108DOM0.03280.01290.01280.01160.3REC0.16720.02680.0260.0131ADD0.01870.01040.01010.0112DOM0.05050.0170.0170.00920.4REC0.14540.02250.02250.0118ADD−0.01490.00760.00740.0083DOM0.07160.01570.01530.00920.5REC0.10370.01970.02010.0128ADD−0.00470.00740.00770.0103DOM0.09190.01740.01750.0081TS uses means as weights, TSM uses medians as weights, and TSC is the proposed test with the correlations. Three genetic models, which are only used for constructing the optimal LRTs (for TS and TSM) and optimal-trend tests (for TSC), are considered. Open table in a new tab Table 2Simulated Null Correlations of the Two p Values of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar and the Asymptotic Type I Errors with Nominal Level 0.05MAFModelSimulated Null CorrelationsTSTSMTSC0.1REC0.27020.08410.08320.0436ADD−0.00490.03940.04030.0477DOM0.02380.04090.04080.04620.2REC0.23270.07650.07720.0476ADD0.00180.04430.04410.0557DOM0.03280.05130.05240.04720.3REC0.16720.07270.0720.0482ADD0.01870.04380.04380.0468DOM0.05050.05190.05240.05310.4REC0.14540.06950.06780.0458ADD−0.01490.05190.05160.0521DOM0.07160.05740.05690.04940.5REC0.10370.07220.07040.0481ADD−0.00470.04740.04830.0509DOM0.09190.06470.06330.0475TS uses means as weights, TSM uses medians as weights, and TSC is the proposed test with the correlations. Three genetic models, which are only used for constructing the optimal LRTs (for TS and TSM) and optimal-trend tests (for TSC), are considered. Open table in a new tab TS uses means as weights, TSM uses medians as weights, and TSC is the proposed test with the correlations. Three genetic models, which are only used for constructing the optimal LRTs (for TS and TSM) and optimal-trend tests (for TSC), are considered. TS uses means as weights, TSM uses medians as weights, and TSC is the proposed test with the correlations. Three genetic models, which are only used for constructing the optimal LRTs (for TS and TSM) and optimal-trend tests (for TSC), are considered. The results in Table 1, Table 2 follow similar patterns. Thus, we focus on Table 1. The simulated null correlations between the p value of LRT and the p value of the exact HWE test in cases indicate that the null correlations are not zero when the LRT is optimal for the REC or DOM models, but they are close to zero for the ADD (MUL) model. Hence, the type I errors of the TS and TSM of Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar are under control when the LRT is optimal for the ADD (MUL) model but are largely inflated when the LRTs are optimal for the REC and DOM models, especially for the REC model. Note that Wang and Shete1Wang J. Shete S. A test for genetic association that incorporates information about deviation from Hardy-Weinberg proportions in cases.Am. J. Hum. Genet. 2008; 83: 53-63Abstract Full Text Full Text PDF PubMed Scopus (32) Google Scholar only considered the LRT optimal for the ADD model. Therefore, their type I errors were under control. On the other hand, the type I errors of TSC, which takes care of the correlations, are close to the nominal level regardless of the targeted genetic models. We also conducted simulations to compare the powers of the TS, TSM, and TSC. For the TS and TSM, the correlations between the two p values were not incorporated. Thus, on the basis of results in Table 1, Table 2, their powers could be inflated under the REC and DOM models, but not under the ADD and MUL models. The powers are presented for the TS, TSM, and TSC (from left to right) under the REC model (Figure 1) and ADD model (Figure 2). The plots for the MUL and DOM models can be found in the Supplemental Data available online (Figures S1 and S2, respectively). The parameter values of the simulations under the alternative hypotheses are similar to those in Table 1, Table 2, except that the genotype relative risk (gamma2, which is defined as the ratio of penetrances with two risk alleles to those with zero risk alleles) ranges from 1 to 2, and the MAF is fixed at 0.3. The “asymptotic” and “simulated” powers in the figures were based on the critical values obtained from 10,000 parametric bootstrap samples and 10,000 permutations, respectively.Figure 2The Asymptotic and Simulated Powers under the ADD ModelShow full captionThe tests from left to right are TS, TSM, and TSC. Gamma2 is the ratio of penetrances with two risk alleles to zero risk allele.View Large Image Figure ViewerDownload Hi-res image Download (PPT) The tests from left to right are TS, TSM, and TSC. Gamma2 is the ratio of penetrances with two risk alleles to zero risk allele. Figure 1 (under the REC model) shows that TS and TSM have similar powers and are more powerful than TSC. This could be due to the fact that TS and TSM had inflated type I errors as shown in Table 1, Table 2. On the other hand, Figure 2 shows that the powers of TS, TSM, and TSC are similar under the ADD model because the three statistics had similar type I errors. For the TS and TSM, the bootstrap and permutation procedures yield similar powers under the ADD, MUL, and DOM models but have slightly different powers under the REC model. We also studied empirical powers of the TSC, the optimal trend test, a robust test MAX311Freidlin B. Zheng G. Li Z. Gastwirth J.L. Trend tests for case-control studies of genetic markers: power, sample size and robustness.Hum. Hered. 2002; 53: 146-152Crossref PubMed Scopus (264) Google Scholar, and classical Pearson's test for association under the four genetic models. The description and summary of our findings are given in Appendix B. The results show that the TSC has moderate power improvement under the REC model but loses significant power under the ADD and MUL models. This can be explained by the fact that testing HWE has little power under the ADD and MUL models. In summary, the tail strength may improve power under some specific genetic models after correction for the correlation. However, when the underlying genetic model is unknown, the robust statistics are more preferable.6Zheng G. Ng H.K.T. Genetic model selection in two-phase analysis for case-control association studies.Biostatistics. 2008; 9: 391-399Crossref PubMed Scopus (69) Google Scholar, 11Freidlin B. Zheng G. Li Z. Gastwirth J.L. Trend tests for case-control studies of genetic markers: power, sample size and robustness.Hum. Hered. 2002; 53: 146-152Crossref PubMed Scopus (264) Google Scholar We would like to thank Yaning Yang for some helpful discussions that brought our attention to the tail strength. The work of Y. Zang and W.K. Fung were partially supported by The Croucher Foundation and China Natural Science Foundation (no. 10701067). Denote the HWDTT by Z∗, which is a statistic testing HWE between cases and control and was proposed by Song and Elston.7Song K. Elston R.C. A powerful method of combining measures of association and Hardy-Weinberg disequilibrium for fine-mapping in case-control studies.Stat. Med. 2006; 25: 105-126Crossref PubMed Scopus (84) Google Scholar Denote the trend test as Zx, where x = 0, 0.5, and 1 for the REC, ADD (MUL), and DOM models, respectively.11Freidlin B. Zheng G. Li Z. Gastwirth J.L. Trend tests for case-control studies of genetic markers: power, sample size and robustness.Hum. Hered. 2002; 53: 146-152Crossref PubMed Scopus (264) Google Scholar, 12Sasieni P.D. From genotype to genes: doubling the sample size.Biometrics. 1997; 53: 1253-1261Crossref PubMed Scopus (751) Google Scholar, 13Zheng G. Freidlin B. Li Z. Gastwirth J.L. Choice of scores in trend tests for case-control studies of candidate-gene associations.Biometrical J. 2003; 45: 335-348Crossref Scopus (59) Google Scholar Under the null hypothesis H0, (Z∗, Zx) follows the bivariate normal distribution N(0, Σ1) with the density function f1, where Σ1=(1ρxρx1), and (−Z∗, Zx) follows the bivariate distribution N(0, Σ2) with the density function f2, where Σ2=(1−ρx−ρx1). The expressions for ρx were given in Zheng and Ng for different x values.6Zheng G. Ng H.K.T. Genetic model selection in two-phase analysis for case-control association studies.Biostatistics. 2008; 9: 391-399Crossref PubMed Scopus (69) Google Scholar The following derivation can be modified to the tail strength of any two correlated p values. The p value of Z∗ is P∗ = 2Φ(− |z∗|), and the p value of Zx is Px = 2Φ(− |zx|), where Φ is the cumulative distribution function of the standard normal N(0, 1), and z∗ and zx are observed statistics. Then the joint distribution of P∗ and Px is:F(x1,x2)=Pr(P∗Φ−1(1−x12),|Zx|>Φ−1(1−x22))=Pr(Z∗<Φ−1(x12),Zx<Φ−1(x22))+Pr(Z∗<Φ−1(x12),−Zx<Φ−1(x22))+Pr(−Z∗<Φ−1(x12),Zx<Φ−1(x22))+Pr(−Z∗<Φ−1(x12),−Zx<Φ−1(x22))=2∫−∞Φ−1(x12)∫−∞Φ−1(x22)f1(y1,y2)dy1dy2+2∫−∞Φ−1(x12)∫−∞Φ−1(x22)f2(y1,y2)dy1dy2. Thus, its density function can be written asf(x1,x2)=∂F(x1,x2)∂x1∂x2=∑i=01exp[−{Φ−1(x12)}2+{Φ−1(x22)}2+(−1)i2ρx{Φ−1(x12)}{Φ−1(x22)}2(1−ρx2)]×{21−ρx2exp[−{Φ−1(x12)}2+{Φ−1(x22)}22]}−1. Therefore, the ordered p values have the cumulative function given byF(x(1),x(2))=Pr(P(1)

Differentiating the Cochran‐Armitage Trend Test and Pearson's Χ2 Test: Location and Dispersion

Decomposing Pearson's Χ2 Test: A Linear Regression and Its Departure from Linearity

A Nonparametric Alternative to the Cochran-Armitage Trend Test in Genetic Case-Control Association Studies: the Jonckheere-Terpstra Trend Test

Comparison of Two‐phase Analyses for Case–control Genetic Association Studies

Exploring Case-Control Genetic Association Tests Using Phase Diagrams.

Pearson's Test, Trend Test, and MAX Are All Trend Tests with Different Types of Scores

Single-Marker Analysis For Unmatched Case-Control Data

A Shrinkage Method for Testing the Hardy–Weinberg Equilibrium in Case‐Control Studies

Robust Genomic Control and Robust Delta Centralization Tests for Case-Control Association Studies

Power Analysis of Principal Components Regression in Genetic Association Studies.

Is the Tail-Strength Measure More Powerful in Tests of Genetic Association?

A Model-Embedded Trend Test with Incorporating Hardy-Weinberg Equilibrium Information

Single-Locus Genetic Association Analysis By Ordinal Tests

Tail Strength to Combine Two P Values: Their Correlation Cannot Be Ignored

Powerful Test of Heterogeneity in Two‐Sample Summary‐Data Mendelian Randomization

Pearson Chi-squared Conditional Randomization Test

Retrospective Versus Prospective Score Tests For Genetic Association With Case-Control Data

A consideration of the chi-square test of Hardy--Weinberg equilibrium in a non-multinomial situation

Enhancing the Power to Detect Low-Frequency Variants in Genome-Wide Screens

A Robust Test for Multi-Ordered 2 × $J$ Ordinal Contingency Tables

Population Admixture: Detection by Hardy-Weinberg Test and Its Quantitative Effects on Linkage-Disequilibrium Methods for Localizing Genes Underlying Complex Traits.