Thanks for the response, However I do not understand on 'add the observation-level sampling' part. Assuming the following is my model, and I draw 10 samples, what should follow next?
```
mod <- inla(
formula,
family = "binomial",
Ntrials = numtrials,
data = inla.stack.data(stk.full),
control.family = control.family(link = 'logit'),
control.predictor = list(
compute = TRUE,
link = 1,
A = inla.stack.A(stk.full)
),
control.compute = list(config = TRUE, return.marginals.predictor = TRUE),
control.inla = list(tolerance = 1e-7),
num.threads = 7,
verbose = F
)
nsamp <- 10
#Posterior sampling
ps <- inla.posterior.sample(nsamp, mod, num.threads = 7, parallel.configs = T)
```