Gene ranking and biomarker discovery under correlation
Open Access
- 30 July 2009
- journal article
- research article
- Published by Oxford University Press (OUP) in Bioinformatics
- Vol. 25 (20), 2700-2707
- https://doi.org/10.1093/bioinformatics/btp460
Abstract
Motivation: Biomarker discovery and gene ranking is a standard task in genomic high-throughput analysis. Typically, the ordering of markers is based on a stabilized variant of the t-score, such as the moderated t or the SAM statistic. However, these procedures ignore gene–gene correlations, which may have a profound impact on the gene orderings and on the power of the subsequent tests. Results: We propose a simple procedure that adjusts gene-wise t-statistics to take account of correlations among genes. The resulting correlation-adjusted t-scores (‘cat’ scores) are derived from a predictive perspective, i.e. as a score for variable selection to discriminate group membership in two-class linear discriminant analysis. In the absence of correlation the cat score reduces to the standard t-score. Moreover, using the cat score it is straightforward to evaluate groups of features (i.e. gene sets). For computation of the cat score from small sample data, we propose a shrinkage procedure. In a comparative study comprising six different synthetic and empirical correlation structures, we show that the cat score improves estimation of gene orderings and leads to higher power for fixed true discovery rate, and vice versa. Finally, we also illustrate the cat score by analyzing metabolomic data. Availability: The shrinkage cat score is implemented in the R package ‘st’, which is freely available under the terms of the GNU General Public License (version 3 or later) from CRAN (http://cran.r-project.org/web/packages/st/). Contact:strimmer@uni-leipzig.deOther Versions
This publication has 31 references indexed in Scilit:
- A general modular framework for gene set enrichment analysisBMC Bioinformatics, 2009
- Metabolomic profiles delineate potential role for sarcosine in prostate cancer progressionNature, 2009
- Testing significance relative to a fold-change threshold is a TREATBioinformatics, 2009
- High-dimensional classification using features annealed independence rulesThe Annals of Statistics, 2008
- A unified approach to false discovery rate estimationBMC Bioinformatics, 2008
- Simultaneous Regression Shrinkage, Variable Selection, and Supervised Clustering of Predictors with OSCARBiometrics, 2008
- Significance levels for studies with correlated test statisticsBiostatistics, 2007
- A multivariate approach for integrating genome-wide expression data and biological knowledgeBioinformatics, 2006
- Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray ExperimentsStatistical Applications in Genetics and Molecular Biology, 2004
- Significance analysis of microarrays applied to the ionizing radiation responseProceedings of the National Academy of Sciences of the United States of America, 2001