process_replications function

Extract posterior predictive distribution