posterior_predict function

Generic function for drawing from the posterior predictive distribution