model 
   {
      for (i in 1 : N) {
         y[i] ~ dbin(P[state1[i]], t[i])
         state[i] ~ dbern(theta)
         state1[i] <- state[i] + 1 
         t[i] <- x[i] + y[i] 
         prop[i] <- P[state1[i]] 
      }
      P[1] <- p 
      P[2] <- 0
      logit(p) <- alpha
      alpha ~ dnorm(0,1.0E-4)
      beta <- exp(alpha) 
      logit(theta) <- delta 
      delta ~ dnorm(0, 1.0E-4)
   }