<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:14pt"><div><span><br></span></div><div>Muito obrigada, Elias.</div><div>Abraços<br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 14pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Elias T. Krainski <eliaskrainski@yahoo.com.br><br> <b><span style="font-weight: bold;">Para:</span></b> "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quarta-feira, 21 de Novembro de 2012 14:07<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] INLA Geoestatística<br> </font> </div> <br>Oi Taynãna, <br><br>DIC negativo pode ocorrer com distribuições contínuas sem problemas. <br>O mesmo
vale para log-Verossimilhança. <br>Veja sobre DIC em <br><a href="http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/dicpage.shtml" target="_blank">http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/dicpage.shtml</a><br><br>### exemplo de log-verossimilhança negativa <br> y1 <- rnorm(100, 0, 100)<br>### log-verossimilhança em 0, 100<br> sum(dnorm(y1, 0, 100, log=TRUE))<br><br>### exemplo de log-verossimilhança positiva<br> y2 <- rnorm(100, 0, 0.01)<br>### log-verossimilhança em 0, 0.01<br>sum(dnorm(y2, 0, 0.01, log=TRUE))<br><br><br><br>### para distribuições discretas, o <br><br>### log-verossimilhança será sempre negativa<br>y3 <- rpois(100, 3)<br>### log lik<br>sum(dpois(y3, 3, log=TRUE))<br><br><br>att. <br>Elias T. Krainski<br>><br>>Ei, Elias.<br>>Agora sim, consegui rodar, mas o DIC está dando negativo, com qualquer das estratégias utilizadas (gaussiana, laplace ou laplace simplificada)...<br>>Percebeu
isso?<br>> <br>>Abraço, Taynãna<br>><br>><br>>Taynãna, fiz <br>> inla.upgrade(testing=TRUE)<br>>atualizou para a versao de 5 de novembro<br>>e consegui executar todo o script SwissRainfall.R sem problemas...<br>><br>><br>><br>> Att.<br>>Elias T. Krainski<br>><br>>>Oi, Elias.<br>>><br>>><br>>>O Havard Rue respondeu dizendo p tentar o comando inla.upgrade(testing=TRUE), mas não deu certo.<br>>>A versão que estou usando é de 5 de novembro.<br>>><br>>><br>>>Abraço, Taynãna<br>>><br>>>Oi Taynãna, <br>>><br>>>qual a versão do seu INLA.<br>>> inla.version()<br>>>Eu consegui executar sem problemas na versão de 18 de agosto<br>>><br>>><br>>><br>>>Elias<br>T. Krainski<br>>><br>>>><br>>>><br>>>>Pessoal,<br>>>><br>>>>Estou
tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo <br>>>><br>>>>de Blangiardo et al (2012). O código em R está disponível e o exemplo 3<br>(SwissRainfall.R).<br>>>><br>>>>Ao rodar o modelo aparece a seguinte mensagem:<br>>>><br>>>>> mod = inla(formula,<br>>>>+ data=inla.stack.data(stack, spde=spde),family="gaussian",<br>>>>+ control.predictor=list(A=inla.stack.A(stack), compute=TRUE),<br>>>>+ <br>>>control.compute=list(cpo=TRUE, dic=TRUE),<br>>>>+ keep=FALSE,
verbose=TRUE)<br>>>>Erro em .local(x, i, j, ..., value) : <br>>>> not-yet-implemented 'Matrix[<-' method<br>>>>Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local<br>>>><br>>>>Alguém tem idéia do quê possa estar<br>acontecendo?<br>>>><br>>>>Att,<br>>>>Taynãna<br>>>>Pesquisadora em Saúde Pública - ENSP/FIOCRUZ <br>><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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><br> </div>
</div> </div></body></html>