у меня ок. Данные за 40 лет из группы, в которой я хочу рассчитать и построить кумулятивную частоту возникновения события (DM) в пяти группах (Clustnumner) и показать их на одном графике. Я сделал начальную кривую выживания, используя код 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% индивидуумов было преобладающее событие на исходном уровне, когда началось наблюдение, как мне go с ним получить аналогичные графики? Буду признателен за помощь, так как я действительно борюсь с ней: (* 1016 *