Phase transition and regularized bootstrap in large-scale t-tests with false discovery rate control

From MaRDI portal
Publication:480982

DOI10.1214/14-AOS1249zbMATH Open1305.62213arXiv1310.4371OpenAlexW2963682607MaRDI QIDQ480982FDOQ480982

Q. M. Shao, Weidong Liu

Publication date: 12 December 2014

Published in: The Annals of Statistics (Search for Journal in Brave)

Abstract: Applying Benjamini and Hochberg (B-H) method to multiple Student's t tests is a popular technique in gene selection in microarray data analysis. Because of the non-normality of the population, the true p-values of the hypothesis tests are typically unknown. Hence, it is common to use the standard normal distribution N(0,1), Student's t distribution tn1 or the bootstrap method to estimate the p-values. In this paper, we first study N(0,1) and tn1 calibrations. We prove that, when the population has the finite 4-th moment and the dimension m and the sample size n satisfy logm=o(n1/3), B-H method controls the false discovery rate (FDR) at a given level alpha asymptotically with p-values estimated from N(0,1) or tn1 distribution. However, a phase transition phenomenon occurs when logmgeqc0n1/3. In this case, the FDR of B-H method may be larger than alpha or even tends to one. In contrast, the bootstrap calibration is accurate for logm=o(n1/2) as long as the underlying distribution has the sub-Gaussian tails. However, such light tailed condition can not be weakened in general. The simulation study shows that for the heavy tailed distributions, the bootstrap calibration is very conservative. In order to solve this problem, a regularized bootstrap correction is proposed and is shown to be robust to the tails of the distributions. The simulation study shows that the regularized bootstrap method performs better than the usual bootstrap method.


Full work available at URL: https://arxiv.org/abs/1310.4371




Recommendations





Cited In (28)

Uses Software





This page was built for publication: Phase transition and regularized bootstrap in large-scale \(t\)-tests with false discovery rate control

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q480982)