4.5 State space form of ETS
One of the main advantages of the ETS model is its state space form, which gives it the flexibility. We would need to revert to linear algebra in this section in order to understand how any ETS model can be presented in a compact state space form.
Hyndman et al. (2008) use the following general formulation of the model with the first equation called “measurement equation” and the second one “transition equation”: \[\begin{equation} \begin{aligned} {y}_{t} = &w(\mathbf{v}_{t-1}) + r(\mathbf{v}_{t-1}) \epsilon_t \\ \mathbf{v}_{t} = &f(\mathbf{v}_{t-1}) + g(\mathbf{v}_{t-1}) \epsilon_t \end{aligned}, \tag{4.21} \end{equation}\] where \(\mathbf{v}_t\) is the state vector, containing the components of series (level, trend and seasonal), \(w(\cdot)\) is the measurement,\(r(\cdot)\) is the error, \(f(\cdot)\) is the transition and \(g(\cdot)\) is the persistence functions. Depending on the types of components these functions can have different values:
- Depending on the types of trend and seasonality \(w(\mathbf{v}_{t-1})\) will be equal either to the addition or multiplication of components. The special cases were presented in tables 3.1 and 3.2 in the ETS Taxonomy section. For example, in case of ETS(M,M,M) it is: \(w(\mathbf{v}_{t-1}) = l_{t-1} b_{t-1} s_{t-m}\);
- If the error is additive, then \(r(\mathbf{v}_{t-1})=1\), otherwise (in case of multiplicative error) it is \(r(\mathbf{v}_{t-1})=w(\mathbf{v}_{t-1})\);
- The transition function \(f(\cdot)\) will produce values depending on the types of trend and seasonality and will correspond to the first parts in the tables 3.1 and 3.2 of the transition equations (dropping the error term). This function records how components interact with each other and how they change from one observation to another (thus the term “transition”). An example is the ETS(M,M,M) model, for which the transition function will produce three values: \(l_{t-1}b_{t-1}\), \(b_{t-1}\) and \(s_{t-m}\) respectively for the level, trend and seasonal components. So, if we drop the persistence function \(g(\cdot)\) and the error term \(\epsilon_t\) for a moment, the second equation in (4.21) will be: \[\begin{equation} \begin{aligned} {l}_{t} = &l_{t-1}b_{t-1} \\ b_t = &b_{t-1} \\ s_t = &s_{t-m} \end{aligned}, \tag{4.22} \end{equation}\]
- Finally, the persistence function will differ from one model to another, but in some special cases it can either be: \(g(\mathbf{v}_{t-1})=\mathbf{g}\), if the error term is additive and \(g(\mathbf{v}_{t-1})=f(\mathbf{v}_{t-1})\mathbf{g}\) if it is multiplicative. \(\mathbf{g}\) is the vector of smoothing parameters, called in the ETS context the “persistence vector.” An example of persistence function is the ETS(M,M,M) model, for which it is: \(l_{t-1}b_{t-1}\alpha\), \(b_{t-1}\beta\) and \(s_{t-m}\gamma\) respectively for the level, trend and seasonal components. Uniting this with the transition function (4.22) we get the equation from the table 3.2: \[\begin{equation} \begin{aligned} {l}_{t} = &l_{t-1}b_{t-1} (1+\alpha\epsilon_t)\\ b_t = &b_{t-1} (1+\beta\epsilon_t)\\ s_t = &s_{t-m} (1+\gamma\epsilon_t) \end{aligned}, \tag{4.23} \end{equation}\]
The compact form (4.21) is thus comfortable to work with and underlies all the 30 ETS models discussed in the sections 5.1 and 5.2. Unfortunately, they cannot be used directly for the derivation of conditional values, so they are needed just for the general understanding of ETS.