climpred.metrics._msess_murphy

Contents

climpred.metrics._msess_murphy#

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

Murphy’s Mean Square Error Skill Score (MSESS).

\[MSESS_{Murphy} = r_{fo}^2 - [\text{conditional bias}]^2 -\ [\frac{\text{(unconditional) bias}}{\sigma_o}]^2,\]

where \(r_{fo}^{2}\) represents the Pearson product-moment correlation coefficient between the forecast and verification data and \(\sigma_{o}\) represents the standard deviation of the verification data over the experimental period. See conditional_bias and unconditional_bias for their respective formulations.

Parameters:

Notes

minimum

-∞

maximum

1.0

perfect

1.0

orientation

positive

See also

References

Example

>>> HindcastEnsemble = HindcastEnsemble.remove_bias(alignment="same_verifs")
>>> HindcastEnsemble.verify(
...     metric="msess_murphy",
...     comparison="e2o",
...     dim="init",
...     alignment="same_verifs",
... )
<xarray.Dataset> Size: 164B
Dimensions:  (lead: 10)
Coordinates:
  * lead     (lead) int32 40B 1 2 3 4 5 6 7 8 9 10
    skill    <U11 44B 'initialized'
Data variables:
    SST      (lead) float64 80B 0.8306 0.8351 0.8295 ... 0.8471 0.813 0.8097
Attributes:
    prediction_skill_software:     climpred https://climpred.readthedocs.io/
    skill_calculated_by_function:  HindcastEnsemble.verify()
    number_of_initializations:     64
    number_of_members:             10
    alignment:                     same_verifs
    metric:                        msess_murphy
    comparison:                    e2o
    dim:                           init
    reference:                     []
Contents