Re: [R] calculating quantiles

2013-10-30 Thread Terry Therneau 



On 10/30/2013 06:00 AM, r-help-requ...@r-project.org wrote:

Hi,

I'm having the following loop:

result<- vector("list",100)
for (i in 1:max(dat$simNumber))
{  
result[[i]]<-survfit(Surv(dat[dat$simNumber==i,]$TAFD,dat[dat$simNumber==i,]$DV)~1)

}

In a next step, I would like to calculate the mean, 5% and 95% PI of the 
Kaplan-Meier estimates of the 100 simulated problems stored under result in the 
result[[i]]$surv.

Can anyone help me with this?

Tx!



I don't know what a "95% PI" is, but for the mean survival I can help.

options(survfit.rmean="common")  #defaults to "none" in the routine
means <- lapply(result, function(x) summary(x)$table[5:6])

This gives a list, each element of which is the estimated mean and se(mean) for 
that curve.

Terry Therneau

__
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.

[R] calculating quantiles

2013-10-29 Thread Rossenu, Stefaan 
Hi,

I'm having the following loop:

result <- vector("list",100)
for (i in 1:max(dat$simNumber))
{  
result[[i]]<-survfit(Surv(dat[dat$simNumber==i,]$TAFD,dat[dat$simNumber==i,]$DV)~1)

}

In a next step, I would like to calculate the mean, 5% and 95% PI of the 
Kaplan-Meier estimates of the 100 simulated problems stored under result in the 
result[[i]]$surv.

Can anyone help me with this?

Tx!



Notice:  This e-mail message, together with any attachme...{{dropped:14}}

__
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.