model { for( i in 1 : N ) { Y[i] ~ dnorm(mu, prec) } mu ~ dunif (-1000, 1000) prec ~ dunif (0, 1000) s2 <- 1/prec }