如何从生存数据中绘制危害函数？

Question

目标

我正在尝试在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生存数据的患者队列的危害函数？