# Jags model file for measuresments within islands analysis based on multi-level radon model in Gelman & Hill p. 350, save to working directory for access model { for (i in 1:N) { y[i] ~ dnorm (y.hat[i], tau.y) y.hat[i] <- a[island[i]] + b0[m[i]] #m is constant for direct/indirect date e.y[i]<- y[i] - y.hat[i] #data-level errors } b0[1] ~dnorm (0, 0.0001) b0[2] ~dnorm (0, 0.0001) b0[3] ~dnorm (0, 0.0001) tau.y <-pow(sigma.y, -2) sigma.y ~dunif (0, 50000) for (j in 1:J) { #J = no. islands a[j] ~ dnorm (a.hat[j], tau.a) a.hat[j] <- mu.a e.a[j] <- a[j] - a.hat[j] #group-level errors } mu.a ~ dnorm (0, 0.0001) tau.a <- pow(sigma.a, -2) sigma.a ~ dunif (0, 50000) #one sigma for all islands }