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

[Mauro Sznelwar]

Alguém poderia passar o código anterior, pois eu apaguei o que chegou antes.



vc pode fazer assim:

SMR2 <- function(a,w=1)
     sum(w*a$Hosp_Death)/sum(w*a$SAPS3Pro2)

SMR2(a) ### SMR observada nos dados

require(boot)

plot(boo <- boot(a,SMR2,999,stype='i'))
boot.ci(boo)

.


---
Este email foi escaneado pelo Avast antivírus.
https://www.avast.com/antivirus