{"id":1305,"date":"2017-04-27T22:27:12","date_gmt":"2017-04-27T22:27:12","guid":{"rendered":"http:\/\/kourentzes.com\/forecasting\/?p=1305"},"modified":"2017-07-22T17:08:03","modified_gmt":"2017-07-22T17:08:03","slug":"multiple-temporal-aggregation-the-story-so-far-part-i","status":"publish","type":"post","link":"https:\/\/kourentzes.com\/forecasting\/2017\/04\/27\/multiple-temporal-aggregation-the-story-so-far-part-i\/","title":{"rendered":"Multiple temporal aggregation: the story so far. Part I"},"content":{"rendered":"<p style=\"text-align: justify;\">Over the last years I have been working (with my co-authors!) on the idea of Multiple Temporal Aggregation (MTA) for time series forecasting. A number of papers have been published introducing and developing the idea further, or testing its effectiveness for forecasting.<\/p>\n<p style=\"text-align: justify;\">In this series of blog posts I will try to summarise the progress so far, and highlight ways that you can use it. This first post will summarise the papers so far and give an overview of the main findings. Later posts will focus on explaining how MTA works.<\/p>\n<p style=\"text-align: justify;\">The key points behind MTA are the following:<\/p>\n<ul>\n<li style=\"text-align: justify;\">It is a radically different approach to time series modelling, recognising that the data sampling frequency may not be the best for a given modelling purpose.<\/li>\n<li style=\"text-align: justify;\">A time series is modelled simultaneously at multiple temporal aggregation levels that can be easily generated from the original data. At each level an appropriate model is fit, focusing on the components of the series that are strengthened by temporal aggregation.<\/li>\n<li style=\"text-align: justify;\">If forecasting is the objective, then the produced forecast reconciles the information from all these models. This makes the forecast robust to modelling uncertainty and lessens the importance of model selection.<\/li>\n<li style=\"text-align: justify;\">The resulting forecasts have been shown to be reliable and typically outperform the conventional modelling approach.<\/li>\n<\/ul>\n<p style=\"text-align: justify;\">Table 1 summarises our contributions on MTA so far (follow the links to access the papers). We have also released two R packages that implement MTA: <a href=\"https:\/\/cran.r-project.org\/package=MAPA\" target=\"_blank\" rel=\"noopener noreferrer\">MAPA<\/a> and <a href=\"https:\/\/cran.r-project.org\/package=thief\" target=\"_blank\" rel=\"noopener noreferrer\">thief<\/a>. The former implements, as the name suggests, MAPA, while the latter provides code to use Temporal Hierarchies.<\/p>\n<div class=\"table-responsive\"><table  style=\"width:100%;  margin-left:auto;margin-right:auto\"  class=\"easy-table easy-table-default \" >\n<caption>Table 1. Papers on MTA<\/caption>\n<thead>\r\n<tr><th  style=\"width:220px;text-align:left\" >Paper<\/th>\n<th  style=\"width:300px;text-align:left\" >Summary<\/th>\n<\/tr>\n<\/thead>\n<tbody>\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2014\/04\/19\/improving-forecasting-by-estimating-time-series-structural-components-across-multiple-frequencies\/\">Kourentzes et al. 2014. Improving forecasting by estimating time series structural components across multiple frequencies.<\/a><\/td>\n<td  style=\"text-align:left\" > The initial paper on MTA modelling. It introduces the Multiple Aggregation Prediction Algorithm (MAPA) and demonstrates its superior performance on the well-known M3 competition.<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2014\/05\/13\/forecast-combinations-for-intermittent-demand\/\">Petropoulos and Kourentzes 2014. Forecast combinations for intermittent demand.<\/a><\/td>\n<td  style=\"text-align:left\" > Expands MAPA for the case of intermittent demand.<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2015\/09\/08\/forecasting-with-multivariate-temporal-aggregation-the-case-of-promotional-modelling\/\">Kourentzes and Petropoulos 2016. Forecasting with multivariate temporal aggregation: The case of promotional modelling.<\/a><\/td>\n<td  style=\"text-align:left\" > Expands MAPA for promotional modelling purposes at Stock Keeping Unit level.<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2016\/03\/30\/distributions-of-forecasting-errors-of-forecast-combinations-implications-for-inventory-management\/\">Barrow and Kourentzes 2016. Distributions of forecasting errors of forecast combinations: implications for inventory management.<\/a><\/td>\n<td  style=\"text-align:left\" >Provides evidence of very strong performance of MAPA over established benchmarks for demand forecasting and inventory management purposes.<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/02\/27\/forecasting-with-temporal-hierarchies-3\/\">Athanasopoulos et al. 2017. Forecasting with temporal hierarchies.<\/a><\/td>\n<td  style=\"text-align:left\" >Introduces a general framework for MTA: Temporal Hierarchies that allows use of any model\/method to produce forecasts at each level.<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" ><a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/04\/22\/demand-forecasting-by-temporal-aggregation-using-optimal-or-multiple-aggregation-levels\/\">Kourentzes et al. 2017. Demand forecasting by temporal aggregation: using optimal or multiple aggregation levels?<\/a><\/td>\n<td  style=\"text-align:left\" >Demonstrates that MTA modelling is more robust to uncertainty than modelling either using the original data or using a single (optimal) temporal aggregation level.<\/td>\n<\/tr>\n<\/tbody><\/table><\/div>\n<p style=\"text-align: justify;\">To give you an idea of the reported improvements, I have collated some of the results from the papers above. The best forecast in each column, in all tables, is highlighted in boldface. Table 2 provides a summary for the quarterly and monthly M3 datasets, using as benchmarks the Exponential Smoothing (ETS) family of models, with automatic model selection (via AICc), and Theta, the best performing method on the original M3 competition \u2013 a position it held for almost 15 years! In this case both MAPA and Temporal Hierarchies make use of the ETS family of models, so you can get a feeling of the improvement provided by MTA over conventional time series forecasting, as the results are directly comparable with the ETS row.<\/p>\n<p style=\"text-align: justify;\">Tables 3 and 4 provide results for a number of real datasets. Table 4 also provides results on a variety of simulated ARIMA series. The detailed results can be found in the respective papers. In all cases MAPA is better, or at least as good, compared to the various benchmarks. Table 5 provides results on real series that have promoted periods. There are two comparisons: forecasts without and with promotional information. In both cases MTA based forecasts (MAPA) are on average the most accurate.<\/p>\n<p><div class=\"table-responsive\"><table  style=\"width:90%;  margin-left:auto;margin-right:auto\"  class=\"easy-table easy-table-default \" >\n<caption>Table 2. sMAPE results on M3 quarterly and monthly data<sup>1<\/sup><\/caption>\n<thead>\r\n<tr><th  style=\"width:100px;text-align:left\" >Forecast<\/th>\n<th  style=\"width:50px;text-align:center\" >Quarterly set<\/th>\n<th  style=\"width:50px;text-align:center\" > Monthly set<\/th>\n<\/tr>\n<\/thead>\n<tbody>\r\n<tr><td  style=\"text-align:left\" >Exponential Smoothing (ETS)<\/td>\n<td  style=\"text-align:center\" >9.94%<\/td>\n<td  style=\"text-align:center\" >14.45%<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Theta (<a href=\"http:\/\/www.forecasters.org\/pdfs\/ijf\/Makridakia-The%20M3%20Competition.pdf\" target=\"_blank\" rel=\"noopener noreferrer\">M3 competition<\/a>)<sup>2<\/sup><\/td>\n<td  style=\"text-align:center\" ><strong>8.96%<\/strong><\/td>\n<td  style=\"text-align:center\" >13.85%<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >MAPA (<a href=\"http:\/\/kourentzes.com\/forecasting\/2014\/04\/19\/improving-forecasting-by-estimating-time-series-structural-components-across-multiple-frequencies\/\">Kourentzes et al. 2014<\/a>)<\/td>\n<td  style=\"text-align:center\" >9.58%<\/td>\n<td  style=\"text-align:center\" >13.69%<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Temporal Hierarchies (<a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/02\/27\/forecasting-with-temporal-hierarchies-3\/\">Athanasopoulos et al. 2017<\/a>)<\/td>\n<td  style=\"text-align:center\" >9.70%<\/td>\n<td  style=\"text-align:center\" ><strong>13.61%<\/strong><\/td>\n<\/tr>\n<\/tbody><\/table><\/div><br \/>\n<sup>1<\/sup> Papers provide results on more robust metrics!<br \/>\n<sup>2<\/sup> Best performance in the original M3 competition.<\/p>\n<div class=\"table-responsive\"><table  style=\"width:90%;  margin-left:auto;margin-right:auto\"  class=\"easy-table easy-table-default \" >\n<caption>Table 3. Scaled RMSE results on Fast Moving Consumer Goods sales (Barrow and Kourentzes, 2016)<\/caption>\n<thead>\r\n<tr><th  style=\"width:100px;text-align:left\" >Forecast<\/th>\n<th  style=\"width:50px;text-align:center\" >1-step ahead<\/th>\n<th  style=\"width:50px;text-align:center\" >3-steps ahead<\/th>\n<th  style=\"width:50px;text-align:center\" >5-steps ahead<\/th>\n<\/tr>\n<\/thead>\n<tbody>\r\n<tr><td  style=\"text-align:left\" >Naive<\/td>\n<td  style=\"text-align:center\" >0.882<\/td>\n<td  style=\"text-align:center\" >0.900<\/td>\n<td  style=\"text-align:center\" >0.919<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >ETS<\/td>\n<td  style=\"text-align:center\" > 0.677<\/td>\n<td  style=\"text-align:center\" >0.688<\/td>\n<td  style=\"text-align:center\" >0.711<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >AR<\/td>\n<td  style=\"text-align:center\" > 0.707<\/td>\n<td  style=\"text-align:center\" >0.719<\/td>\n<td  style=\"text-align:center\" >0.737<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >ARIMA<\/td>\n<td  style=\"text-align:center\" > 1.446<\/td>\n<td  style=\"text-align:center\" >0.701<\/td>\n<td  style=\"text-align:center\" >0.721<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Theta<\/td>\n<td  style=\"text-align:center\" > 0.674<\/td>\n<td  style=\"text-align:center\" >0.685<\/td>\n<td  style=\"text-align:center\" >0.705<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >MAPA<\/td>\n<td  style=\"text-align:center\" > <strong>0.668<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.670<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.687<\/strong><\/td>\n<\/tr>\n<\/tbody><\/table><\/div>\n<div class=\"table-responsive\"><table  style=\"width:90%;  margin-left:auto;margin-right:auto\"  class=\"easy-table easy-table-default \" >\n<caption>Table 4. Average Relative MAE on simulated and real data (Kourentzes et al., 2017)<\/caption>\n<thead>\r\n<tr><th  style=\"width:100px;text-align:left\" >Forecast<\/th>\n<th  style=\"width:50px;text-align:center\" >Simulated ARIMA<\/th>\n<th  style=\"width:50px;text-align:center\" > Manaufacturing<\/th>\n<th  style=\"width:50px;text-align:center\" > Call centre<\/th>\n<\/tr>\n<\/thead>\n<tbody>\r\n<tr><td  style=\"text-align:left\" >Single Exponential Smoothing (SES)<\/td>\n<td  style=\"text-align:center\" >1.000<\/td>\n<td  style=\"text-align:center\" >1.000<\/td>\n<td  style=\"text-align:center\" >1.000<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Exponential Smoothing (ETS)<\/td>\n<td  style=\"text-align:center\" > 0.985<\/td>\n<td  style=\"text-align:center\" >1.011<\/td>\n<td  style=\"text-align:center\" >1.005<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Optimal Temporal Aggregation &amp; SES<\/td>\n<td  style=\"text-align:center\" > 0.974<\/td>\n<td  style=\"text-align:center\" > 0.999<\/td>\n<td  style=\"text-align:center\" >1.080<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >MAPA<\/td>\n<td  style=\"text-align:center\" > <strong>0.971<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.994<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.979<\/strong><\/td>\n<\/tr>\n<\/tbody><\/table><\/div>\n<div class=\"table-responsive\"><table  style=\"width:90%;  margin-left:auto;margin-right:auto\"  class=\"easy-table easy-table-default \" >\n<caption>Table 5. Scaled MAE results on SKUs with promotions (Kourentzes and Petropoulos, 2016)<\/caption>\n<thead>\r\n<tr><th  style=\"width:100px;text-align:left\" >Forecast<\/th>\n<th  style=\"width:50px;text-align:center\" >4-steps ahead<\/th>\n<th  style=\"width:50px;text-align:center\" >8-steps ahead<\/th>\n<th  style=\"width:50px;text-align:center\" >12-steps ahead<\/th>\n<\/tr>\n<\/thead>\n<tbody>\r\n<tr><td  style=\"text-align:left\" >Naive<\/td>\n<td  style=\"text-align:center\" > 0.743<\/td>\n<td  style=\"text-align:center\" > 0.818<\/td>\n<td  style=\"text-align:center\" > 0.704<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >ETS<\/td>\n<td  style=\"text-align:center\" > 0.704<\/td>\n<td  style=\"text-align:center\" > 0.774<\/td>\n<td  style=\"text-align:center\" > 0.701<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >MAPA<\/td>\n<td  style=\"text-align:center\" > 0.679<\/td>\n<td  style=\"text-align:center\" > 0.754<\/td>\n<td  style=\"text-align:center\" > 0.736<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >Regression + Promotional<\/td>\n<td  style=\"text-align:center\" > 0.611<\/td>\n<td  style=\"text-align:center\" > 0.659<\/td>\n<td  style=\"text-align:center\" > 0.714<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >ETS + Promotional<\/td>\n<td  style=\"text-align:center\" > 0.642<\/td>\n<td  style=\"text-align:center\" > 0.627<\/td>\n<td  style=\"text-align:center\" > 0.543<\/td>\n<\/tr>\n\r\n<tr><td  style=\"text-align:left\" >MAPA + Promotional<\/td>\n<td  style=\"text-align:center\" > <strong>0.525<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.521<\/strong><\/td>\n<td  style=\"text-align:center\" ><strong>0.515<\/strong><\/td>\n<\/tr>\n<\/tbody><\/table><\/div>\n<p style=\"text-align: justify;\">The main argument in all papers is that MTA helps to improve forecast accuracy due to the way it mitigates modelling uncertainty. As we will see this comes at no additional data cost and relatively limited additional computations. An added benefit, which is not very evident from the summarised tables provided here, is that the MTA forecasts are reliable both for short and long term forecasting, providing a way to reconcile operational, tactical and strategic planning.<\/p>\n<p style=\"text-align: justify;\">Unpublished results on different applications provide a similar picture in terms of accuracy. There is also evidence that MTA can strengthen statistical tests, as the initial results of <a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/03\/30\/can-you-spot-a-trend-in-a-time-series\/\">this experiment<\/a> show. However, all this is ongoing research, so until a full analysis is conducted and the results are peer reviewed, I would add a pinch of salt to these!<\/p>\n<p style=\"text-align: justify;\">In following blog posts I will explain how MTA works and elaborate more on results from the various papers.<\/p>\n<p style=\"text-align: justify;\">Multiple Temporal Aggregation: the story so far: Part I; <a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/05\/04\/multiple-temporal-aggregation-the-story-so-far-part-ii\/\">Part II<\/a>; <a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/07\/04\/multiple-temporal-aggregation-the-story-so-far-part-iii\/\">Part III<\/a>; <a href=\"http:\/\/kourentzes.com\/forecasting\/2017\/07\/22\/multiple-temporal-aggregation-the-story-so-far-part-iv-temporal-hierarchies\/\">Part IV<\/a>.<\/p>\n<div class=\"SPOSTARBUST-Related-Posts\"><H3>Related Posts<\/H3><ul class=\"entry-meta\"><li class=\"SPOSTARBUST-Related-Post\"><a title=\"Stochastic Coherency in Forecast Reconciliation\" href=\"https:\/\/kourentzes.com\/forecasting\/2021\/07\/09\/stochastic-coherency-in-forecast-reconciliation\/\" rel=\"bookmark\">Stochastic Coherency in Forecast Reconciliation<\/a><\/li>\n<li class=\"SPOSTARBUST-Related-Post\"><a title=\"Visitor Arrivals Forecasts amid COVID-19: A Perspective from the Africa Team\" href=\"https:\/\/kourentzes.com\/forecasting\/2021\/07\/09\/visitor-arrivals-forecasts-amid-covid-19-a-perspective-from-the-africa-team\/\" rel=\"bookmark\">Visitor Arrivals Forecasts amid COVID-19: A Perspective from the Africa Team<\/a><\/li>\n<li class=\"SPOSTARBUST-Related-Post\"><a title=\"ISF 2020: A geometry inspired hierarchical forecasting methodology\" href=\"https:\/\/kourentzes.com\/forecasting\/2020\/10\/25\/isf-2020-using-information-from-the-business-environment-to-improve-forecasting\/\" rel=\"bookmark\">ISF 2020: A geometry inspired hierarchical forecasting methodology<\/a><\/li>\n<\/ul><\/div><!-- AddThis Advanced Settings generic via filter on the_content --><!-- AddThis Share Buttons generic via filter on the_content -->","protected":false},"excerpt":{"rendered":"<p>Over the last years I have been working (with my co-authors!) on the idea of Multiple Temporal Aggregation (MTA) for time series forecasting. A number of papers have been published introducing and developing the idea further, or testing its effectiveness for forecasting. In this series of blog posts I will try to summarise the progress\u2026 <span class=\"read-more\"><a href=\"https:\/\/kourentzes.com\/forecasting\/2017\/04\/27\/multiple-temporal-aggregation-the-story-so-far-part-i\/\">Read More &raquo;<\/a><\/span><!-- AddThis Advanced Settings generic via filter on get_the_excerpt --><!-- AddThis Share Buttons generic via filter on get_the_excerpt --><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":[],"categories":[41],"tags":[32,76,36,68],"_links":{"self":[{"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/posts\/1305"}],"collection":[{"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/comments?post=1305"}],"version-history":[{"count":22,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/posts\/1305\/revisions"}],"predecessor-version":[{"id":1397,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/posts\/1305\/revisions\/1397"}],"wp:attachment":[{"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/media?parent=1305"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/categories?post=1305"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/kourentzes.com\/forecasting\/wp-json\/wp\/v2\/tags?post=1305"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}<!-- WP Super Cache is installed but broken. The constant WPCACHEHOME must be set in the file wp-config.php and point at the WP Super Cache plugin directory. -->