Statistical inference after variable selection in Cox models: a neutral simulation study
- Choosing relevant predictors is central to the analysis of biomedical time-to-event data. Classical frequentist inference, however, presumes that the set of covariates is fixed in advance and does not account for data-driven variable selection. As a consequence, naive post-selection inference may be biased and misleading. In right-censored survival settings, these issues may be further exacerbated by the additional uncertainty induced by censoring. We investigate several inference procedures applied after variable selection for the coefficients of the Lasso and its extension, the adaptive Lasso, in the context of the Cox model. The methods considered include sample splitting, post-selection inference procedures that condition explicitly on the Lasso selection event, and the debiased Lasso. Because these methods address different inferential targets, we distinguish selected-submodel targets from full-model targets and interpret empirical coverage, interval width, power, and type I errorChoosing relevant predictors is central to the analysis of biomedical time-to-event data. Classical frequentist inference, however, presumes that the set of covariates is fixed in advance and does not account for data-driven variable selection. As a consequence, naive post-selection inference may be biased and misleading. In right-censored survival settings, these issues may be further exacerbated by the additional uncertainty induced by censoring. We investigate several inference procedures applied after variable selection for the coefficients of the Lasso and its extension, the adaptive Lasso, in the context of the Cox model. The methods considered include sample splitting, post-selection inference procedures that condition explicitly on the Lasso selection event, and the debiased Lasso. Because these methods address different inferential targets, we distinguish selected-submodel targets from full-model targets and interpret empirical coverage, interval width, power, and type I error accordingly. Their performance is examined in a neutral simulation study reflecting realistic covariate structures and censoring rates commonly encountered in biomedical applications. The primary focus is post-selection inference after Cox-Lasso variable selection, not a comprehensive benchmark of very-high-dimensional variable-selection performance. To complement the simulation results, we illustrate the practical behavior of these procedures in an applied example using a publicly available survival dataset.…
