Gene ranking and biomarker discovery under correlation

Abstract

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. 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. The shrinkage cat score is implemented in the R package "st" available from URL http://cran.r-project.org/web/packages/st/

18 pages, 5 figures, 1 table

References in corpus (5)

• In silico prediction of protein-protein interactions in human macrophages
• Classifier Technology and the Illusion of Progress
• High-dimensional classification using features annealed independence rules
• Feature selection in omics prediction problems using cat scores and false nondiscovery rate control
• Correlation-sharing for detection of differential gene expression

Cited by in corpus (10)

• Optimal whitening and decorrelation
• Feature selection in omics prediction problems using cat scores and false nondiscovery rate control
• A method for generating realistic correlation matrices
• Signal Identification for Rare and Weak Features: Higher Criticism or False Discovery Rates?
• Optimal feature selection for sparse linear discriminant analysis and its applications in gene expression data
• A novel algorithm for simultaneous SNP selection in high-dimensional genome-wide association studies
• The Limitations of Simple Gene Set Enrichment Analysis Assuming Gene Independence
• Analysis of Gene Interaction Graphs as Prior Knowledge for Machine Learning Models
• Classification of high-dimensional data with spiked covariance matrix structure
• Variable Selection Using Relative Importance Rankings