	model
	{
		for(i in 1 : batches) {
			mu[i] ~ dnorm(theta, tau.btw)
			for(j in 1 : samples) {
				y[i , j] ~ dnorm(mu[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)
	}
