Elias,

Parece que agora funcionou. Obrigado. Não tenho certeza, mas me parece que o atributo probability é a amplitude do intervalo de 95%, não é? Bom agora estou aqui pensando ... com os dados origianis que levou quase dez minutos reamostrando, talvez seja interessante um ajuste possivelmente com funções do pacote paralelo como parSapply pra fazer mais rápido. Mas isso não deve ser dificil. 

> SMR2 <- function(a){length(a$Hosp_Death[which(a$Hosp_Death==1)])/(length(a$SAPS3Pro2) * mean(a$SAPS3Pro2))}
> smr.boot <- sapply(1:9999, function(x) SMR2(a[sample(1:nrow(a), replace=TRUE),]))
> require(coda)
Carregando pacotes exigidos: coda
> hpdic <- HPDinterval(as.mcmc(smr.boot), 0.95)
> 
> hist(smr.boot, breaks=pretty(smr.boot, 50))
> abline(v=hpdic)
> hpdic
         lower    upper
var1 0.9836823 1.016965
attr(,"Probability")
[1] 0.949995
> 

Abraço forte, 

Pedro Brasil

Em 23 de setembro de 2015 11:33, Elias Teixeira Krainski <eliaskrainski@yahoo.com.br> escreveu:

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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo mnimo reproduzvel.