model;
{
   for( j in 1 : r ) {
      t[j] ~ dnorm( 0.0,1.0E-6)
   }
   for( i in 1 : n ) {
      y[i] ~ dnorm(mu[i],psi)
   }
   for( i in 1 : n ) {
      mu[i] <- t[treat[i]]
   }
   psi ~ dnorm( 0.0,1.0E-6)
   phi <- 1/psi
}

data;
list(n=32,r=7,
   treat=c(1,1,1,1,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,6,7,7,7,7),
   y=c(12, 30, 10, 18, 24, 32, 29, 26,
        9,  9, 16,  4,
       30,  7, 21,  9,
       16, 10, 18, 18,
       18, 24, 12, 19,
       10,  4,  4,  5,
       17,  7, 16, 17)