r - 如何计算和绘制 30 年数据中的累积发病率（发病率）并包括 R 中的基线发生率

Question

我正在努力寻找 R 包和代码来计算和绘制 R 中的累积发病率/IR（95% CI）。我有两个数据集 D1 和 D2，每个数据集有八个事件（E1...E8）和五个组(G1...G5) 每个和；

在 D1 中，我想计算 IR（CI 为 95%），然后为每个组（G1...G5）绘制一个单独的图，显示每个组中的 E1 到 E8 IR。如下所示；

在 D2 中，我想做同样的计算，除了我的情节不应该从每个事件的零开始，因为我在研究开始日期对每个事件都有普遍的案例。它应该如下所示；

我应该从生存分析开始吗？合适吗？

Answer 1

试试下面的 Kaplan-Meier 图

library(survival)
fit = survfit(Surv(Time_variable, Censor_Variable) ~ as.factor(Stratify_Variable), data=ds)
plot(fit, conf.int=FALSE, lty=1:length(unique(Stratify_Variable)))