check_obs_against_fcst: Observation error check against forecast

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

Observation error check against forecast

Description

For each stratification the standard deviation of the forecast over all forecast cycles and models (mname), and in the case of an eps forecast, the ensemble members is computed. The difference between the observation and the forecast is expected to be smaller than a number of multiples of the standard deviation. The number of multiples of the standard deviation can be supplied or a default value used depending on the parameter.

Usage

check_obs_against_fcst(
  .fcst,
  parameter,
  num_sd_allowed = NULL,
  stratification = c("SID", "quarter_day")
)

Arguments

.fcst	A data frame of class harp_point_forecast_obs.
parameter	The forecast parameter.
num_sd_allowed	The number of standard deviations of the forecast that the difference between the forecast and the observation must be smaller than.
stratification	The columns to stratify the data by when computing the allowed tolerance. In most cases the column must exist in the input data, but "quarter_day" can be passed to divide the observations into classes of [0, 6), [6, 12), [12, 18) and [18, 24) hour of day. The default behaviour is to stratify by station ("SID") and "quarter_day".

Value

A data frame of class harp_point_forecast_obs with an attribute named removed_cases containing a data frame of the removed cases.

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