model{ for (i in 1:n) { y[i] ~ dbern(p[i]) logit(p[i]) <- alpha + beta*x[i] } alpha~dnorm(0,0.001) beta~dnorm(0,0.001) }