model{
  for (i in 1:N){
    y[i] ~ dbin(pi[i], n[i])
    pi[i] ~ dbeta(a,b)
    }

a ~ dgamma(0.02, 0.01)
b ~ dgamma(0.02, 0.01)

}