251. M-N-A

SES + add. seasonality, mult. errors

ETS(𝑀,𝑁,𝐴)𝑥𝑡=(𝑙𝑡−1+𝑠𝑡−𝑚)(1+𝜀𝑡)𝑙𝑡=𝑙𝑡−1+𝛼(𝑙𝑡−1+𝑠𝑡−𝑚)𝜀𝑡𝑠𝑡=𝑠𝑡−𝑚+𝛾(𝑙𝑡−1+𝑠𝑡−𝑚)𝜀𝑡𝑥̂𝑡+ℎ|𝑡=𝑙𝑡+𝑠𝑡+ℎ−𝑚ℎ+
Example: ETS(𝑀,𝑁,𝐴)

Given

  • Smoothing parameters: 𝛼=0.5, 𝛾=0.2
  • Initial states: 𝑙0=12, (𝑠−3,𝑠−2,𝑠−1,𝑠0)=(2,0,−3,1), seasonal period 𝑚=4
  • Data:
𝑡12345678910111213141516
x𝑡121081114129131614111518161317

Step 1 — formula

Observation:

𝑥𝑡=(𝑙𝑡−1+𝑠𝑡−𝑚)(1+𝜀𝑡)

Conditional mean (one-step-ahead forecast 𝑥̂𝑡|𝑡−1=𝜇𝑡):

𝜇𝑡=𝑙𝑡−1+𝑠𝑡−𝑚

Innovation:

𝜀𝑡=(𝑥𝑡−𝜇𝑡)/𝜇𝑡

State updates:

𝑙𝑡=𝑙𝑡−1+𝛼𝜇𝑡𝜀𝑡𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜇𝑡𝜀𝑡

Forecast ℎ steps ahead from time 𝑡 (using current-period states):

𝑥̂𝑡+ℎ|𝑡=𝑙𝑡+𝑠𝑡+ℎ−𝑚ℎ+

where 𝑚ℎ+=((ℎ−1)mod𝑚)+1 picks the right seasonal slot for the period ℎ steps ahead (cycles through 1,2,…,𝑚).

Step 2 — apply at 𝑡=1

𝜇1=12+2=14𝜀1=(𝑥1−𝜇1)/𝜇1=(12−14)/14=−0.1429𝑙1=12+0.5⋅14⋅(−0.1429)=11𝑠1=2+0.2⋅14⋅(−0.1429)=1.6

Step 3 — iterate

Each column header is the equation that produced its values. Values rounded to 4 decimal places; arithmetic performed at full precision.

𝑡𝑥𝑡𝜇𝑡=𝑙𝑡−1+𝑠𝑡−𝑚𝜀𝑡𝑙𝑡=𝑙𝑡−1+𝛼𝜇𝑡𝜀𝑡𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜇𝑡𝜀𝑡
11214−0.1429111.6
21011−0.090910.5−0.2
387.50.066710.75−2.9
41111.75−0.063810.3750.85
51411.9750.169111.38752.005
61211.18750.072611.7937−0.0375
798.89370.011911.8469−2.8787
81312.69690.023911.99840.9106
91614.00340.142612.99672.4043
101412.95920.080313.51710.1707
111110.63840.03413.6979−2.8064
121514.60860.026813.89370.9889
131816.2980.104414.74472.7447
141614.91530.072715.2870.3876
151312.48060.041615.5467−2.7025
161716.53560.028115.77891.0818