我有一个大约。来自一个队列的 40 年随访数据，我想计算和绘制五组 (Clustnumner) 中事件 (DM) 的累积发生率，并希望在同一个图中显示它们。我使用代码制作了初始生存曲线
fit = survfit(Surv(Folowup_time, DM_inc) ~ as.factor(Clustnumber), data=Co_followUp，然后
plot(fit, conf.int=F, xlab = "Time in years", ylab = "Survival probability")得到以下生存曲线。每条线代表一组。
我使用代码将其转换为累积发生率图plot(fit, conf.int=F, fun = function(x) 1-x, xlab = "Time in years", ylab = "Cumulative incidence")并得到以下图
问题我的问题是，如果对于所述事件“DM”显示事件随时间的发生率，我还有另一列 (DM_B) 也显示它在基线时的流行率（在随访开始日期），我想表明图中的流行率，例如，假设我不希望我的累积发病率图从零开始，而是希望一条线从 0.3 开始，以显示 30% 的个体在随访开始时在基线发生流行事件，如何做我用它来获得类似的图表吗？帮助将不胜感激，因为我真的很努力:(