Abstract:High-dimensional data applications often entail the use of various statistical and machine-learning algorithms to identify an optimal signature based on biomarkers and other patient characteristics that predicts the desired clinical outcome in biomedical research. Both the composition and predictive performance of such biomarker signatures are critical in various biomedical research applications. In the presence of a large number of features, however, a conventional regression analysis approach fails to yield a good prediction model. A widely used remedy is to introduce regularization in fitting the relevant regression model. In particular, a L1 penalty on the regression coefficients is extremely useful, and very efficient numerical algorithms have been developed for fitting such models with different types of responses. This L1-based regularization tends to generate a parsimonious prediction model with promising prediction performance, i.e., feature selection is achieved along with construction of the prediction model. The variable selection, and hence the composition of the signature, as well as the prediction performance of the model depend on the choice of the penalty parameter used in the L1 regularization. The penalty parameter is often chosen by K-fold cross-validation. However, such an algorithm tends to be unstable and may yield very different choices of the penalty parameter across multiple runs on the same dataset. In addition, the predictive performance estimates from the internal cross-validation procedure in this algorithm tend to be inflated. In this paper, we propose a Monte Carlo approach to improve the robustness of regularization parameter selection, along with an additional cross-validation wrapper for objectively evaluating the predictive performance of the final model. We demonstrate the improvements via simulations and illustrate the application via a real dataset.

The effect of regularization in high dimensional Cox regression

NETWORK-REGULARIZED HIGH-DIMENSIONAL COX REGRESSION FOR ANALYSIS OF GENOMIC DATA.

Regularization paths for Cox’s proportional hazards model via coordinate descent

On the Maximum Penalized Full Likelihood Approach for Cox Model with Extreme Value for Heavily Censored Survival Data

Bayesian regularization for flexible baseline hazard functions in Cox survival models

The Impact of Regularization on High-dimensional Logistic Regression

Oracle Inequalities for the Lasso in the Cox Model

Univariate shrinkage in the Cox model for high dimensional data

PENALIZED VARIABLE SELECTION PROCEDURE FOR COX MODELS WITH SEMIPARAMETRIC RELATIVE RISK

The L1/2 regularization method for variable selection in the Cox model

Improving the Robustness of Variable Selection and Predictive Performance of Regularized Generalized Linear Models and Cox Proportional Hazard Models

Cox's proportional hazards model with a high-dimensional and sparse regression parameter

The L1/2 regularization method for variable selection in the Cox model.

Iterative L1/2 Regularization Algorithm for Variable Selection in the Cox Proportional Hazards Model

Non-asymptotic Oracle Inequalities for the High-Dimensional Cox Regression via Lasso

Adjusting For High-Dimensional Covariates In Sparse Precision Matrix Estimation By L(1)-Penalization

Asymptotic behaviour of penalized robust estimators in logistic regression when dimension increases

Diaconis-Ylvisaker prior penalized likelihood for $p/n \to κ\in (0,1)$ logistic regression

High-Dimensional Sparse Additive Hazards Regression

Adjusting for high-dimensional covariates in sparse precision matrix estimation by ℓ1-penalization