256. M-Ad-N

Damped linear, mult. errors

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

Given

  • Smoothing parameters: 𝛼=0.5, 𝛽=0.4, 𝜑=0.8
  • Initial states: 𝑙0=12, 𝑏0=0.5
  • Data:
𝑡12345678910111213141516
x𝑡121081114129131614111518161317

Step 1 — formula

Observation:

𝑥𝑡=(𝑙𝑡−1+𝜑𝑏𝑡−1)(1+𝜀𝑡)

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

𝜇𝑡=𝑙𝑡−1+𝜑𝑏𝑡−1

Innovation:

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

State updates:

𝑙𝑡=(𝑙𝑡−1+𝜑𝑏𝑡−1)(1+𝛼𝜀𝑡)𝑏𝑡=𝜑𝑏𝑡−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).

Step 2 — apply at 𝑡=1

𝜇1=12+0.8⋅0.5=12.4𝜀1=(𝑥1−𝜇1)/𝜇1=(12−12.4)/12.4=−0.0323𝑙1=(12+0.8⋅0.5)(1+0.5⋅(−0.0323))=12.2𝑏1=0.8⋅0.5+0.4(12+0.8⋅0.5)⋅(−0.0323)=0.24

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+𝛽(𝑙𝑡−1+𝜑𝑏𝑡−1)𝜀𝑡
11212.4−0.032312.20.24
21012.392−0.19311.196−0.7648
3810.5842−0.24429.2921−1.6455
4117.97570.37929.4878−0.1067
5149.40250.48911.70121.7537
61213.1042−0.084312.55210.9613
7913.3211−0.324411.1605−0.9594
81310.3930.250811.69650.2753
91611.91670.342713.95841.8535
101415.4412−0.093314.72060.9063
111115.4457−0.287813.2228−1.0532
121512.38030.211613.69010.2053
131813.85440.299215.92721.8225
141617.3852−0.079716.69260.9039
151317.4157−0.253515.2079−1.0432
161714.37330.182715.68670.2161