model;
{
  for (i in 1:8){
    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(6,13,18,28,52,53,61,60),
   n=c(59,60,62,56,63,59,62,60),
  x=c(1.6907,1.7242,1.7552,1.7842,1.8113,1.8369,1.8610,1.8839))
      
inits;
list(beta0=0,beta1=1)