{

  theta ~ dnorm( 0.0, 1.0E-6 )
  sigma.H ~ dunif( 0.0, 5.0 )
  sigma.P ~ dunif( 0.0, 50.0 )

  for ( i in 1:10 ) {

    a.H[ i ] ~ dnorm( 0.0, tau.H )
    mu[ i ] <- theta + a.H[ i ]

    for ( j in 1:150 ) {

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

    }

  }

  tau.H <- 1.0 / ( sigma.H * sigma.H )
  tau.P <- 1.0 / ( sigma.P * sigma.P )

}