model
{
for (i in 1:N) {
	for (j in 1:m) {
		x[i,j] ~ dpois( theta[i] )
	}
	theta[i] ~ dgamma( alpha, beta )
}

alpha ~ dexp( 0.1 )
beta ~ dgamma(0.1, 1)
}


inits
list( alpha = 1, beta = 1)
list( alpha = 5, beta = 2)
list( alpha = 2, beta = 5)


#data
list( N = 4, m = 3 )

#more data
x[,1] x[,2] x[,3] 
1        1        0
3        2        4
4        0        1
1        7        3
END