Learn R Programming

SpecsVerification (version 0.4-1)

DressCrpss: Calculate continous ranked probability skill score (CRPSS), comparing two dressed ensemble forecasts.

Description

Calculate the average continuous ranked probability skill score (CRPSS) for a collection of forecast distributions generated from two competing ensemble forecasts dressed with Gaussian kernels. The standard deviation of the skill score approximated by propagation of uncertainty is provided. The higher the skill score, the higher is the improvement of `dressed.ens` over `dressed.ens.ref`.

Usage

DressCrpss(dressed.ens, dressed.ens.ref, obs)

Arguments

dressed.ens
An object of class `dressed.ens`. The ensemble forecast to be analyzed. See ?DressEnsemble for details.
dressed.ens.ref
An object of class `dressed.ens`. The reference forecast, that predicts the same targets as `dressed.ens` and to which `dressed.ens` is compared. See ?DressEnsemble for details.
obs
vector of length nrow(dressed.ens$ens). The same verifying observations for both forecasts.

Value

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

Examples

Run this code
  # Example:
  ens <- matrix(rnorm(100),20,5)
  ens.ref <- ens + 0.2
  obs <- rnorm(20)
  d.ens <- DressEnsemble(ens)
  d.ens.ref <- DressEnsemble(ens.ref)
  DressCrpss(d.ens, d.ens.ref, obs)

Run the code above in your browser using DataLab