SpecsVerification (version 0.4-1)

EnsRpss: Calculate ranked probability skill score for competing ensemble forecasts.

Description

Calculate the ranked probability skill score (RPSS) between two competing ensemble forecasts for the same observation. An approximate standard deviation of the skill score is provided. The higher the skill score, the higher the improvement of `ens` over `ens.ref`.

Usage

EnsRpss(ens, ens.ref, obs)

Arguments

ens

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

ens.ref

N*K matrix, similar to ens

obs

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

Value

A list with the following elements:"rpss": The value of the skill score."sigma.rpss": The standard deviation of the skill score, approximated by propagation of uncertainty.

Examples

  # Example:
  ens <- matrix(c(1,2,3, 1,3,2), byrow=TRUE, nrow=2)
  ens.ref <- matrix(c(2,2,2, 2,2,2), byrow=TRUE, nrow=2)
  obs <- matrix(c(0,0,1, 0,1,0), byrow=TRUE, nrow=2)
  EnsRpss(ens, ens.ref, obs)