[R] Survival-Analysis: How to get numerical values from survfit (and not just a plot)?

Bernhard Reinhardt bernhard.reinhardt at dlr.de
Tue Feb 17 16:37:27 CET 2009


Hi!

I came across R just a few days ago since I was looking for a toolbox 
for cox-regression.

I´ve read
"Cox Proportional-Hazards Regression for Survival Data
Appendix to An R and S-PLUS Companion to Applied Regression" from John Fox.

As described therein plotting survival-functions works well 
(plot(survfit(model))). But I´d like to do some manipulation with the 
survival-functions before plotting them e.g. dividing one 
survival-function by another.

survfit() only returns an object of the following structure

      n events median 0.9LCL 0.9UCL
55.000 55.000  1.033  0.696  1.637

Can you tell me how I can calculate a survival- or baseline-function out 
of these values and how I extract the values from the object? I´m sure 
the calculation is done by the corresponding plot-routine, but I 
couldn´t find that one either.

Regards

Bernhard



More information about the R-help mailing list