model { for (i in 1:n) { y[i] ~ dbern(p) } p ~ dbeta(1,1) }