我有两个非常相似的风险评分,分别将患者分为三个风险组。将它们与对数似然比方法进行比较是否合适?
mod1 <- coxph(Surv(data$Time, data$Status==1) ~ data$Risk_1, x=T)
mod2 <- coxph(Surv(data$Time, data$Status==1) ~ data$Risk_2, x=T)
lrtest(mod1, mod2)
我的结果说,logLik(mod1)=-130,logLik(mod2)=-125(p <0.001)。因此,我得出结论:与mod1相比,mod2更适合预测生存率?
谢谢!
将它们与对数似然比进行比较是否合适方法?
我想说的太简单了。正如this论文中详述的那样,您可能想要更严格的内容,该论文已被引用30多次。在该论文中,他们比较了许多性能指标,例如简单的危险比和R [2],两组危险比和brier score,对数秩检验的p值,萨默斯等级相关性D xy,ROC曲线和曲线下面积(AUC),随机检验等。[其他相关文章可能是this,甚至可能是this引文。