model;
{
   for( i in 1 : m ) {
      x[i] ~ dnorm(lambda,psi)
   }
   for( j in 1 : n ) {
      y[j] ~ dnorm(mu,psi)
   }
   mu ~ dnorm( 0.0,1.0E-6)
   nu ~ dnorm( 0.0,1.0E-6)
   psi ~ dgamma(0.01,0.001)
   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))

inits;
list(lambda=0,mu=0,psi=1)