climpred.metrics._less

Contents

climpred.metrics._less#

climpred.metrics._less(forecast: Dataset, verif: Dataset, dim: str | List[str] | None = None, **metric_kwargs: Any) → Dataset | DataArray[source]#

Logarithmic Ensemble Spread Score.

\[LESS = ln(\frac{variance}{MSE})= ln(\frac{\sigma^2_f}{\sigma^2_o})\]
Parameters:

  • forecast – Forecasts.

  • verif – Verification.

  • dim – The dimension(s) over which to aggregate. Defaults to None, meaning aggregation over all dims other than lead.

Returns:

less – reduced by dimensions dim

Notes

maximum

∞

positive

overdisperive / underconfident

perfect

0

negative

underdisperive / overconfident

minimum

-∞

orientation

None

Example

>>> # better detrend before
>>> from climpred.stats import rm_poly
>>> HindcastEnsemble.map(rm_poly, dim="init_or_time", deg=2).verify(
...     metric="less",
...     comparison="m2o",
...     dim=["member", "init"],
...     alignment="same_verifs",
... ).SST
<xarray.DataArray 'SST' (lead: 10)> Size: 80B
array([ 0.12633668, -0.12707633, -0.26143178, -0.25096534, -0.29267364,
       -0.29057248, -0.43579506, -0.33774944, -0.46008436, -0.61010384])
Coordinates:
  * lead     (lead) int32 40B 1 2 3 4 5 6 7 8 9 10
    skill    <U11 44B 'initialized'
Attributes:
    units:    None

References

Kadow et al. [2016]

Contents