model{ b[1] ~ dnorm(0.9, 0.2 ^ -2) b[2] ~ dnorm(1.6, 0.8 ^ -2) for(k in 3:K){ b[k] ~ dlogis(0, 1) } for(i in 1:N){ y[i] ~ dbern(p[i]) logit(p[i]) <- b[1] + b[2] * age[i] + b[3] * year[i] + b[4] * weight[i] } }