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"), ...)
Value
A data table with the scores
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.