model
{
                for (i in 1 : N) {
                        theta[i] ~ dgamma(alpha, beta)
                        lambda[i] <- theta[i] * t[i]
                        x[i] ~ dpois(lambda[i])
                }
                alpha ~ dexp(1)
                beta ~ dgamma(0.1, 1.0)
}