Как рассчитать и построить кумулятивный коэффициент заболеваемости (уровень заболеваемости) за данные за 30 лет и включить исходный уровень в R - PullRequest
Новая
       21

Как рассчитать и построить кумулятивный коэффициент заболеваемости (уровень заболеваемости) за данные за 30 лет и включить исходный уровень в R

0 голосов
/ 09 января 2020

Я изо всех сил пытаюсь найти пакет R и коды для расчета и построения кумулятивного коэффициента заболеваемости / IR (с 95% CI) в R. У меня есть два набора данных D1 и D2 с восемью событиями в каждом (E1 ... E8) и пять групп (G1 ... G5) каждая;

В D1 я хочу рассчитать IR (с 95% -ным доверительным интервалом), а затем иметь отдельный график для каждой группы (G1 ... G5), показывая E1 до E8 IR в каждом. Примерно так:

enter image description here

И в D2 я хочу делать те же вычисления, за исключением того, что мои графики не должны начинаться с нуля для каждого события, так как я Распространенные случаи для каждого события на дату начала исследования. Это должно быть примерно так:

enter image description here

Должен ли я начать с анализа выживания? это подходит?

1 Ответ

0 голосов
/ 09 января 2020

Попробуйте следующий график Каплана-Мейера 

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