[R] simple save question

David Winsemius dwinsemius at comcast.net
Wed Jul 13 18:22:32 CEST 2011


On Jul 13, 2011, at 11:16 AM, Tom La Bone wrote:

>
> No cigar. This is what I get and my session info. Any suggestions?

First think ... try this

sfit <- summary(fit, rmean = "individual")
sfit$table[5]

Hmmm. Second think. I suppose it's also possible that some of my  
unsuccessful efforts last evening are having durable effects. In  
particular, after reading the code, I tried at one point setting two  
options:

options(survfit.rmean = "individual")
options(survfit.print.rmean = "individual")

I hate writing this but since I don't know the causes for our  
different results ... Good luck.
-- 
David.

>>
>> library(survival)
>> library(ISwR)
>> dat.s <- Surv(melanom$days,melanom$status==1)
>> fit <- survfit(dat.s~1)
>> print(fit, print.rmean=TRUE)
> Call: survfit(formula = dat.s ~ 1)
>
>   records      n.max    n.start     events     *rmean *se(rmean)      
> median
>       205        205        205         57       4125         
> 161         NA
>   0.95LCL    0.95UCL
>        NA         NA
>    * restricted mean with upper limit =  5565
>> sfit <- summary(fit)
>> sfit$table[5]
> median
>    NA
>
>
>
>> sessionInfo()
> R version 2.13.1 (2011-07-08)
> Platform: i386-pc-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=English_United States.1252
> [2] LC_CTYPE=English_United States.1252
> [3] LC_MONETARY=English_United States.1252
> [4] LC_NUMERIC=C
> [5] LC_TIME=English_United States.1252
>
> attached base packages:
> [1] splines   stats     graphics  grDevices utils     datasets   
> methods
> [8] base
>
> other attached packages:
> [1] ISwR_2.0-5      survival_2.36-9 rj_0.5.5-4
>
> loaded via a namespace (and not attached):
> [1] tools_2.13.1
>
>
>
>
> --
> View this message in context: http://r.789695.n4.nabble.com/simple-save-question-tp3429148p3665289.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at 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.

David Winsemius, MD
West Hartford, CT



More information about the R-help mailing list