Volume 6 Supplement 1
Genetic Analysis Workshop 14: Microsatellite and singlenucleotide polymorphism
Power and type I error rate of false discovery rate approaches in genomewide association studies
 Qiong Yang^{1}Email author,
 Jing Cui^{2},
 Irmarie Chazaro^{3},
 L Adrienne Cupples^{1} and
 Serkalem Demissie^{1}
DOI: 10.1186/147121566S1S134
© Yang et al; licensee BioMed Central Ltd 2005
Published: 30 December 2005
Abstract
In genomewide genetic studies with a large number of markers, balancing the type I error rate and power is a challenging issue. Recently proposed false discovery rate (FDR) approaches are promising solutions to this problem. Using the 100 simulated datasets of a genomewide marker map spaced about 3 cM and phenotypes from the Genetic Analysis Workshop 14, we studied the type I error rate and power of Storey's FDR approach, and compared it to the traditional Bonferroni procedure. We confirmed that Storey's FDR approach had a strong control of FDR. We found that Storey's FDR approach only provided weak control of familywise error rate (FWER). For these simulated datasets, Storey's FDR approach only had slightly higher power than the Bonferroni procedure. In conclusion, Storey's FDR approach is more powerful than the Bonferroni procedure if strong control of FDR or weak control of FWER is desired. Storey's FDR approach has little power advantage over the Bonferroni procedure if there is low linkage disequilibrium among the markers. Further evaluation of the type I error rate and power of the FDR approaches for higher linkage disequilibrium and for haplotype analyses is warranted.
Background
Singlenucleotide polymorphisms (SNPs) are the most frequent types of polymorphisms and are commonly used in association mapping of candidate genomic regions. With the completion of the whole human genome sequence and the reduction of costs in SNP genotyping, genomewide studies with a dense SNP map consisting of 200,000 to 1 million SNP markers will become available in the near future [1]. How to efficiently control the false positive rate, or type I error rate, when a large number of tests are conducted in a genomewide study is a challenging problem.
For multiple testing problems, the chance to have a false significant test is higher than the nominal level even if each test is controlled at that nominal level. It is therefore of importance to control the probability of having one or more false significant tests. This probability is commonly referred as the familywise error rate (FWER). That can be written as P(V > 0), where V is the number of rejections of truly null hypotheses. There can be different types of controls for FWER: weak, exact, and strong, corresponding to conditioning on A = all null hypotheses are true (P(V > 0A)), B = the exact set of truly null hypotheses (P(V > 0B)), and C = any subset of null hypotheses are true (P(V > 0C)), respectively. While it is most desirable to have exact control of FWER, it is impossible to calculate a pvalue because the exact set of truly null hypotheses is unknown. Weak control and strong control are possible to evaluate, but the former may not be sufficient and the latter may be too conservative.
The simplest solution to the multiple testing problems is to apply the Bonferroni procedure, where each test is controlled at α/m, with m being the total number of tests. The Bonferroni procedure yields a strong control of FWER at α, but it could be too conservative when many tests are dependent. Resampling procedures, for example, minP [2], is less conservative than the Bonferroni procedure. The minP procedure generally provides weak control of FWER, and under certain conditions, also renders a strong control of FWER. However it is not always possible to examine the conditions and the resampling could be computationally demanding when there are a large number of tests.
Recently proposed false discovery rate (FDR) procedures by Benjamini and Hochberg [3] and by Storey [4], do not require resampling. Instead of controlling FWER, the FDR procedures control for a less restrictive measure, FDR, i.e., the expected proportion of truly null hypotheses among all the rejected null hypotheses. FDR can be written as E(V/R), where R is the total number of rejected null hypotheses. Both approaches give strong control of FDR and weak control of FWER. A distinction between Storey's and Benjamini's approaches is that the former incorporates an estimated proportion of truly null hypotheses while the latter assumes all hypotheses are truly null. As a result, the FDR approach of Storey is more powerful as a result [4].
There are few published simulation studies that compare FDR procedures with existing approaches for genomewide association studies. In this work, we will evaluate the power and type I error rate of FDR approach in the context of genomewide association tests with closely spaced markers using the Genetic Analysis Workshop 14 (GAW14) simulated datasets.
Methods
Subjects
We used all 100 simulated replicate datasets of the Danacaa population of GAW14. In each dataset, there were 100 nuclear families, each with an average of 6 members in which at least 2 were affected. There were a total of 1,469 genetic markers used in this study, including 416 microsatellite markers spaced approximately 7 cM apart, 917 SNPs spaced approximately 3 cM apart, and an additional 136 SNPs spaced approximately 0.3 cM apart in the region of disease loci D1, D2, D3, and D4.
In the Danacaa population, only disease loci D1, D2, and D5 (D5 was a marker) were related to the disease status.
Association analyses
A familybased association test (FBAT) was carried out at each genetic marker using the FBAT program http://www.biostat.harvard.edu/~fbat/default.html. In FBAT, a score statistic capturing the covariance between the phenotype and marker genotype is computed and standardized by its mean and variance conditional on the parental marker genotype data to adjust for population admixture. The score statistic is distributed as chisquare with degree of freedom equal to the number of different alleles at the marker minus one.
FDRadjusted pvalue: qvalue
To adjust for multiple testing, for each dataset all the pvalues from FBAT tests for each dataset were input into the QVALUE program http://faculty.washington.edu/~jstorey/qvalue/index.html to compute FDR adjusted pvalues, called qvalues. The qvalue of the j^{th} test with pvalue p_{ j }is defined as , and is the estimated proportion of truly null hypotheses; m is the total number of hypothesis tests. Because the pvalues of truly null hypotheses are uniformly distributed, one can imagine that most of the large pvalues, say, p > λ, are corresponding to truly nulls. Thus, π_{0} can be estimated as , where λ is chosen using a smoother method or a bootstrap method to balance the bias and variance of [4]. A null hypothesis will be rejected if the corresponding qvalue is less than or equal to α.
Estimation of FDR, FWER, and power
FDR, FWER, and power were estimated using the results of the 100 simulated datasets. To evaluate FWER when all the null hypotheses are true (i.e., m_{0} = m, where m_{0} = π_{0}·m is the number of truly null hypotheses), we used an arbitrary phenotype, that was a random 0–1 variable not associated with any marker. To evaluate the power and FWER when not all null hypotheses are true (i.e., m_{0} <m), we used disease status as the phenotype in the analyses.
Results
FDR of Storey's FDR approach
Estimated FDR for Storey's FDR approach. Estimated FDR when none of the SNPs associated with the phenotype (m_{0} = m) and when some of the SNPs associated with the phenotype (m_{0} <m) based on the 100 simulated datasets.
α  FDR^{a}  

m_{0} = m  m_{0} <m  
0.01  0.01  0.01 
0.05  0.08  0.09 
0.1  0.15  0.15 
0.15  0.20  0.19 
FWER of Bonferroni and Storey's FDR approaches
Estimated FWER of Bonferroni and Storey's FDR approaches. Estimated FWER when none of the SNPs associated with the phenotype (m_{0} = m) and when some of the SNPs associated with the phenotype (m_{0} <m) based on the 100 simulated datasets.
α  FWER (m_{0} = m)^{a}  FWER (m_{0} <m)  

Bonferroni  Storey's FDR  Bonferroni  Storey's FDR  
0.01  0  0  0.03  0.03 
0.05  0.01  0.01  0.08  0.21 
0.10  0.03  0.03  0.15  0.37 
0.15  0.06  0.06  0.21  0.46 
Power of Bonferroni and Storey's FDR approaches
Power of Bonferroni and Storey's FDR approaches. Two measures of power, E(T) and P(T > 0), were estimated. The E(T) was the expected number of detected true associations and P(T > 0) was the probability that at least one true associations was detected.
α  E(T)  P(T > 0)  

Bonferroni  Storey's FDR  Bonferroni  Storey's FDR  
0.01  1.0  1.0  0.76  0.78 
0.05  1.5  1.7  0.90  0.92 
0.10  1.6  2.0  0.93  0.93 
0.15  1.7  2.1  0.93  0.95 
Discussion
In this study, we have evaluated the FDR, FWER, and power of Storey's FDR approach, for a genomewide association study with a closely spaced marker map. We found FDR was slightly inflated for Storey's FDR approach based on 1,469 tests that are mostly independent using 100 simulated datasets. We noted that Storey's FDR approach only yielded weak control of FWER while the Bonferroni method provided close to strong control of FWER. For these simulated dataset, the power of Bonferroni and Storey's FDR approaches was comparable, although the latter yielded slightly higher power at the price of more false positives at markers that were not associated with the disease status.
Most markers in this study were spaced at approximately 3 cM apart and thus there was low LD among them. That may be one of the reasons that the power of Storey's FDR approach is not much higher than the Bonferroni method as expected. Another possible reason is that Storey's FDR approach was originally proposed in the context of microarray data analyses where many null hypotheses were not true. For current study, only three loci contributed to the disease, and only a few surrounding SNPs may be associated with the disease loci.
Unlike single testing problem in which type I error is uniquely defined as the probability that the pvalue ≤ α given that the null hypothesis is true, there is no unique definition for type I error rate for multiple testing problem. The FWER may be the most commonly used type I error measure for multiple testing problems. But even for FWER, there can be different types of control. The FDR is a valuable addition to the measures for type I error rate for multiple testing problems. Understanding which measure a multiple testing procedure can control is fundamental for appropriately using the procedures and interpreting the results.
For dependent tests, Storey's FDR approach is only valid under weak dependence and with a large number of tests. Weak dependence can be described as a form of dependence that almost does not affect the validity of the approach as the number of tests increases to infinity. An example of weak dependence is dependence in finite blocks. As a rule of thumb, several thousand tests are needed to obtain reasonable estimate of qvalues and thus control of FDR [5].
For tests that do not satisfy weak dependence, the properties of FDR approaches may be unclear. In this case, resampling methods such as the improved stepdown minP algorithm by Ge et al. [6] may be more promising to balance type I error and power. One such scenario is haplotype analyses in genomewide genetic studies, where SNPs are analyzed by group according to a moving window or LD blocks. Because the groups may have some SNPs in common, many of the tests can be dependent and the dependence among the tests can be very complicated.
Conclusion
In conclusion, Storey's FDR approach is more powerful than the Bonferroni procedure if strong control of FDR or weak control of FWER is desired. Storey's FDR approach has little power advantage over the Bonferroni procedure if there is low LD among the markers. Further evaluation of the type I error rate and power of the FDR approaches for higher LD and for haplotype analyses is warranted.
Abbreviations
 FBAT:

Familybased association test
 FDR:

False discovery rate
 FWER:

Familywise error rate
 GAW14:

Genetic Analysis Workshop 14
 LD:

Linkage disequilibrium
 SNP:

Singlenucleotide polymorphism
Declarations
Acknowledgements
This work is supported in part by contract N01HC25195.
Authors’ Affiliations
References
 Carlson CS, Eberle MA, Kruglyak L, Nickerson DA: Mapping complex disease loci in wholegenome association studies. Nature. 2004, 429: 446452. 10.1038/nature02623.View ArticlePubMedGoogle Scholar
 Westfall PH, Young SS: ResamplingBased Multiple Testing: Examples and Methods for PValue Adjustment. 1993, New York: WileyGoogle Scholar
 Benjamini Y, Hochberg Y: Controlling the false discovery rate – a practical and powerful approach to multiple testing. J Roy Stat Soc Ser B Meth. 1995, 57: 289300.Google Scholar
 Storey JD: A direct approach to false discovery rates. J Roy Stat Soc Ser B Meth. 2002, 64: 479498. 10.1111/14679868.00346.View ArticleGoogle Scholar
 Storey JD, Taylor JE, Siegmund D: Strong control, conservative point estimation and simultaneous conservative consistency of false discovery rates: a unified approach. J Roy Stat Soc Ser B Meth. 2004, 66: 187205. 10.1111/j.14679868.2004.00439.x.View ArticleGoogle Scholar
 Ge Y, Dudoit S, Speed TP: Resamplingbased multiple testing for microarray data analysis. Test. 2003, 12: 177. 10.1007/BF02595811.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.