model 
{
for (i in 1:N) {
   s[i] ~ dbin(g[i], M[i])   
   logit(g[i]) <- c[prov[i]]*log(encp[i]-1)+b[regimenum[i]]*incumb[i] +rho[regimenum[i]]*lnvdp[i] 
               }

for (k in 1:K) {
            rho[k]~dexp(mu.rho[k])
		log(mu.rho[k]) <- mu.rho.n
		b[k]~dnorm(mu.b[k],tau.b)
		mu.b[k] <- r.b*refelec[k]
	         }

for (j in 1:24) {
		c[j] ~ dnorm(0,tau.c)
	} 
   r.b ~ dnorm(0,.001)
   mu.rho.n ~ dnorm(0,.001)
   tau.b <- pow(sigma.b,-2)
   sigma.b ~ dunif(0,1000)
   tau.c <- pow(sigma.c,-2)
   sigma.c ~ dunif(0,1000)

}