model { 
	for (i in 1:T) {
		totals[i] ~ dbin(p[i],N[i]);
		logit(p[i]) <- z[i]; 
		z[i] ~ dnorm(mu,tau);
		}
	mu ~ dnorm(0,0.001); # Prior distribution for mu 
	sigma~dunif(0,20);
	tau <-pow(sigma,-2);	
pr<-p[];
	#tau ~ dgamma(0.001,0.001);	# Prior distribution for tau	
	#y ~ dnorm(mu, tau) ; # Predictive distribution for rate
	#sigma <- 1/sqrt(tau);	# SD on the logit scale 
	#w <- exp(y)/(1+exp(y)); # Predictive dist back on p-scale
}