sim_pred_gamma_pois function

Simulate from a gamma-Poisson posterior predictive distribution