model {
for (i in 1:11) {
y[i] ~ dbin(theta[i], n[i])
logit(theta[i]) <- mu + eta[i]*xi
eta[i] ~ dnorm(0, inv.omega.eta.squared)
}
inv.omega.eta.squared <- 1/pow(omega.eta, 2)
omega.eta <- omega/abs(xi)
xi ~ dnorm(0, 1)
omega ~ dunif(0, 100)
mu ~ dunif(-100, 100)
}
Data:
list(y=c(25,24,23,25,42,24,53,26,25,58,31),
n=c(187,323,122,164,405,239,482,195,177,581,301))
Inits:
list(eta = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
mu = -2, omega = 0.1, xi = 1)