combineSLP: Spread-Adjusted Linear Pool (SLP) of two Normals

Description Usage Arguments Details Value Author(s) References Examples

Description

Computes the weight and scale parameters of the SLP of two normals from a rolling training period.

Usage

1
2
3
combineSLP(x, par_one = list(mean = NULL, sd = NULL),
  par_two = par_one, train = 90, weight_grid = seq(0, 1, 0.1),
  scale_grid = seq(0.6, 1.4, 0.1))

Arguments

x

A vector of observations of a weather quantity.

par_one

A list with two vectors mean and sd of the same length as x respectively, providing parameters of the normal predictive distribution.

par_two

A list with two vectors mean and sd of the same length as x respectively, providing parameters of the second normal predictive distribution.

train

The length of the training period.

weight_grid

The possible values of the SLP weight.

scale_grid

The possible values of the SLP scale parameter.

Details

For each forecast date (all dates except for the first train dates), the optimal combination of SLP weight and scale parameters is found on a grid determined by weight_grid and scale_grid, such that the corresponding (predictive distribution, observation)-pair minimizes the average CRPS with respect to the rolling training period of length train.

Value

A data frame with the observation, predictive mean and standard deviation of the two normals, and SLP weight and scale parameter for the forecast period.

Author(s)

J. Gross, A. Moeller.

References

Gneiting T., Ranjan R. 2013. Combining predictive distributions. Electronic Journal of Statistics, 7, 1747–1782.

Examples

1
combineSLP(16.1, list(mean = 15, sd = 0.8), list(mean = 18, sd = 1), 0)

JuGross/ensAR documentation built on May 10, 2019, 8:23 a.m.