4.1 Overview
4.1.1 Naive model
Let’s start with a simple example, an exponential smoothing model with \(\alpha=1\). This is called the Naive model: \[\hat{x}_{t} = x_{t-1}\] For the naive model, our forecast is simply the value in the previous time step. For example, a naive forecast of the anchovy landings in 1988 is the anchovy landings in 1987. \[\hat{x}_{1988} = x_{1987}\] This is the same as saying that we put 100% of the ‘weight’ on the most recent value and no weight on any value prior to that. \[\hat{x}_{1988} = 1 \times x_{1987} + 0 \times x_{1986} + 0 \times x_{1985} + \dots\] Past values in the time series have information about the current state, but only the most recent past value.
We can fit this with forecast::Arima()
.
fit.rwf <- forecast::Arima(anchovy87ts, order=c(0,1,0))
fr.rwf <- forecast::forecast(fit.rwf, h=5)
Alternatively we can fit with rwf()
or naive()
which are shortcuts for the above lines. All fit the same model.
fr.rwf <- forecast::rwf(anchovy87ts, h=5)
fr.rwf <- forecast::naive(anchovy87ts, h=5)
A plot of the forecast shows the forecast and the prediction intervals.
plot(fr.rwf)
4.1.2 Exponential smoothing
The naive model is a bit extreme. Often the values prior to the last value also have some information about future states. But the ‘information content’ should decrease the farther in the past that we go.
A smoother is another word for a filter, which in time series parlance means a weighted sum of sequential values in a time series: \[w_1 x_t + w_2 x_{t-1} + w_3 x_{t-2} + \dots\] An exponential smoother is a filter (weighted sum) where the weights decline exponentially (Figure @ref(fig:ets.alpha)).
4.1.3 Exponential smoothing model
A simple exponential smoothing model is like the naive model that just uses the last value to make the forecast, but instead of only using the last value it will use values farther in the past also. The weighting function falls off exponentially as shown above. Our goal when fitting an exponential smoothing model is to find the the \(\alpha\), which determines the shape of the weighting function (Figure @ref(fig:ets.alpha2)), that minimizes the forecast errors.