目标
我正在尝试在 R 中绘制生存数据的危险函数。
试图
survival我正在使用R 包中的 NCCTG 肺癌数据集。
据我了解，危险函数具有以下形式：
h(t) = h0(t) * e^(b1 * x1 + ... + bn * xn)
x1，...，xn 的值是从数据自身导出的自变量，使用以下代码，我可以获得系数 b1，...，bn 的值。
# Sample code
library(survival)
lung$status <- lung$status - 1
coxph_mdl <- coxph( Surv(time, status ) ~ age, data = lung )
summary(coxph_mdl)
问题
如何在 R 中绘制一组具有生存数据的患者的风险函数？