# SIMlogit02.bug model{ for(i in 1:n){ xb[i] <- alpha + beta[1] * x[i,1] + beta[2] * x[i,2] + beta[3] * x[i,3] theta[i] <- 1 / (1 + exp(-xb[i])) y[i] ~ dbern(theta[i]) } alpha ~ dnorm(0, .001) for(i in 1:k){ beta[i] ~ dnorm(0, .001) } }