>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))


雷达卡



京公网安备 11010802022788号







