SpecsVerification (version 0.4-1)

FairBrier: Fair Brier Score for ensemble forecasts.

Description

Calculate the Fair Brier Score for ensemble forecasts and their corresponding observations

Usage

FairBrier(ens, obs, tau)

Arguments

ens

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

obs

vector of length N. The corresponding verifying observations.

tau

The threshold value whose exceedance defines the binary event. If ensemble members and observations are coded as 0's and 1's already, setting tau=0.5 (the default) will produce the desired result.

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(FairBrier(ens, obs, tau=1))