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

data;
list(k=10, Y=c(5,1,5,14,3,19,1,1,4,22),
      t=c(94.320,15.720,62.880,125.760,5.240,
      31.440,1.048,1.048,2.096,10.480))

inits;
list(a=c(1,1,1,1,1,1,1,1,1,1), nu=2,S=2)