model {
  for (j in 1:q) {
    for (i in 1:n[j]) {
      y[i, j] ~ dnorm(mu[j], sigma.prec[j])
    }
  }
  # Priors 
  for (j in 1:q) {
    mu[j] ~ dnorm(0, 0.00001)
    sigma[j] <- 1 / sqrt(sigma.prec[j])
    sigma.prec[j] ~ dgamma(0.0001,0.0001)
  }
} 