model{
for(i in 1:n){
ones[i] <-1
ones[i] ~ dbern(p[i])
term1[i] <- ct*pow(tau,1/2)*exp(-1/2*tau*pow(y[i]-mu[i],2))
term2[i] <- phi(-mu[i]*pow(tau,1/2))
p[i] <- pow(term1[i],ind[i])*pow(term2[i],1-ind[i])/K
mu[i] <- b1*X[i,1]+b2*X[i,2]+b3*X[i,3]}
b1 ~ dnorm(0, 0.002)
b2 ~ dnorm(0, 0.002)
b3 ~ dnorm(0, 0.002)
tau ~ dgamma(0.001, 0.001)
}
