Efficient computation of predictive probabilities in probit models via expectation propagation