{

    for ( j in 1:6 ) {

      beta[ j ] ~ dflat( )

  }

    tau.e ~ dgamma( 0.001, 0.001 )

    for ( i in 1:N ) {

      y[ i ] ~ dbin( p[ i ], n[ i ] )
      logit( p[ i ] ) <- x[ i, 1 ] * beta[ 1 ] + x[ i, 2 ] * beta[ 2 ]  +
        x[ i, 3 ] * beta[ 3 ] + x[ i, 4 ] * beta[ 4 ]  + 
        x[ i, 5 ] * beta[ 5 ] + x[ i, 6 ] * beta[ 6 ]  + e[ i ]
      e[ i ] ~ dnorm( 0.0, tau.e )

    }

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

  beta.0 <- inprod( c1[ ], beta[ ] )

  alphaR.1 <- inprod( c2[ ], beta[ ] )

  route.effect <- exp( 2 * alphaR.1 )

  alphaT.1 <- inprod( c3[ ], beta[ ] )

  alphaT.2 <- inprod( c4[ ], beta[ ] )

  alphaRT.1 <- inprod( c5[ ], beta[ ] )

  alphaRT.2 <- inprod( c6[ ], beta[ ] )

  p.33 <- p[ 33 ]

  y33.new ~ dbin( p[ 33 ], n[ 33 ] )

}