Feature Request: The original formulation for Brier score inherently supports multiclass classification [source]. This is currently absent in scikit-learn, which restricts Brier score to binary classification.

I have found it useful to have multiclass brier_score_loss when comparing calibration, and have a custom function for the same. I would be happy to work on this and port the code to the scikit-learn API.