
如图,在做生存分析,画生存曲线的时候,根据文献,画出的图应该是入图2所示(十字表示删失),但是我画出的确实图1,不存在删失,这是为什么呢?执行代码如下,用的是survival包里内置的数据集ovarian,请求大神指导。
>library(survival)
> cph.ovarian<-coxph(Surv(futime, fustat)~rx+age , ovarian)
> strata.fit<-survfit(cph.ovarian.str,newdata=data.frame(rx=c(1,2)))
> plot(strata.fit,xlab="Survival time",ylab="Survival",lty=c(1,1,2,2),col=c(1,2,1,2))
> legend(850,0.2,legend=c("rx=1,age<60","rx=2,age<60","rx=1,age>60","rx=2,age>60"),col=c(1,2,1,2),lty=c(1,1,2,2))