MSES: Mean Square Error Skill score

In SeaVal: Validation of Seasonal Weather Forecasts

Mean Square Error Skill score

Description

Function for taking MSE skill scores of ensemble forecasts stored in long data tables. Can also handle point forecasts. The skill score needs a climatological forecast as reference. This is so far always based on the leave-one-year-out climatology.

Usage

MSES(dt, f, o = "obs", by = by_cols_ens_fc_score(), pool = c("year"), ...)

Arguments

dt	Data table containing the predictions.
f	column name of the prediction.
o	column name of the observations.
by	column names of grouping variables, all of which need to be columns in dt. A separate MSE is computed for each value of the grouping variables. Default is to group by all instances of month, season, lon, lat, system and lead_time that are columns in dt.
pool	column name(s) for the variable(s) along which is averaged. Needs to contain 'year' since the reference climatology forecast is leave-one-year-out.
...	passed on to MSE

Value

A data table with the scores

Examples

dt = data.table(fc = 1:4,obs = c(4,4,7,7),member = c(1,2,1,2),year = c(1999,1999,2000,2000))
MSES(dt,f = 'fc')

SeaVal documentation built on June 24, 2024, 5:15 p.m.