sim_pred_beta_binom function

Simulate from a beta-binomial posterior predictive distribution