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])
				for (kk in 1:(nGroups[3])){
					suma[i,j,k,kk]<-inprod2(tS3[kk,,i],structure23[(nGroups[2]*(k-1)+j),((kk-1)*nGroups[2]+1):((kk-1)*nGroups[2]+nGroups[2])])
				}
				S123[i,j,k]<-sum(suma[i,j,k,1:(nGroups[3])])
			}
			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]*nGroups[3])){
		for (j in 1:(nGroups[2]*nGroups[3])){
			structure23[i,j] ~ dnorm(0,prec) 
		}
	}
	prec<-pow(sdstruct,-2)
	sdstruct~dunif(0,100)
}