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)
}