model
{
  for (i in 1:k) {
    a[i] ~ dchisqr(nu);
    theta[i] <- a[i]/S;
    lambda[i] <- theta[i]*t[i];
    Y[i] ~ dpois(lambda[i]);
    }
    nu ~ dexp(1.0);
    S ~ dchisqr(1.0);
}