r - 在 R 中使用 Cox Proportional Hazard 模型计算生存预测

Question

我正在尝试使用 R 中的 Cox Proportional Hazard 模型计算生存预测。

    library(survival)
    data(lung)
    model<-coxph(Surv(time,status ==2)~age + sex + ph.karno + wt.loss, data=lung)
    predict(model, data=lung, type ="expected")

当我使用上面的代码时，我得到了与公式对应的累积风险的预测

    h^i(t)=h^0(t)exp(x′iβ^)

但我关心的是预测对应于公式的生存，

    S^i(t)=S^0(t)exp(x′iβ^)

我如何预测 R 中的生存？
提前致谢。

Answer 1

您可以使用 predict 或 survfit 。使用 predict 您需要为 newdata 参数提供一个包含模型中所有变量值的列表:

predict(model, 
      newdata=list(time=100,status=1,age=60,sex=1, ph.karno=60,wt.loss=15),
      type ="expected")
[1] 0.2007497

survfit 对象有一个绘图方法:

?survreg
png(); plot(survfit(model)); dev.off()