SpecsVerification (version 0.4-1)

FairCrps: Fair continuously ranked probability score (CRPS) for ensemble forecasts.

Description

Calculate the fair continuously ranked probability score (CRPS) for ensemble forecasts and their corresponding observations.

Usage

FairCrps(ens, obs)

Arguments

ens

N*K matrix. N cases of a K-member ensemble forecast.

obs

vector of length N. The corresponding verifying observations.

Value

The function returns a vector of length N with the individual score values of the N ensemble forecasts.

References

Ferro, C.A.T. (2013). Fair scores for ensemble forecasts. Q.J.R. Meteorol. Soc. DOI: 10.1002/qj.2270.

Examples

  # Example:
  ens <- matrix(rnorm(500),100,5)
  obs <- rnorm(100)
  mean(FairCrps(ens, obs))