model;
{
  for (i in 1:6){
    y[i] ~ dbin(pi[i],n[i])
    logit(pi[i]) <- beta0 + beta1*x[i]
  }
  beta0 ~ dnorm(0,0.001)
  beta1 ~ dgamma(0.001,0.001)
}

data;
list(y=c(0,9,21,47,60,63),
   n=c(70,70,70,70,70,70).
   x=c(0,1,2,3,4,5))

inits;
list(beta0=0,beta1=1)