221. A-N-A

SES + additive seasonality

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

Given

  • Smoothing parameters: 𝛼=0.5, 𝛾=0.2
  • Initial states: 𝑙0=12, (𝑠3,𝑠2,𝑠1,𝑠0)=(2,0,3,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+𝛼𝜀𝑡𝑠𝑡=𝑠𝑡𝑚+𝛾𝜀𝑡

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+2=14𝜀1=𝑥1𝜇1=1214=2𝑙1=12+0.5(2)=11𝑠1=2+0.2(2)=1.6

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+𝛼𝜀𝑡𝑠𝑡=𝑠𝑡𝑚+𝛾𝜀𝑡
112142111.6
21011110.50.2
387.50.510.752.9
41111.750.7510.3750.85
51411.9752.02511.38752.005
61211.18750.812511.79370.0375
798.89370.106311.84692.8787
81312.69690.303111.99840.9106
91614.00341.996612.99672.4043
101412.95921.040813.51710.1707
111110.63840.361613.69792.8064
121514.60860.391413.89370.9889
131816.2981.70214.74472.7447
141614.91531.084715.2870.3876
151312.48060.519415.54672.7025
161716.53560.464415.77891.0818