{

  alpha ~ dnorm( 0.0, 1.0E-6 )
  beta ~ dnorm( 0.0, 1.0E-6 )
  tau.theta ~ dgamma( 0.001, 0.001 )

  x.bar <- mean( x[ ] )

  for ( i in 1:n ) {

    mu[ i] <- alpha + beta * ( x[ i ] - x.bar )
    theta[ i ] ~ dnorm( mu[ i ], tau.theta )
    y[ i ] ~ dnorm( theta[ i ], tau[ i ] )

  }

  sigma.theta <- 1.0 / sqrt( tau.theta )

}