227. 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.80.5+2=14.4𝜀1=𝑥1𝜇1=1214.4=2.4𝑙1=12+0.80.5+0.5(2.4)=11.2𝑏1=0.80.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.42.411.20.561.52
21010.7520.75210.3760.74880.1504
386.7771.22310.38850.10982.7554
41111.30060.300610.15030.20810.9399
51411.50382.496211.23190.8322.0192
61211.74710.252912.02390.76670.0998
799.8820.88212.19640.26062.9318
81313.34470.344712.23250.07060.8709
91614.30821.691813.13490.73322.3576
101413.62160.378413.91060.73790.0241
111111.56920.569214.21640.36273.0456
121515.37740.377414.31780.13920.7954
131816.78671.213315.03580.59662.6003
141615.48890.511115.76860.68170.0781
151313.26840.268416.17980.4383.0993
161717.32570.325716.36740.22020.7303