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)
    }
} 