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