model
{
	for (j in 1:N.site) {
		Y[j] ~ dpois(mean[j])
		log(mean[j]) <- beta + r[j]
		r[j] ~ dnorm(0.0, tau)
	}
	beta ~ dnorm(0, 1.0E-4)
	tau <- 1 / (s * s)
	s ~ dunif(0, 1.0E+4)
}