Be the first to rate this file! 8 Downloads (last 30 days) File Size: 3.74 KB File ID: #27423

Two-stage Benjamini, Krieger, & Yekutieli FDR procedure

by

 

30 Apr 2010 (Updated )

Two-stage procedure for controlling the false discovery rate in a family of hypothesis tests

| Watch this File

File Information
Description

     Executes the "two-stage" Benjamini, Krieger, & Yekutieli (2006) procedure for controlling the false discovery rate (FDR) of a family of hypothesis tests. FDR is the expected proportion of rejected hypotheses that are mistakenly rejected (i.e., the null hypothesis is actually true for those tests). FDR is a somewhat less conservative/more powerful method for correcting for multiple comparisons than procedures like Bonferroni correction that provide strong control of the family-wise error rate (i.e., the probability that one or more null hypotheses are mistakenly rejected).
     The procedure implemented by this function is more powerful than the original Benjamini & Hochberg (1995) procedure when a considerable percentage of the hypotheses in the family are false. It is only slightly less powerful than the original procedure when there are very few false hypotheses. To the best of my knowledge, this procedure is only guaranteed to control FDR if the tests are independent. However, simulations suggest that it can control FDR even when the tests are positively correlated (Benjamini et al., 2006).

References:
  Benjamini, Y., Krieger, A.M., & Yekutieli, D. (2006) Adaptive linear step-up procedures that control the false discovery rate. Biometrika. 93(3), 491-507.

  Benjamini, Y. & Hochberg, Y. (1995) Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, Series B (Methodological). 57(1), 289-300.

MATLAB release MATLAB 7.10 (R2010a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
03 May 2010 Rob Campbell

Perhaps these files should be grouped together into the same zip file or even integrated into the same m file? At least you could provide linkso that people would know about the existence of both.

Updates
17 Sep 2012

Comments updated

Contact us