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

data;
list(m=12,n=7,x=c(134,146,104,119,124,161,107,83,113,129,97,123),
              y=c(70,118,101,85,107,132,94),psi=0.002186,omega=0.00235)