222. 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=121.2=14.4𝜀1=𝑥1𝜇1=1214.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.42.4111.16
21011110.50.9818
388.40.410.250.7924
41110.250.7510.6251.0146
51412.3251.67511.3471.1915
61211.14070.859311.78460.997
799.33790.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