model{
	for(k in 1:(nGroups[3])){#Sex
		for(j in 1:(nGroups[2])){#Disease
			for(i in 1:(nGroups[1])){#Municipality
				Obs[i,j,k]~dpois(mu[i,j,k])
				log(mu[i,j,k])<-log(Expect[i,j,k])+alpha[j,k]+S123[i,j,k]
				RR[i,j,k]<-exp(alpha[j,k]+S123[i,j,k])
				S13[i,j,k]<-inprod2(tS3[,j,i],structure3[,k])
				S123[i,j,k]<-inprod2(S13[i,,k],structure2[,j,k])
			}
			alpha[j,k]~dflat()
			tS3[k,j,1:(nGroups[1])]~car.proper(ceros[],C[],adj[],num[],M[],1,gamma)
			
		}
	}

	for(l in 1:(nGroups[1])){ceros[l]<-0}
	gamma.inf<-min.bound(C[],adj[],num[],M[])
	gamma.sup<-max.bound(C[],adj[],num[],M[])
	gamma~dunif(gamma.inf,gamma.sup)

	#Dimension 2 (Disease): 
	for (i in 1:(nGroups[2])){
		for (j in 1:(nGroups[2])){
			for (k in 1:(nGroups[3])){
				structure2[i,j,k] ~ dnorm(0,prec) 
			}
		}
	}

	#Dimension 3 (Sex): 
	for (i in 1:(nGroups[3])){
		for (j in 1:(nGroups[3])){
			structure3[i,j] ~ dnorm(0,prec)
		}
	}
	prec<-pow(sdstruct,-2)
	sdstruct~dunif(0,100)
}