Dear all,

I run the example of "censboot" contained in "boot" package.  But, I can't
find the confidence interval of the resulted "censboot" object.  Any idea ?

> aml.fun <- function(data) {
+           surv <- survfit(Surv(time, cens)~group, data=data)
+           out <- NULL
+           st <- 1
+           for (s in 1:length(surv$strata)) {
+                inds <- st:(st+surv$strata[s]-1)
+                md <- min(surv$time[inds[1-surv$surv[inds]>=0.5]])
+                st <- st+surv$strata[s]
+                out <- c(out,md)
+           }
+           out
+      }

> aml.case <- censboot(aml,aml.fun,R=499,strata=aml$group)

> aml.case

STRATIFIED CASE RESAMPLING BOOTSTRAP FOR CENSORED DATA


Call:
censboot(data = aml, statistic = aml.fun, R = 499, strata = aml$group)


Bootstrap Statistics :
    original    bias    std. error
t1*       31       Inf         NaN
t2*       23 -2.501002     8.67887


> aml.ci <- boot.ci(aml.case,conf=c(0.90,0.95),type="all")
Error in if (stype == "w") type <- "inf" else if (!is.null(boot.out) &&  :
        argument is of length zero


regards,

abdul kudus
INSPEM - UPM
Malaysia

	[[alternative HTML version deleted]]


