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@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

______________________________________________
R-help@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.

Reply via email to