Model	{ 
	for (i in 1:Nareas) { 
            for (j in 1:Ntimeperiods) {
            for (k in 1:Ndiseases) {
	             O[i,j,k]  ~ dpois(mu[i,j,k]) 
                  log(mu[i,j,k])