Dear all,
I have applied the following inla model. In my model, I have 11 covariates and a GMRF spatial random field. I also included a linear combination predictor (mycomb).
M2 <- inla(f2,
family = "gaussian",
lincomb =mycomb,
data = inla.stack.data(mystack),
control.compute = list(dic = TRUE, waic = TRUE,
openmp.strategy="huge",
config = TRUE),
control.predictor = list(
A = inla.stack.A(mystack)),
control.inla = list (lincomb.derived.only = FALSE,
strategy = "gaussian", int.strategy = "eb"),
verbose=TRUE
)
Now I would like to sample from the joint posterior distribution:
mysample <- inla.posterior.sample(n=100, result=M2)
However, I have got the following error:
Error in dimnames(x) <- dn :
length of 'dimnames' [1] not equal to array extent
Can someone help me with where i should check the error is?
Thanks
Chun