242. A-Ad-A

Damped additive Holt-Winters

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

Given

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

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

𝑥̂𝑡+ℎ|𝑡=𝑙𝑡+(𝜑+𝜑2+…+𝜑ℎ)𝑏𝑡+𝑠𝑡+ℎ−𝑚ℎ+

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+0.8⋅0.5+2=14.4𝜀1=𝑥1−𝜇1=12−14.4=−2.4𝑙1=12+0.8⋅0.5+0.5⋅(−2.4)=11.2𝑏1=0.8⋅0.5+0.4⋅(−2.4)=−0.56𝑠1=2+0.2⋅(−2.4)=1.52

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+𝛽𝜀𝑡𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜀𝑡
11214.4−2.411.2−0.561.52
21010.752−0.75210.376−0.7488−0.1504
386.7771.22310.3885−0.1098−2.7554
41111.3006−0.300610.1503−0.20810.9399
51411.50382.496211.23190.8322.0192
61211.74710.252912.02390.7667−0.0998
799.882−0.88212.19640.2606−2.9318
81313.3447−0.344712.23250.07060.8709
91614.30821.691813.13490.73322.3576
101413.62160.378413.91060.7379−0.0241
111111.5692−0.569214.21640.3627−3.0456
121515.3774−0.377414.31780.13920.7954
131816.78671.213315.03580.59662.6003
141615.48890.511115.76860.68170.0781
151313.2684−0.268416.17980.438−3.0993
161717.3257−0.325716.36740.22020.7303