model {
  for (j in 1:q) {
    for (i in 1:n[j]) {
      y[i, j] ~ dt(mu[j], sigma.prec[j], nu[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)
    nu[j] <- nu.prec[j] + 1
    nu.prec[j] ~ dexp(1 / 29)
  }
} 