model{
for(i in 1:n){
y[i] ~ dbin(p[i],1)
p[i] <- phi(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)}
