model;
{
   for( k in 1 : K ) {
      for( m in 1 : 2 ) {
         logit(tau[k , m]) <- gamma[k , m]
         gamma[k , m] ~ dnorm( 0.0,0.45)
      } 
   }

   alpha[1] ~ dnorm(0,0.45)
   pi[1] <- exp(alpha[1]) / (1 + exp(alpha[1]))
   pi[2] <- 1-pi[1]

for(i in 1:N) {
	eta[i] ~ dcat(pi[1:2])
	for(k in 1:K) {
	    logit(p[i , k]) <- gamma[k , eta[i]]
         y[k,i] ~ dbin(p[i , k],1)
      }

   }
}