I think I would prefer to use the output of decision_function. This way, we will be able to compute the loss value of linear classifiers without predict_proba support. As indicated in Eq. (4.104), we just need to apply softmax to the output of decision_function.