#priors on precisions rather than icc

model
{
for( i in 1 : H ) {
mu[i] ~ dnorm(theta, tau.btw)
for( j in 1 : Y ) {
y[i , j] ~ dnorm(mu[i], tau.with)
}
}
theta ~ dnorm(0.0, 1.0E-2)
# prior for within-variation

sigma.with~dunif(0,20)
sigma.btw~dunif(0,20)
tau.with <-pow(sigma.with,-2)
tau.btw <-pow(sigma.btw,-2)

sigma2.with <-pow(sigma.with,2)
sigma2.btw <-pow(sigma.btw,2)
ICC<- sigma2.btw/(sigma2.btw+sigma2.with)
}