model { prec ~ dgamma(0.001, 0.001) for (i in 1:N) { y[i] <- points[i] y[i] ~ dnorm(mu[i], prec) I(0,) mu[i] <- hca*location[i] + off[o1[i]] + off[o2[i]] + off[o3[i]] + off[o4[i]] + off[o5[i]] + def[d1[i]] + def[d2[i]] + def[d3[i]] + def[d4[i]] + def[d5[i]] } hca ~ dnorm(0,1.0E-6) I(0,) for (i in 1:P) { off[i] ~ dnorm(0,1.0E-6) I(0,) def[i] ~ dnorm(0,1.0E-6) I(0,) } }