jitter_fcst: Jitter a forecast to account for observation errors.

In andrew-MET/harpPoint: Point verifition for NWP forecasts

Jitter a forecast to account for observation errors.

Description

To account for observation errors in ensemble verification the forecast values can be perturbed by sampling from a specified error distribution. Jittering the forecast is only likely to have an effect on the ensemble spread and rank histograms.

Usage

jitter_fcst(.fcst, jitter_function, fcst_suffix = "_mbr")

Arguments

.fcst	An object of class 'harp_fcst'.
jitter_function	The function to be applied to the forecast values.
fcst_suffix	The suffix that identifies forecast values in the harp_fcst object. This will normally be "_mbr".

Details

Note that the function is applied rowwise so for sampling functions n should always be 1.

Value

An object of class harp_fcst with jittered forecast values.

Examples

jitter_fcst(fcst, function(x) x + rnorm(1, ))

andrew-MET/harpPoint documentation built on Feb. 23, 2023, 1:06 a.m.