abdul kudus
2005-Dec-19 11:48 UTC
[R] Package "boot": How to construct CI from censboot object?
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.caseSTRATIFIED 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]]