我使用cmprsk绘制CIF生存曲线:
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条曲线:存在/不存在糖尿病的事件和存在/不存在糖尿病的竞争事件。
我想只有两条曲线:在没有竞争事件的情况下存在/没有糖尿病的事件&#39;曲线。
> 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:糖尿病的存在
我该如何解决? 谢谢大家