[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