Abstract:Background: An important application of microarrays is to discover genomic biomarkers, among tens of thousands of genes assayed, for disease diagnosis and prognosis. Thus it is of interest to develop efficient statistical methods that can simultaneously identify important biomarkers from such high-throughput genomic data and construct appropriate classification rules. It is also of interest to develop methods for evaluation of classification performance and ranking of identified biomarkers.Results: The ROC (receiver operating characteristic) technique has been widely used in disease classification with low dimensional biomarkers. Compared with the empirical ROC approach, the binormal ROC is computationally more affordable and robust in small sample size cases. We propose using the binormal AUC (area under the ROC curve) as the objective function for two-sample classification, and the scaled threshold gradient directed regularization method for regularized estimation and biomarker selection. Tuning parameter selection is based on V-fold cross validation. We develop Monte Carlo based methods for evaluating the stability of individual biomarkers and overall prediction performance. Extensive simulation studies show that the proposed approach can generate parsimonious models with excellent classification and prediction performance, under most simulated scenarios including model mis-specification. Application of the method to two cancer studies shows that the identified genes are reasonably stable with satisfactory prediction performance and biologically sound implications. The overall classification performance is satisfactory, with small classification errors and large AUCs.Conclusion: In comparison to existing methods, the proposed approach is computationally more affordable without losing the optimality possessed by the standard ROC method.

Nonparametric statistical inference method for partial areas under receiver operating characteristic curves, with application to genomic studies.

A Nonparametric Approach for Partial Areas under ROC Curves and Ordinal Dominance Curves

A Non-Parametric Method for the Comparison of Partial Areas under ROC Curves and Its Application to Large Health Care Data Sets.

Semi-Parametric Inference for the Partial Area under the ROC Curve

Nonparametric Covariate Adjustment for Receiver Operating Characteristic Curves

Nonparametric receiver operating characteristic curve analysis with an imperfect gold standard

Simultaneous inference for partial areas under receiver operating curves -- with a view towards efficiency

Joint Statistical Inference for the Area under the ROC Curve and Youden Index under a Density Ratio Model

Interval Estimation for the Difference in Paired Areas under the ROC Curves in the Absence of a Gold Standard Test

Statistical Inference for Box–Cox based Receiver Operating Characteristic Curves

Nonparametric ROC Summary Statistics for Correlated Diagnostic Marker Data

Semiparametric empirical likelihood confidence intervals for the difference of areas under two correlated ROC curves under density ratio model

Small-sample precision of ROC-related estimates

Bayesian nonparametric inference for the covariate-adjusted ROC curve

Nonparametric Estimation of the ROC Curve Based on the Bernstein Polynomial

A comparison of confidence/credible interval methods for the area under the ROC curve for continuous diagnostic tests with small sample size

Regularized Binormal ROC Method in Disease Classificationusing Microarray Data.

A Marginal Model Approach for Analysis of Multi-Reader Multi-Test Receiver Operating Characteristic (ROC) Data

A non-inferiority test of areas under two parametric ROC curves

A Non-Parametric Method for Building Predictive Genetic Tests on High-Dimensional Data

ROC Curve Estimation under Test-Result-dependent Sampling.