用 bshazard 包 + 生命表中的危险比绘制危险函数答案

【问题标题】：Plotting Hazard Function with bshazard package + Hazard ratios in life table用 bshazard 包 + 生命表中的危险比绘制危险函数
【发布时间】：2021-12-01 19:55:23
【问题描述】：

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

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

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

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

函数“overall=FALSE”给了我两条平滑的危险曲线，两者都不包括置信区间，我需要从图中推断结果。这是我从代码中获得的绘图图像：

如果有人知道一种方法来获取时间表中的危险率（具有上下置信区间）以便了解每个时间间隔的这些值，那将对我有很大帮助。

感谢任何可以提供帮助的人！

【问题讨论】：

标签： r survival-analysis confidence-interval survival hazard

【解决方案1】：

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

-for 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)

-for 级别 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)

【讨论】：