Forecasting time series with complex seasonal patterns using exponential smoothing

Alysha M De Liv­era, Rob J Hyndman and Ralph D Snyder

Journal of the Amer­ican Stat­ist­ical Asso­ci­ation (2011) 106(496), 1513–1527.

Abstract
A new innov­a­tions state space mod­el­ing frame­work, incor­por­at­ing Box-​​Cox trans­form­a­tions, Four­ier series with time vary­ing coef­fi­cients and ARMA error cor­rec­tion, is intro­duced for fore­cast­ing com­plex sea­sonal time series that can­not be handled using exist­ing fore­cast­ing mod­els. Such com­plex time series include time series with mul­tiple sea­sonal peri­ods, high fre­quency sea­son­al­ity, non-​​integer sea­son­al­ity and dual-​​calendar effects. Our new mod­el­ling frame­work provides an altern­at­ive to exist­ing expo­nen­tial smooth­ing mod­els, and is shown to have many advant­ages. The meth­ods for ini­tial­iz­a­tion and estim­a­tion, includ­ing like­li­hood eval­u­ation, are presen­ted, and ana­lyt­ical expres­sions for point fore­casts and inter­val pre­dic­tions under the assump­tion of Gaus­sian errors are derived, lead­ing to a simple, com­pre­hens­ible approach to fore­cast­ing com­plex sea­sonal time series. Our tri­go­no­met­ric for­mu­la­tion is also presen­ted as a means of decom­pos­ing com­plex sea­sonal time series, which can­not be decom­posed using any of the exist­ing decom­pos­i­tion meth­ods. The approach is use­ful in a broad range of applic­a­tions, and we illus­trate its ver­sat­il­ity in three empir­ical stud­ies where it demon­strates excel­lent fore­cast­ing per­form­ance over a range of pre­dic­tion hori­zons. In addi­tion, we show that our tri­go­no­met­ric decom­pos­i­tion leads to the iden­ti­fic­a­tion and extrac­tion of sea­sonal com­pon­ents, which are oth­er­wise not appar­ent in the time series plot itself.

Keywords: expo­nen­tial smooth­ing, Four­ier series, pre­dic­tion inter­vals, sea­son­al­ity, state space mod­els, time series decomposition.

Pre-​​publication work­ing paper

Pub­lished paper

Data

Errata

  • p.1517. \bm{\theta}=(\theta_1,\theta_2,\dots,\theta_q).