	
	model
	{
		for( i in 1 : batches ) {
			m[i] ~ dnorm(theta, tau.btw)
			for( j in 1 : samples ) {
				y[i , j] ~ dnorm(m[i], tau.with)
			}
		}	
		sigma2.with <- 1 / tau.with
		sigma2.btw <- 1 / tau.btw
		tau.with ~ dgamma(0.001, 0.001)
		tau.btw ~ dgamma(0.001, 0.001)
		theta ~ dnorm(0.0, 1.0E-10)
	}
