model{
 for(i in 1:23){
   y[i] ~ dbin(p[i],1)
   logit(p[i]) <- beta[1] + beta[2]*tau[i]
   }
   beta[1] ~ dnorm(10,.001)
   beta[2] ~ dnorm(0,.004)
}