
23 Set
2015
23 Set
'15
15:03
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. >