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