SpecsVerification (version 0.4-1)

FairRps: Fair Ranked Probability Score for ensemble forecasts.

Description

Calculate the Fair Ranked Probability Score for ensemble forecasts and their corresponding observations

Usage

FairRps(ens, obs)

Arguments

ens

N*K matrix. ens[i,j] is the number of ensemble members that predict category j at time i.

obs

N*K matrix. obs[i,j] = 1 if category j is observed at time i, 0 otherwise.

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(c(1,2,3, 2,2,2), byrow=TRUE, nrow=2)
  obs <- matrix(c(0,0,1, 0,1,0), byrow=TRUE, nrow=2)
  FairRps(ens,obs)