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

[Pedro Emmanuel Alvarenga Americano do Brasil]

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 em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150923/78d8ac98/attachment.html>