def _crps_gaussian(forecast, mu, sig, **metric_kwargs):
"""CRPS assuming a gaussian distribution. Helper function for CRPSS.
See also:
* properscoring.crps_gaussian
* xskillscore.crps_gaussian
"""
return crps_gaussian(forecast, mu, sig)
def _crps_gaussian(forecast, mu, sig, dim='svd', comparison=None):
return crps_gaussian(forecast, mu, sig).mean(dim)
def crpss(o, f_prob, dim=None):
return 1 - xs.crps_ensemble(o, f_prob, dim=dim) / xs.crps_gaussian(
o, o.mean("time"), o.std("time"), dim=[])