Re: [R-br] RES: Intervalos de confiança com boot.ci

O argumento 'w' e' mudado internamente na funcao boot(), ou seja, nao e' usado como sendo igual a 1 internamente. Essa alteracao e' o truque usado para o processo de reamostragem considerando estratos. Creio que vc nao precisa usar o pacote boot() para fazer reamostragem. E pode fazer reamostragem simples, pois nao vi situacao de haver estratos nos seus dados. Assim, voce poderia fazer smr.fun <- function(x) sum(x$Hosp_Death)/sum(x$SAPS3Pro2) smr0 <- smr.fun(a) smr.boot <- sapply(1:9999, function(x) smr.fun(a[sample(1:nrow(a), replace=TRUE),])) ### HPD interval (Intervalo de mais alta densidade) require(coda) hpdic <- HPDinterval(as.mcmc(smr.boot), 0.95) hist(smr.boot, breaks=pretty(smr.boot, 50)) abline(v=hpdic) Elias