Совокупная заболеваемость, как показать только один фактор? - PullRequest
0 голосов
/ 02 декабря 2018

Я загрузил пакет "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

Но как я могу отобразить только графики "TRM" или "REL ", так что только красные или только черные, а не все вместе?

Спасибо за помощь!

1 Ответ

0 голосов
/ 03 декабря 2018

Я сделал улучшение.
При использовании термина lty = c(0, 0, 1, 2) графики исчезают.Но, к сожалению, описание до сих пор в легенде.У кого-нибудь есть идеи как его убрать?

...