predict_history function

Computes predicted cumulative history using posterior predictive distribution.