model;
{
   for( i in 1 : m ) {
      x[i] ~ dnorm(lambda,psi)
   }
   lambda ~ dnorm( 0.0,1.0E-6)
   mu ~ dnorm( 0.0,1.0E-6)
   delta <- lambda-mu
   psi ~ dgamma(0.001,0.001)
   omega ~ dgamma(0.001,0.001)
   for( j in 1 : n ) {
      y[j] ~ dnorm(mu,omega)
   }
}

#   Bartlett quoted only means and s.d.s, so data below is
#   random normal deviates given a linear transofrmation 
#   to have required sample means and variances.

data;
list(m=7,n=7,x=c(430.1,465.9,428.6,450.4,460.1,434.6,470.5),
             y=c(399.5,364.7,451.9,438.5,360.7,427.3,417.7))

data;
list(lambda=0,mu=0,psi=1,omega=1)