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