for(i in (ns.a + 1):(ns.a + ns.r2)) {
	m[i,2] ~ dnorm(delta[i,2], prec[i,2])
	prec[i,2] <- pow(e[i,2], -2)
}
