# Model Selection with Competing risks in Cox regression

When doing cox proportional hazards regression one often has competing risks. The typical approach for this is to fit separate cox proportional hazards models for each risk, censoring the competing risk.

How do you do model selection in this scenario? Should the two models be the same? Do I build a model based on one of the risks and use the same model for the other risk? Do I build a model based off a completely different response (e.g. overall survival)?

• You need to define what "build a model" means. If it means estimating regression coefficients then your approach to fit a series of models with different target outcomes seems sound. If you mean "variable selection" then you have to question why you would spend time doing variable selection and not using subject matter expertise to specify the model. Nov 27, 2014 at 13:36
• @FrankHarrell I was referring to variable selection. This is for a project in school on a real dataset and we don't have access to the "subject matter experts". However, I think "you should only consult the experts" is kind of a cop out answer because that belittles the entire field of model selection. Nov 27, 2014 at 14:15
• The entire field of model selection is full of false promises. Please read elsewhere on this site about the many pitfalls and high probability of getting a misleading result. Or just look at stata.com/support/faqs/statistics/stepwise-regression-problems Nov 27, 2014 at 17:23
• I didn't suggest stepwise procedures either, but rather AIC/BIC type of criteria, but you have a valid argument and I can appreciate your experience. Happy Thanksgiving :) Nov 27, 2014 at 23:32
• BIC will underfit. AIC works well if you have say 3 competing models to choose from. AIC is just a (smarter) restatement of $P$-values so doesn't solve the fundamental problem if the number of competing models exceeds a few. Nov 28, 2014 at 12:24

However, you wondered whether you should consider overall survival in modeling competing events. This is another approach, based on Fine and Gray (1999) cause-specific ${subdistribution}$ (i.e. the cumulative incidence function). Here, the drawback is that they consider in the risk set for a given event also individuals that may not have the event itself (because the competing one has already occurred previously), thus, their estimates are not based on probability in a proper sense. However, their approach has a direct interpretation in terms of survival probabilities: it allows to calculate the probability of each event. For example, a Fine and Gray approach may lead us to conclude that, in a given sample, with a given covariate pattern, $25\%$ of patients will die for cardiological reasons, $20\%$ of patients for other reasons, and $55\%$ will survive until the end of the follow-up.