<div><span style="color: #0000ff; font-size: medium;">Muito obrigado pela resposta! Vou tentar instalá-lo.</span></div>
<div><br />O INLA nóa est[a dispon[ivel no CRAN<br /><a href="http://www.r-inla.org" target="_blank">www.r-inla.org</a><br /><br /><br /><br />On Thu, 11 Jun 2015, Mauro Sznelwar wrote:<br /><br />> Tentei instalar esta biblioteca INLA e não consegui, por acaso é exclusiva<br />> do LINUX ou não está disponível no CRAN?<br />><br />><br />> Como os dados são contínuos seria necessário considerar uma distribuição<br />> contínua com mistura Bernoulli, veja uma proposta aqui:<br />> <a href="http://jupiter.est.ufmg.br/~posgrad/mestrado/dissertacao_mestrado_zaida%20_q" target="_blank">http://jupiter.est.ufmg.br/~posgrad/mestrado/dissertacao_mestrado_zaida%20_q</a><br />> uiroz.pdf<br />> e um exemplo no script abaixo:<br />><br />> ## hurdle model for zero-gamma (zero inflation and continuous outcome)<br />> require(INLA)<br />><br />> ### rainfall data from Parana' state in Brazil<br />> data(PRprec)<br />><br />> ### analize data from one location<br />> PRprec[301,1:3]<br />> y ><br />> summary(y)<br />><br />> plot.ts(y)<br />> if (TRUE)<br />> dev.off()<br />><br />> ### organize the dataset to build a model for<br />> ### smoothing the rainfall amount over the year<br />> ### accounting for a smoothed zero probability as well<br />> n > dat > c((y>0)+0, rep(NA,n)),<br />> c(rep(NA,n), ifelse(y>0,y,NA))),<br />> t0=c(1:n, rep(NA,n)),<br />> t1=c(rep(NA,n), 1:n))<br />><br />> ### two choices to smooth the rainfall amount<br />> form1 > f(t0, model='rw2', cyclic=TRUE) +<br />> f(inla.group(t1, n=73), model='rw1', cyclic=TRUE)<br />> form2 > f(t0, model='rw2', cyclic=TRUE) +<br />> f(t1, model='rw2', cyclic=TRUE)<br />> inla.setOption(scale.model=TRUE)<br />><br />> res1 > res2 ><br />> day > par(mfrow=c(2,1), mar=c(3,3,1,1), mgp=c(2,1,0))<br />> plot(day, y>0)<br />> for (j in 4:6) {<br />> lines(day, 1/(1+exp(-res1$summary.ran$t0[,j])), lty=2)<br />> lines(day, 1/(1+exp(-res2$summary.ran$t0[,j])), lty=2, col=2)<br />> }<br />> plot(day, y, type='l')<br />> for (j in 4:6) {<br />> lines(day[seq(3,n,5)], exp(res1$summary.ran[[2]][,j]), lty=2)<br />> lines(day, exp(res2$summary.ran[[2]][,j]), lty=2, col=2)<br />> }<br />><br />><br />> Elias<br />> _______________________________________________<br />> R-br mailing list<br />> R-br@listas.c3sl.ufpr.br<br />> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br />> mínimo reproduzível.<br />><br />><br />> ---<br />> Este email foi escaneado pelo Avast antivírus.<br />> <a href="https://www.avast.com/antivirus" target="_blank">https://www.avast.com/antivirus</a><br />><br />> _______________________________________________<br />> R-br mailing list<br />> R-br@listas.c3sl.ufpr.br<br />> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br />>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div>