RV: help

Jani Lobo <lobo <at> grupocomar.com> writes:

I want to estimate the survival mean of a few specific teams. I'm trying to
calculate it through a Kaplan Meier estimator. For doing so, I load the
"survival" package and run the following instructions: 

"options(survfit.print.mean=TRUE)"      allows showing the mean and mean
standard error 

"KM=survfit(Surv(Dias,Censura))"          runs the model

"KM"                                     shows the model

It is in this very moment when it shows the mean value and the mean standard
error.

I'd like to know if there is any instruction that returns the mean value and
its standard error, something like:

"KM$rmean"

It's always good to use an example that can be run by other people, 
and it is easy: just use the example from the documentation

library(survival)
fit <- survfit(Surv(time, status) ~ x, data=aml)
# str tells you what fit has to offer
str(fit)

### Too bad, no summary data available here, looks
## like print.survfit does the job
getAnywhere(print.survfit)

Checking the code: there are a lot of print statements in here, 
and the summary results are not returned. Now you have two choices: 
either use the code from print.survfit, copy it, and get the values
you need from a changed function. Can be some work.

Or check if someone else has done the work for you, and bingo:

summary.survfit in package Design looks like it is doing the job.

Dieter
















summary.survfit {Design}