predict_scores function

Compute predictive Bayesian scores