logitp[i] <- alpha + beta * voc[i] p[i] <- exp(logitp[i]) / (1 + exp(logitp[i])) observed[i] ~ dbin(p[i], 1)
For next edition, rewrite all JAGS code so that the order/format matches how we write down equations describing our models. E.g.:
response variable ~ statistical distribution(parameters) transformation(parameters) <- linear predictor