An ensembleData object including ensemble forecasts with
the corresponding verifying observations and their dates.
Missing values (indicated by NA) are allowed.
trainingDays
An integer giving the number of time steps (e.g. days)
in the training period. There is no default.
consecutive
If TRUE then the sequence of dates in the training set are
treated as consecutive, i.e. date gaps are ignored.
dates
The dates for which EMOS forecasting models are desired.
By default, this will be all dates in ensembleData
for which modeling is allowed given the training rule.
control
A list of control values for the fitting functions. The defaults are
given by the function controlMOSnormal.
warmStart
If TRUE, then starting values for parameters in optimization
are set to the estimates of the preceding date's fit.
exchangeable
A numeric or character vector or factor indicating groups of
ensemble members that are exchangeable (indistinguishable).
The modeling will have equal parameters within each group.
The default determines exchangeability from ensembleData.
Details
Given an ensemble forecast of size m: X_1, … , X_m, the
following Gaussian predictive distribution is fit by
ensembleMOSnormal
Y sim mathcal{N}
≤ft( a + b_1X_1 + cdots + b_mX_m , c + dS^2
ight)
B is a vector of fitted regression coefficients: b_1,
… ,b_m. Specifically, a,b_1,… , b_m, c,d are
fitted to optimize
control$scoringRule over the specified training period using
optim with method = "BFGS".
Value
A list with the following output components:
training
A list containing information on the training length and lag and
the number of instances used for training for each modeling date.
a
A vector of fitted EMOS intercept parameters for each date.
B
A matrix of fitted EMOS coefficients for each date.
c,d
Vectors of the fitted variance parameters for each date, see details.
References
T. Gneiting, A. E. Raftery, A. H. Westveld and T. Goldman,
Calibrated probabilistic forecasting using ensemble model output
statistics and minimum CRPS estimation.
Monthly Weather Review 133:1098–1118, 2005.