model 
{
     for (i in 1:N) {
               sancion[i] ~ dbern(g[i])   
               logit(g[i]) <- balpha[congress[i]] + bpjpresi*pjpresi[i] + binicsen*inicsen[i]+ binicses*inicses[i]+ binicsen.pjpres*(inicsen[i]*pjpresi[i]) + bhoney*honeymoo[i] + bcomdip*comdip[i] + bcomsen*comsen[i] + bdistrib[distrib[i]] + breach[reach[i]]  
     		   e.y[i] <- sancion[i] - g[i]
            }

	for (k in 1:K) {
            balpha[k]~dnorm(mu.alpha[k], tau[15])
            mu.alpha[k] <- cincumb*incumb[k]+cnbills*nbills[k]+cippoleco*ippoleco[k]+ cratdec*ratdec[k]+ cgpprov*log(gpprov[k]/gpnac[k])
				
		a[k] <- xi.a*balpha[k]
		a.hat[k] <- xi.a*mu.alpha[k]
		e.a[k] <- a[k]-a.hat[k]

      	} 


   bcomdip~dnorm(0,tau[1])
   bcomsen~dnorm(0,tau[2])
   bpjpresi~dnorm(0,tau[3])
   binicsen~dnorm(0,tau[4])
   binicsen.pjpres~dnorm(0,tau[5])
   bhoney~dnorm(0,tau[6])
   binicses~dnorm(0,tau[12])

   cincumb~dnorm(0,tau[7])
   cnbills~dnorm(0,tau[8])
   cippoleco~dnorm(0,tau[9])
   cratdec~dnorm(0,tau[10])
   cgpprov~dnorm(0,tau[11])

   xi.a ~ dnorm(0,tau[16])

  for (t in 1:4) {breach[t]~dnorm(0,tau[13])}
  for (j in 1:3) {bdistrib[j]~dnorm(0,tau[14])}

  for (p in 1:16) {
   tau[p] <- pow(sigma[p],-2)
   sigma[p] ~ dunif(0,1000)
         }
 
}