Блог

Я новичок в R и пытаюсь построить survfit кривые выживания.

Поиграв с survfit объектом, я обнаружил, что получаю 2 разных графика для следующего:

 library(survival)

#example survfit object
mysurvfit <- survfit(Surv(time, status)~1, data=aml)

#default survfit plot, survival curve with upper amp; lower conf intervals
plot(mysurvfit, mark.time=FALSE, conf.int=TRUE)

#create another curve by accessing surv, upper, lower
#(I'd expect this to produce the same as above, but it doesn't)
lines(mysurvfit$surv, col="blue",lty=1)
lines(mysurvfit$upper, col="blue",lty=2)
lines(mysurvfit$lower, col="blue",lty=2)
  

Почему эти кривые отличаются? Чего я не понимаю в объекте survfit?

Вам не хватает time переменной

Попробуйте

 plot(mysurvfit, mark.time=FALSE, conf.int=TRUE)
lines(mysurvfit$surv ~ mysurvfit$time, col="blue",lty=1)
lines(mysurvfit$upper ~ mysurvfit$time, col="blue",lty=2)
lines(mysurvfit$lower ~ mysurvfit$time, col="blue",lty=2)
  

который выглядит как