237. A-N-M

SES + multiplicative seasonality

ETS(𝐴,𝑁,𝑀)𝑥𝑡=𝑙𝑡−1𝑠𝑡−𝑚+𝜀𝑡𝑙𝑡=𝑙𝑡−1+𝛼𝜀𝑡/𝑠𝑡−𝑚𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜀𝑡/𝑙𝑡−1𝑥̂𝑡+ℎ|𝑡=𝑙𝑡𝑠𝑡+ℎ−𝑚ℎ+
Example: ETS(𝐴,𝑁,𝑀)

Given

  • Smoothing parameters: 𝛼=0.5, 𝛾=0.2
  • Initial states: 𝑙0=12, (𝑠−3,𝑠−2,𝑠−1,𝑠0)=(1.2,1,0.8,1), seasonal period 𝑚=4
  • Data:
𝑡12345678910111213141516
x𝑡121081114129131614111518161317

Step 1 — formula

Observation:

𝑥𝑡=𝑙𝑡−1𝑠𝑡−𝑚+𝜀𝑡

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

𝜇𝑡=𝑙𝑡−1𝑠𝑡−𝑚

Innovation:

𝜀𝑡=𝑥𝑡−𝜇𝑡

State updates:

𝑙𝑡=𝑙𝑡−1+𝛼𝜀𝑡/𝑠𝑡−𝑚𝑠𝑡=𝑠𝑡−𝑚+𝛾𝜀𝑡/𝑙𝑡−1

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

𝑥̂𝑡+ℎ|𝑡=𝑙𝑡𝑠𝑡+ℎ−𝑚ℎ+

where 𝑚ℎ+=((ℎ−1)mod𝑚)+1 picks the right seasonal slot for the period ℎ steps ahead (cycles through 1,2,…,𝑚).

Step 2 — apply at 𝑡=1

𝜇1=12⋅1.2=14.4𝜀1=𝑥1−𝜇1=12−14.4=−2.4𝑙1=12+0.5⋅(−2.4)/1.2=11𝑠1=1.2+0.2⋅(−2.4)/12=1.16

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
11214.4−2.4111.16
21011−110.50.9818
388.4−0.410.250.7924
41110.250.7510.6251.0146
51412.3251.67511.3471.1915
61211.14070.859311.78460.997
799.3379−0.337911.57140.7866
81311.74071.259312.19191.0364
91614.52711.472912.811.2157
101412.77111.228913.42631.0162
111110.56180.438213.70490.7932
121514.20370.796314.0891.048
131817.12790.872114.44771.2281
141614.6811.31915.09671.0344
151311.97431.025715.74330.8068
161716.49930.500715.98221.0544