model { for(i in 1 : N) { scores[i] ~ dnorm(theta, tau) } theta ~ dnorm(mu, xi) }