目标
我正在尝试在R中绘制生存数据的危险函数。
尝试
我正在使用R survival软件包中的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生存数据的患者队列的危害函数？