Я загрузил пакет "cmprsk" и "Cumincidence.R" (http://www.stat.unipg.it/luca/misc/CumIncidence.R) для создания кумулятивных графиков инцидентности.
Это сработало нормально. Я использовал следующие команды:
status=factor(status, levels=c(0:2), labels=c("censored", "TRM", "REL"))
abc=factor(abc, levels=c(0,1), labels=c("A", "B"))
fit=CumIncidence(ftime, status, abc, cencode= "censored", xlab = "years")
и я получаю прикрепленный график.
![enter image description here](https://i.stack.imgur.com/CILXo.png)
Но как я могу отобразить только графики "TRM" или "REL ", так что только красные или только черные, а не все вместе?
Спасибо за помощь!