for(i in (ns.a + ns.r2 + 1):(ns.a + ns.r2 + ns.rm)) {
	for (k in 1:(na[i]-1)) {
		for (j in 1:(na[i]-1)) {
			Sigma[i,j,k] <- pow(e[i,1], 2) * (1 - equals(j,k)) + pow(e[i,k+1], 2) * equals(j,k)
		}
	}
	Omega[i,1:(na[i]-1),1:(na[i]-1)] <- inverse(Sigma[i,1:(na[i]-1),1:(na[i]-1)])
	m[i,2:na[i]] ~ dmnorm(delta[i,2:na[i]],Omega[i,1:(na[i]-1),1:(na[i]-1)])
}
