PDF(608 KB)
Group MCP for Cox Models with Time-Varying Coefficients
Xiaodong XIE, Shaozhi ZHENG
Journal of Systems Science and Information ›› 2016, Vol. 4 ›› Issue (5) : 476-488.
PDF(608 KB)
PDF(608 KB)
Group MCP for Cox Models with Time-Varying Coefficients
Cox's proportional hazard models with time-varying coefficients have much flexibility for modeling the dynamic of covariate effects. Although many variable selection procedures have been developed for Coxs proportional hazard model, the study of such models with time-varying coefficients appears to be limited. The variable selection methods involving nonconvex penalty function, such as the minimax concave penalty (MCP), introduces numerical challenge, but they still have attractive theoretical properties and were indicated that they are worth to be alternatives of other competitive methods. We propose a group MCP method that uses B-spline basis to expand coefficients and maximizes the log partial likelihood with nonconvex penalties on regression coefficients in groups. A fast, iterative group shooting algorithm is carried out for model selection and estimation. Under some appropriate conditions, the simulated example shows that our method performs competitively with the group lasso method. By comparison, the group MCP method and group lasso select the same amount of important covariates, but group MCP method tends to outperform the group lasso method in selection of unimportant covariates.
group MCP / group lasso / varying coefficient / B-spline / Cox models {{custom_keyword}} /
[1] Cox D R. Regression models and life-tables (with Discussion). Journal of the Royal Statistical Society Series B, 1972, 34(2):187-220.
[2] Cox D R. Partial likelihood. Biometrika, 1975, 62(2):269-76.
[3] Sauerbrei W, Schumacher M. A bootstrap resampling procedure for model building:Application to the Cox regression model. Statistics in Medicine, 1992, 11(16):2093-2109.
[4] Faraggi D, Simon R. Bayesian variable selection method for censored survival data. Biometrics, 1998, 54(4):1475-1485.
[5] Ibrahim J G, Chen M H, Maceachern S N. Bayesian variable selection for proportional hazards models. The Canadian Journal of Statistics, 1999, 27(4):701-717.
[6] Fan J Q, Li R Z. Variable selection for Coxs proportional hazards model and frailty model. The Annals of Statistics, 2002, 30(1):74-99.
[7] Tibshirani R. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society Series B, 1996, 58(1):267-288.
[8] Tibshirani R. The lasso method for variable selection in the Cox model. Statistics in Medicine, 1997, 16(4):385-395.
[9] Fan J Q, Li R Z. Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association, 2001, 96(456):1348-1360.
[10] Zhang H H, LuWB. Adaptive lasso for Coxs proportional hazards model. Biometrika, 2007, 94(3):691-703.
[11] Hastie T, Tibshirani R. Varying-coefficient models. Journal of the Royal Statistical Society Series B, 1993, 55(4):757-796.
[12] Lin Y, Zhang H H. Component selection and smoothing in multivariate nonparametric regression. Annals of Statistics, 2006, 34(5):2272-2297.
[13] Leng C L, Zhang H H. Model selection in nonparametric hazard regression. Journal of Nonparametric Statistics, 2006, 18(7-8):417-429.
[14] Huang J, Horowitz J L, Wei F R. Variable selection in nonparametric additive models. Annals of Statistics, 2010, 38(4):2282-2313.
[15] Yan J, Huang J. Model selection for Cox models with time-varying coefficients. Biometrics, 2012, 68(2):419-428.
[16] Li H Z, Luan Y H. Kernel Cox regression models for linking gene expression profiles to censored survival data. Pacific Symposium of Biocomputing, 2003, 8:65-76.
[17] Li H Z, Gui J. Partial Cox Regression Analysis for High-Dimensional Microarray Gene Expression Data. Bioinformatics, 2004, 20(1):208-215.
[18] Tibshirani R. Univariate shrinkage in the Cox model for high dimensional data. Statistical Applications in Genetics and Molecular Biology, 2009, 8(1):3498-3528.
[19] Bradic J, Fan J Q, Jiang J C. Regularization for Coxs proportional hazards model with NP-dimensionality. Annals of Statistics, 2011, 39(6):3092-3120.
[20] Fan J Q, Feng Y, Wu Y C. High-dimensional variable selection for Coxs proportional hazards model. Institute of Mathematical Statistics, 2010, 6:70-86.
[21] Yang G R, Huang J, Zhou Y. Concave group methods for variable selection and estimation in highdimensional varying coeffcient models. Science China Mathematics, 2014, 57(10):2073-2090.
[22] Zucker D M, Karr A F. Nonparametric survival analysis with time-dependent covariate effects:A penalized partial likelihood approach. The Annals of Statistics, 1990, 18(1):329-353.
[23] Zhang C H. Nearly unbiased variable selection under minimax concave penalty. The Annals of Statistics, 2010, 38(2):894-942.
[24] Breheny P, Huang J. Coordinate descent algorithms for nonconvex penalized regression, with applications to biological feature selection. Annals of Applied Statistics, 2011, 5(1):232-253.
[25] Yuan M, Lin Y. Model selection and estimation in regression with grouped variables. Journal of the Royal Statistical Society Series B, 2006, 68(1):49-67.
[26] R Development Core Team. R:A Language and Environment for Statistical Computing. Vienna:R Foundation for Statistical Computing, 2011.
[27] Breslow N. Covariance analysis of censored survival data. Biometrics, 1974, 30(1):89-99.
[28] Tseng P, Yun S. Block-coordinate gradient descent method for linearly constrained nonsmooth separable optimization. Journal of Optimization Theory and Applications, 2009, 140(3):513-535.
[29] Craven P, Wahba G. Smoothing noisy data with spline functions. Numerische Mathematik, 1978, 31(4):377-403.
[30] Perperoglou A, le Cessie S, van Houwelingen H C. A fast routine for fitting Cox models with time varying effects of the covariates. Biomedicine, 2006, 81(2):154-161.
PDF(608 KB)
/
| 〈 |
|
〉 |
