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