	model
	{
		for( i in 1 : N ) {
			r[i] ~ dbin(p[i],n[i])
			logit(p[i]) <- alpha.star + beta * (x[i] - mean(x[]))
			rhat[i] <- n[i] * p[i]
		}
		alpha <- alpha.star - beta * mean(x[])
		beta ~ dnorm(0.0,0.001)
		alpha.star ~ dnorm(0.0,0.001)	
	}
