The mean posterior prediction at each point in X_new should be:
np.mean(pred_samples['f_pred'], axis=0)
Then you should be able to sum them. So all in one step:
np.sum(np.mean(pred_samples['f_pred'], axis=0))
Is that what you meant?
The mean posterior prediction at each point in X_new should be:
np.mean(pred_samples['f_pred'], axis=0)
Then you should be able to sum them. So all in one step:
np.sum(np.mean(pred_samples['f_pred'], axis=0))
Is that what you meant?