An Efficient LASSO Estimation for Time‑Dependent Cox Models via Adaptive Active‑Set Coordinate Descent

Document Type : Research article

Authors

1 Department of Statistics, Faculty of Mathematical Sciences, Shahrood University of Technology, Shahrood, Iran

2 Department of Statistics, Faculty of Mathematical Sciences, Ferdowsi University of Mashhad, P. O. Box 1159, Mashhad 91775, Iran.

Abstract

The penalized Cox proportional hazard model is a popular analytical approach for survival data with a large set of covariates. Such problems are especially challenging when covariates vary over follow-up time (i.e., the covariates are time-dependent), leading to increased computational complexity and difficulties in efficient variable selection. In this paper, we propose an Adaptive Active-Set Coordinate Descent (AACD) algorithm for LASSO-penalized time-dependent Cox models. The proposed method combines coordinate descent with an adaptive active-set strategy and warm starts along the regularization path, allowing the algorithm to focus computation on relevant variables and better exploit sparsity. Simulation studies demonstrate that AACD consistently outperforms \texttt{glmnet}, achieving uniformly lower mean squared errors with relative efficiency greater than one across all settings, and reducing estimation error by up to 37\% in small samples. The method remains robust as dimensionality increases, improves variable selection by reducing false negatives, and produces more parsimonious models with comparable predictive performance while reducing runtime by approximately 36\% in real data analysis.

Keywords

Main Subjects


 

Article PDF

[1] P. K. Andersen and R. D. Gill, Cox’s regression model for counting processes: A large sample study, The Annals of Statistics, 10(4), 1100–1120, (1982).
[2] J. Bradic, J. Fan and J. Jiang, Regularization for Cox’s proportional hazards model with NP-dimensionality, The Annals of Statistics, 39(6), 3092–3120, (2011).
[3] N. Breslow, Covariance analysis of censored survival data, Biometrics, 30(1), 89–99, (1974).
[4] D. R. Cox, Regression models and life-tables, Journal of the Royal Statistical Society: Series B (Methodological), 34(2), 187–202, (1972).
[5] J. Fan and R. Li, Variable selection via nonconcave penalized likelihood and its oracle properties, Journal of the American Statistical Association, 96(456), 1348–1360, (2001).
[6] J. H. Friedman, T. Hastie and R. Tibshirani, Regularization paths for generalized linear models via coordinate descent, Journal of Statistical Software, 33(1), 1–22, (2010).
[7] J. J. Goeman, L1 penalized estimation in the Cox proportional hazards model, Biometrical Journal, 52(1), 70–84, (2010).
[8] J. Gui and H. Li, Penalized Cox regression analysis in the high-dimensional and low-sample size settings, with applications to microarray gene expression data, Bioinformatics, 21(13), 3001–3008, (2005).
[9] F. E. Harrell, R. M. Califf, D. B. Pryor, K. L. Lee and R. A. Rosati, Evaluating the yield of medical tests, JAMA, 247(18), 2543–2546, (1982).
[10] P. J. Heagerty, T. Lumley and M. S. Pepe, Time-dependent ROC curves for censored survival data and a diagnostic marker, Biometrics, 56(2), 337–344, (2000).
[11] J. Huang, Y. Jiao, B. Jin, J. Liu, X. Lu and C. Yang, A unified primal–dual active set algorithm for nonconvex sparse recovery, Statistical Science, 36(2), 215–238, (2021).
[12] N. Simon, J. Friedman, T. Hastie and R. Tibshirani, Regularization paths for Cox’s proportional hazards model via coordinate descent, Journal of Statistical Software, 39(5), 1–13, (2011).
[13] I. Sohn, J. Kim, S.-H. Jung and C. Park, Gradient lasso for Cox proportional hazards model, Bioinformatics, 25(14), 1775–1781, (2009).
[14] T. M. Therneau and T. Lumley, Package ’survival’, R Package Version, 2.38, 1–37, (2015).
[15] T. M. Therneau, C. S. Crowson and E. Atkinson, Using time-dependent covariates and time-dependent coefficients in the Cox model, Survival Vignettes, 1, 1–25, (2017).
[16] R. Tibshirani, The LASSO method for variable selection in the Cox model, Statistics in Medicine, 16(4), 385–395, (1997).
[17] R. Tibshirani, J. Bien, J. Friedman, T. Hastie, N. Simon, J. Taylor and R. J. Tibshirani, Strong rules for discarding predictors in lasso-type problems, Journal of the Royal Statistical Society: Series B (Statistical Methodology), 74(2), 245–266, (2012).
[18] P. Tseng, Convergence of a block coordinate descent method for nondifferentiable minimization, Journal of Optimization Theory and Applications, 109(3), 475–494, (2001).
[19] Y. Yang and H. Zou, A cocktail algorithm for solving the elastic net penalized Cox’s regression in high dimensions, Statistics and Its Interface, 6(2), 167–173, (2013).
[20] H. Zou and T. Hastie, Regularization and variable selection via the elastic net, Journal of the Royal Statistical Society: Series B (Statistical Methodology), 67(2), 301–320, (2005). 
Volume 11, Issue 2
September 2026
Pages 273-288
  • Receive Date: 17 April 2026
  • Revise Date: 24 May 2026
  • Accept Date: 18 June 2026
  • Publish Date: 04 September 2026