model
			{
			  for(j in 1 : J) {
			     y[j] ~ dbin(p[j], n[j])
			     logit(p[j]) <- theta[1] + theta[2] * X[j]
			     X[j] ~ dnorm(mu[j], tau)
			     mu[j] <- alpha + beta * Z[j]
			  }
			  theta[1] ~ dnorm(0.0, 0.001)
			  theta[2] ~ dnorm(0.0, 0.001)
			}