Hi, Lars. I don't understand well your question. why don't you simply type
pred_leuk Call: survfit(formula = leuk.cox, newdata = leuk_new) records n.max n.start events median 0.95LCL 0.95UCL 1 42 42 42 30 1 1 23 2 42 42 42 30 7 5 NA 3 42 42 42 30 15 10 NA 4 42 42 42 30 NA NA NA 5 42 42 42 30 8 5 NA 6 42 42 42 30 23 22 NA 7 42 42 42 30 8 6 NA 8 42 42 42 30 NA 22 NA 9 42 42 42 30 6 4 NA 10 42 42 42 30 23 17 NA Yes, the quantile method returns the correct medians. The NAs usually appear when the unit belongs to a category for which median survival time has not yet been reached, so there is no way for the model to estimate the median time. You can see the expected survival curve with plot(predl_leuk) ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.