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

Question

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

在 D1 中，我想计算 IR（95% CI），然后为每个组（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)))

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

How to calculate and plot cumulative incidence rate (incidence rate) over a data of 30 years and include baseline occurrence in R

r

rate

plots.jl