245. A-M-A

Multiplicative trend, additive seasonality

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

Given

  • Smoothing parameters: 𝛼=0.5, 𝛽=0.4, 𝛾=0.2
  • Initial states: 𝑙0=12, 𝑏0=1, (𝑠−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𝑏𝑡−1+𝑠𝑡−𝑚

Innovation:

𝜀𝑡=𝑥𝑡−𝜇𝑡

State updates:

𝑙𝑡=𝑙𝑡−1𝑏𝑡−1+𝛼𝜀𝑡𝑏𝑡=𝑏𝑡−1+𝛽𝜀𝑡/𝑙𝑡−1𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜀𝑡

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

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

where ℎ∈{1,2,3,…} is the forecast horizon (how many steps ahead); 𝑚ℎ+=((ℎ−1)mod𝑚)+1 picks the right seasonal slot for the period ℎ steps ahead (cycles through 1,2,…,𝑚).

Step 2 — apply at 𝑡=1

𝜇1=12⋅1+2=14𝜀1=𝑥1−𝜇1=12−14=−2𝑙1=12⋅1+0.5⋅(−2)=11𝑏1=1+0.4⋅(−2)/12=0.9333𝑠1=2+0.2⋅(−2)=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+𝑠𝑡−𝑚𝜀𝑡𝑙𝑡=𝑙𝑡−1𝑏𝑡−1+𝛼𝜀𝑡𝑏𝑡=𝑏𝑡−1+𝛽𝜀𝑡/𝑙𝑡−1𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜀𝑡
11214−2110.93331.6
21010.2667−0.266710.13330.9236−0.0533
386.35951.640510.17980.9884−2.6719
41111.0616−0.061610.03080.9860.9877
51411.49012.509911.1451.08612.102
61212.0509−0.050912.07881.0842−0.0635
7910.4243−1.424312.38411.0371−2.9568
81313.8308−0.830812.42771.01020.8215
91614.65691.343113.22651.05352.3706
101413.87010.129913.99851.0574−0.0375
111111.8452−0.845214.37941.0332−3.1258
121515.6789−0.678914.51791.01440.6857
131817.0970.90315.17791.03922.5512
141615.73590.264115.90551.04620.0153
151313.5145−0.514516.3831.0333−3.2287
161717.6137−0.613716.62111.01830.563