model{
   eta ~ dbeta(0.01,0.01)
   p[1] <- (2+eta)/4
   p[2] <- (1-eta)/4
   p[3] <- (1-eta)/4
   p[4] <- eta/4
   x[1,1:4] ~ dmulti(p[],N)
}

data;
list(N=197,x=c(125,18,20,32))