在 coxph 之后计算概率答案

【问题标题】：Compute probabilities after coxph在 coxph 之后计算概率
【发布时间】：2015-01-29 16:13:55
【问题描述】：

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

我有一个类似的模型

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

而且我有兴趣计算一个活到时间 t 的人会在时间 t+1 之前死去的概率。

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

【问题讨论】：

这听起来更像是统计帮助而不是编码帮助。你能澄清一下你需要什么吗？
公平点@rawr。我会仔细检查我的统计数据，但假设我确实想计算风险函数 h(t)（不是集成/累积风险函数）的估计值，有没有办法做到这一点？我可以使用 predict.coxph 和 arg type = 'risk' 来计算 exp(X*beta)，然后将其乘以 basehaz 的基线危害吗？谢谢！

标签： r survival-analysis

【解决方案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))

【讨论】：