model italy
{
 for (i in 1:N) {
 lambda[i] ~ dgamma(alpha,beta)
 y[i] ~ dpois(lambda[i])
 }
alpha ~ dgamma(1,1)
beta ~ dgamma(1,1)
}


list(y=c(7,9,8,7,7,6,6,5,5,7,9,10,8,8,8,7), N=16)

list(alpha=10, beta=10)