ogrisel, (edited ) Today with @dholzmueller we explored the possibility to reduce a probabilistic regression problem to a classification problem by binning the target variable and interpolating the conditional CDF estimated by
classifier.predict_proba(X_test).cumsum(axis=1)
to the original continuous range.Here is a notebook with the results of my experiments: