model
    {
       for( i in 1 : N ) {
          r[i] ~ dbin(p[i],n[i])
          logit(p[i]) <- alpha.star + beta * (x[i] - mean(x[]))
          rhat[i] <- n[i] * p[i]
       }
       alpha <- alpha.star - beta * mean(x[])
       beta ~ dnorm(0.0,0.001)
       alpha.star ~ dnorm(0.0,0.001)   
    }