Dear INLA community,
I want to validate my model by simulating data from it and compare it with the observations. In inla.doc("gamma") documents I see that the gamma distribution needs the follwing input variables: y = rgamma(n, shape = a, scale = b), whereby
a = prec.par * prec.scale
b = mu / (prec.par * prec.scale)
prec.scale = runif(n, min = 0.5, max = 2)
prec.par = 1.2
I use the "Precision parameter for the Gamma observations" (one of the available hyperpars of the model) for prec.par, but I cannot find what I should use for prec.scale? Is this something specific per model simulation or something else?
Regards, Bram