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

[Elias Teixeira Krainski]

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