Mostly copy-paste from scoringRules::crps_edf
. Adjusted to the data table format, where the observation is a vector of the same length as the ensemble forecast,
but is just repeated (which is why only y[1]
) is used.
crps_aux_esc(y, dat)
vector of length m with m identical entries, the observation
vector of length m containing the m ensemble forecasts