i'm new r , trying plot survfit
survival curves.
in playing around survfit
object, found 2 different plots following:
library(survival) #example survfit object mysurvfit <- survfit(surv(time, status)~1, data=aml) #default survfit plot, survival curve upper & lower conf intervals plot(mysurvfit, mark.time=false, conf.int=true) #create curve accessing surv, upper, lower #(i'd expect produce same above, doesn't) lines(mysurvfit$surv, col="blue",lty=1) lines(mysurvfit$upper, col="blue",lty=2) lines(mysurvfit$lower, col="blue",lty=2)
why these curves different? missing survfit object?
you missing time
variable
try
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)
which looks like
Comments
Post a Comment