model; { r ~ dbin(pi,N) pi ~ dbeta(1,1)I( 0.1, 0.6) } data list(N=14,r=4)