241. 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.80.5=12.4𝜀1=(𝑥1𝜇1)/𝜇1=(1212.4)/12.4=0.0323𝑙1=(12+0.80.5)(1+0.5(0.0323))=12.2𝑏1=0.80.5+0.4(12+0.80.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.40.032312.20.24
21012.3920.19311.1960.7648
3810.58420.24429.29211.6455
4117.97570.37929.48780.1067
5149.40250.48911.70121.7537
61213.10420.084312.55210.9613
7913.32110.324411.16050.9594
81310.3930.250811.69650.2753
91611.91670.342713.95841.8535
101415.44120.093314.72060.9063
111115.44570.287813.22281.0532
121512.38030.211613.69010.2053
131813.85440.299215.92721.8225
141617.38520.079716.69260.9039
151317.41570.253515.20791.0432
161714.37330.182715.68670.2161