Я строю кривые выживаемости CIF, используя cmprsk:
cc<-cuminc(fu_m, event, diabetes, rho=0, cencode=0,
na.action=na.omit)
plot(cc,lty=1:6,color=1, xlab="Months", ylab="Cumulative Incidence")
Таким образом, я получаю 4 кривые: событие при наличии / отсутствии диабета и соревнование при наличии / отсутствии диабета.
Вот сюжет:
Мне бы хотелось иметь только две кривые: событие при наличии / отсутствии диабета без кривых конкурирующих событий.
> str(event)
Factor w/ 3 levels "0","1","2": 1 3 2 1 2 1 1 2 1 1 ...
> str(diabetes)
Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1 1 1 ...
Событие 0: осудившиеся наблюдения,
1: событие,
2: соревнование
Диабет 0: нет диабета,
1: наличие диабета
Как я могу это исправить?
Спасибо всем