在coxph之后计算概率

Question

我试图弄清楚如何使用coxph包中的survival运行生存分析后计算一些概率。我已经阅读了很多帖子，但我似乎无法找到一个回答这个特定问题的帖子。

我有一个类似

的模型

coxph(Surv(time, event) ~ x1 + x2)  

而且我有兴趣计算一个在t之前生活的人在时间t+1之前死亡的可能性。

这可能吗？我想我只需要在时间t找到危险函数的估计值，但我想知道如何计算它。谢谢你的建议。

Answer 1

是的，可以通过以下代码完成：

cox.fit<- coxph(Surv(df$time,df$Status)~X1+X2+X3+X4))
summary(cox.fit)
predicted_cox.fit<-survfit(cox.fit)
Exp_Beta<- predict(predicted_cox.fit,type='lp')
basehaz_cox<-basehaz(cox.fit)

Predicted_probability_at_t1 <- as.data.frame((exp(-basehaz_cox(t1))^(exp(Exp_Beta))