253. M-A-N

Holt’s linear, multiplicative errors

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

Given

  • Smoothing parameters: 𝛼=0.5, 𝛽=0.4
  • 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):

𝑥̂𝑡+ℎ|𝑡=𝑙𝑡+ℎ𝑏𝑡

where ℎ∈{1,2,3,…} is the forecast horizon (how many steps ahead).

Step 2 — apply at 𝑡=1

𝜇1=12+0.5=12.5𝜀1=(𝑥1−𝜇1)/𝜇1=(12−12.5)/12.5=−0.04𝑙1=(12+0.5)(1+0.5⋅(−0.04))=12.25𝑏1=0.5+0.4(12+0.5)⋅(−0.04)=0.3

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.5−0.0412.250.3
21012.55−0.203211.275−0.72
3810.555−0.24219.2775−1.742
4117.53550.45989.2677−0.3562
5148.91160.57111.45581.6792
61213.135−0.086412.56751.2252
7913.7927−0.347511.3963−0.6919
81310.70450.214411.85220.2263
91612.07860.324714.03931.7949
101415.8342−0.115814.91711.0612
111115.9783−0.311613.4892−0.9301
121512.55910.194413.77950.0463
131813.82580.301915.91291.716
141617.6289−0.092416.81441.0644
151317.8788−0.272915.4394−0.8871
161714.55230.168215.77610.092