model; {
theta0 ~ dnorm( 0.0,1.0E-6)
theta1 ~ dnorm( 0.0,1.0E-6)
for( i in 1 : I ) {
mu[i] <- theta0+theta1 * x[i]
}
for( i in 1 : I ) {
logit(p[i]) <- mu[i]
}
for( i in 1 : I ) {
y[i] ~ dbin(p[i],n[i])
}
}
and we can read in the data using
list(I=6,x=c(28,29,30,31,32,33),n=c(6,5,9,9,20,15),
y=c(2,2,7,7,16,14))