使用 bshazard 包绘制风险函数 + 生活中的风险比 table

Question

我正在尝试绘制我正在为我的博士学位所做的生存分析的风险函数，比较两种不同条件下的风险率。

我找不到使代码按预期运行的方法（here 供参考，图 4，第 7 页），以便获得平滑危险线的置信区间预测变量的水平。

我正在添加我的代码以供参考：

fitt<-bshazard(Surv(time,event) ~ session.type,data=data,lambda=10,nbin=60) 
plot(fitt,overall=FALSE, col=1, conf.int = TRUE)

函数“overall=FALSE”为我提供了两条平滑的风险曲线，均不包括置信区间，我需要从图中推断结果。这是我从代码中获得的情节图像：

如果有人知道在一段时间内获得风险率（具有上限和下限置信区间）的方法 table 以便了解每个时间间隔的这些值，那将对我有很大帮助。

感谢任何能提供帮助的人！

Answer 1

一种方法是运行函数 bshazard 按 session.type 的两个级别对数据进行分层。考虑 session.type 有两个级别（例如 0 和 1），您获得风险率（具有上下置信区间）的代码是：

-对于 session.type = 0:

Fitt0 <- bshazard(Surv(time,evento) ~1, data= data [data$session.type ==0,],lambda=10,nbin=60)
plot(fitt0,overall=TRUE, col=1, conf.int = TRUE)

-级别session.type=1

Fitt1 <- bshazard(Surv(time,evento) ~1, data= data [data$session.type ==1,],lambda=10,nbin=60)
plot(fitt1,overall=TRUE, col=1, conf.int = TRUE)

Plotting Hazard Function with bshazard package + Hazard ratios in life table