I have developed a cox-model in my research, and I want to compare it to the conventional staging system for prognostic evaluation. The traditional staging system categorizes patients into risk groups with different prognoses. To compare my new model with the conventional model, I have categorized my model prediction into risk groups with the same size (n) as the traditional model. The Kaplan Meier analysis of the traditional and new models shows that both are significant predictors of prognostic outcome ( Log-rank test p<.0.001), but the new model is better for placing patients with poor prognoses at the high-risk groups. My question is: What is the best statistical approach to compare these models?
My solutions so far have been to 1) fit the model group/predictor in a cox-model and compare them with an ANOVA-test (this yields a p-value < 0.001.
- perform a comparison of the concordance indexes for the two models predictors: compareC(time,event, new.model.predictor,conventional.model.predictor), concordance index:0.86 vs 0.75 with p-value = 0.02.
I guess approach nr.1 compares the model fit, while approach 2 compares the discriminatory abilities of the two models.
Inputs on which approach would be the most informative and appropriate would be greatly appreciated.