model{
   a <- y+x[4]
   b <- x[2]+x[3]
   eta ~ dbeta(a,b)
   pi <- eta/(eta+2)
   y ~ dbin(pi,x[1])
}

data;
list(x=c(125,18,20,34))

inits;
list(y=0,eta=0.5)