Estimation of the proportion of true null hypotheses under sparse dependence: Adaptive FDR controlling in microarray data

Abstract

The proportion of non-differentially expressed genes is an important quantity in microarray data analysis and an appropriate estimate of the same is used to construct adaptive multiple testing procedures. Most of the estimators for the proportion of true null hypotheses based on the thresholding, maximum likelihood and density estimation approaches assume independence among the gene expressions. Usually, sparse dependence structure is natural in modelling associations in microarray gene expression data and hence it is necessary to develop methods for accommodating the sparse dependence well within the framework of existing estimators. We propose a clustering based method to put genes in the same group that are not coexpressed using the estimated high dimensional correlation structure under sparse assumption as dissimilarity matrix. This novel method is applied to three existing estimators for the proportion of true null hypotheses. Extensive simulation study shows that the proposed method improves an existing estimator by making it less conservative and the corresponding adaptive Benjamini-Hochberg algorithm more powerful. The proposed method is applied to a microarray gene expression dataset of colorectal cancer patients and the results show gain in terms of number of differentially expressed genes. The R code is available at https://github.com/aniketstat/Proportiontion-of-true-null-under-sparse-dependence-2021.

Keywords

False discovery rate differentially expressed genes high dimensional data ANOVA Tocher’s method 62P10.

Get full access to this article

View all access options for this article.

References

Benjamini

Hochberg

. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological) 1995; 57: 289–300.

Benjamini

Krieger

Yekutieli

. Adaptive linear step-up procedures that control the false discovery rate. Biometrika 2006; 93: 491–507.

Lim

Tan

Lim

, et al. Compendiums of cancer transcriptomes for machine learning applications. Sci Data 2019a; 6: 1–8.

Lim

Chua

MLK

Yeong

JPS

, et al. Pan-cancer analysis connects tumor matrisome to immune response. NPJ Precision Oncology 2019b; 3: 1–9.

Friguet

Kloareg

Causeur

. A factor model approach to multiple testing under dependence. J Am Stat Assoc 2009; 104: 1406–1415.

Sun

Tony Cai

. Large-scale multiple testing under dependence. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2009; 71: 393–424.

. On false discovery control under dependence. Ann Stat 2008; 36: 364–380.

Friguet

Causeur

. Estimation of the proportion of true null hypotheses in high-dimensional data under dependence. Comput Stat Data Anal 2011; 55: 2665–2676.

Ostrovnaya

Nicolae

. Estimating the proportion of true null hypotheses under dependence. Stat Sin 2012; 22: 1689–1716.

10.

Dobra

Hans

Jones

, et al. Sparse graphical models for exploring gene expression data. J Multivar Anal 2004; 90: 196–212.

11.

Storey

. A direct approach to false discovery rates. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2002; 64: 479–498.

12.

Langaas

Lindqvist

Ferkingstad

. Estimating the proportion of true null hypotheses, with application to DNA microarray data. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2005; 67: 555–572.

13.

Biswas

. Estimating the proportion of true null hypotheses with application in microarray data. Communications in Statistics-Simulation and Computation 2020; 1–15. DOI: 10.1080/03610918.2020.1800036

14.

Pounds

Morris

. Estimating the occurrence of false positives and false negatives in microarray studies by approximating and partitioning the empirical distribution of

p

-values. Bioinformatics 2003; 19: 1236–1242.

15.

Wachter

. The strong limits of random matrix spectra for sample matrices of independent elements. The Annals of Probability 1978; 6: 1–18.

16.

Cai

Ren

Zhou

. Estimating structured high-dimensional covariance and precision matrices: Optimal rates and adaptive estimation. Electron J Stat 2016; 10: 1–59.

17.

Rothman

. Positive definite estimators of large covariance matrices. Biometrika 2012; 99: 733–740.

18.

Bickel

Levina

. Covariance regularization by thresholding. The Annals of Statistics 2008; 36: 2577–2604.

19.

Abu-Jamous

Kelly

. Clust: automatic extraction of optimal co-expressed gene clusters from gene expression data. Genome Biol 2018; 19: 1–11.

20.

Wang

Zhang

Huang

. Generalized gene co-expression analysis via subspace clustering using low-rank representation. BMC Bioinformatics 2019; 20: 17–27.

21.

Rao

. Advanced Statistical Methods in Biometric Research. New York: Wiley, 1952.

22.

Silva

ARD

Dias

CTDS

. A cophenetic correlation coefficient for tocher’s method. Pesqui Agropecu Bras 2013; 48: 589–596.

23.

Cheng

Gao

Tong

. Bias and variance reduction in estimating the proportion of true-null hypotheses. Biostatistics 2015; 16: 189–204.

24.

Guan

Zhao

. Nonparametric estimator of false discovery rate based on bernštĭen polynomials. Stat Sin 2008; 18: 905–923.

25.

Jiang

Doerge

. Estimating the proportion of true null hypotheses for multiple comparisons. Cancer Inform 2008; 6: 25–32.

26.

Lai

. A moment-based method for estimating the proportion of true null hypotheses and its application to microarray gene expression data. Biostatistics 2007; 8: 744–755.

27.

Markitsis

Lai

. A censored beta mixture model for the estimation of the proportion of non-differentially expressed genes. Bioinformatics 2010; 26: 640–646.

28.

Zelterman

. A parametric model to estimate the proportion from true null using a distribution for p-values. Comput Stat Data Anal 2017; 114: 105–118.

29.

Tong

Feng

Hilton

, et al. Estimating the proportion of true null hypotheses using the pattern of observed p-values. J Appl Stat 2013; 40: 1949–1964.

30.

Biswas

Chattopadhyay

Chatterjee

. Bias-corrected estimators for proportion of true null hypotheses: application of adaptive FDR-controlling in segmented failure data. J Appl Stat 2021; 1–23. DOI: 10.1080/02664763.2021.1957790

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.17 MB