[R] Median expected survival
Marco Barbàra
jabbba at gmail.com
Thu Jul 10 13:15:49 CEST 2014
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)
More information about the R-help
mailing list