model; { r ~ dbin(pi,N) pi ~ dbeta(1,1) pi2 ~ dbeta(1,1) diff <- pi2 - pi ratio <- pi2 / pi r2 ~ dbin(pi2,N2) } data list(N=14,r=4,N2=12,r2=6)