model { for (i in 1:N) { X[i] ~ dnorm(mu,tau); } mu ~ dnorm(0.0, 1.0E-4); tau ~ dgamma(1.0E-3, 1.0E-3); sigma <- 1.0/sqrt(tau); y ~ dnorm(mu,tau) }