<div dir="ltr"><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)">Elias,</div><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)"><br></div><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)">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. </div><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)"><br></div><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)"><div class="gmail_default">> SMR2 <- function(a){length(a$Hosp_Death[which(a$Hosp_Death==1)])/(length(a$SAPS3Pro2) * mean(a$SAPS3Pro2))}</div><div class="gmail_default">> smr.boot <- sapply(1:9999, function(x) SMR2(a[sample(1:nrow(a), replace=TRUE),]))</div><div class="gmail_default">> require(coda)</div><div class="gmail_default">Carregando pacotes exigidos: coda</div><div class="gmail_default">> hpdic <- HPDinterval(as.mcmc(smr.boot), 0.95)</div><div class="gmail_default">> </div><div class="gmail_default">> hist(smr.boot, breaks=pretty(smr.boot, 50))</div><div class="gmail_default">> abline(v=hpdic)</div><div class="gmail_default">> hpdic</div><div class="gmail_default">         lower    upper</div><div class="gmail_default">var1 0.9836823 1.016965</div><div class="gmail_default">attr(,"Probability")</div><div class="gmail_default">[1] 0.949995</div><div class="gmail_default">> </div><div class="gmail_default"><br></div><div class="gmail_default">Abraço forte, </div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="'courier new', monospace" style="background-color:rgb(255,255,255)" color="#000066">Pedro Brasil</font></div></div></div></div></div></div>
<br><div class="gmail_quote">Em 23 de setembro de 2015 11:33, Elias Teixeira Krainski <span dir="ltr"><<a href="mailto:eliaskrainski@yahoo.com.br" target="_blank">eliaskrainski@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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.<br>
<br>
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<br>
<br>
smr.fun <- function(x) sum(x$Hosp_Death)/sum(x$SAPS3Pro2)<br>
smr0 <- smr.fun(a)<br>
smr.boot <- sapply(1:9999, function(x) smr.fun(a[sample(1:nrow(a), replace=TRUE),]))<br>
<br>
### HPD interval (Intervalo de mais alta densidade)<br>
require(coda)<br>
hpdic <- HPDinterval(as.mcmc(smr.boot), 0.95)<br>
<br>
hist(smr.boot, breaks=pretty(smr.boot, 50))<br>
abline(v=hpdic)<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
Elias</font></span><div class="HOEnZb"><div class="h5"><br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</div></div></blockquote></div><br></div>