model{
for(i in 1:n){
y[i] ~ dbin(p[i],1)
logit(p[i]) <- b1*X[i,1]+b2*X[i,2]+b3*X[i,3]}
b1 ~ dnorm(0, 0.001)
b2 ~ dnorm(0, 0.001)
b3 ~ dnorm(0,0.001)
}
