{

  mu ~ dnorm( 0.0, 1.0E-4 )
  tau.theta ~ dgamma( 0.0005, 0.000765 )

  for ( i in 1:k ) {

    theta[ i ] ~ dnorm( mu, tau.theta )
    y[ i ] ~ dnorm( theta[ i ], tau.y[ i ] )

  }

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

  positive.effect <- step( mu )

}