model;
{
   mu ~ dnorm( 0.0,1.0E-6)
   psi ~ dgamma(0.01,0.01)
   for( i in 1 : n ) {
      theta[i] ~ dnorm(mu,psi)
   }
   for( i in 1 : n ) {
      x[i] ~ dnorm(theta[i],1)
   }
}

data;
list(n=18,x=c(9.186,8.881,8.571,8.258,7.938,7.938,7.613,7.280,6.938,6.938,
                      6.586,6.586,6.586,6.586,6.586,6.220,5.839,5.439))