Hierarchical Time Series Forecasting with Robust Reconciliation

This paper focuses on forecasting hierarchical time-series data, where each higher-level observation equals the sum of its corresponding lower-level time series. In such contexts, the forecast values should be coherent, meaning that the forecast value of each parent series exactly matches the sum of the forecast values of its child series. Existing hierarchical forecasting methods typically generate base forecasts independently for each series and then apply a reconciliation procedure to adjust them so that the resulting forecast values are coherent across the hierarchy. These methods generally derive an optimal reconciliation, using a covariance matrix of the forecast error. In practice, however, the true covariance matrix is unknown and has to be estimated from finite samples in advance. This gap between the true and estimated covariance matrix may degrade forecast performance. To address this issue, we propose a robust optimization framework for hierarchical reconciliation that accounts for uncertainty in the estimated covariance matrix. We first introduce an uncertainty set for the estimated covariance matrix and formulate a reconciliation problem that minimizes the worst-case average of weighted squared residuals over this uncertainty set. We show that our problem can be cast as a semidefinite optimization problem. Numerical experiments demonstrate that the proposed robust reconciliation method achieved better forecast performance than existing hierarchical forecasting methods, which indicates the effectiveness of integrating uncertainty into the reconciliation process.

💡 Research Summary

The paper addresses the problem of producing coherent forecasts for hierarchical time‑series, where each parent series must equal the sum of its children. Traditional reconciliation methods such as Bottom‑Up, Top‑Down, GLS, and the Minimum‑Trace (MinT) approach rely on an estimated error covariance matrix (W). In practice this matrix is unknown and must be inferred from limited residual data, which introduces estimation error that can degrade forecast accuracy.

To mitigate this vulnerability, the authors formulate a robust optimization framework that treats the inverse covariance matrix (W^{-1}) as an uncertain parameter. They define a box‑type uncertainty set (\mathcal M={M\mid M_{\text{low}}\preceq M\preceq M_{\text{up}},;M\succeq0}), where the lower and upper bound matrices encode element‑wise confidence intervals around the estimated covariance. The robust reconciliation problem then minimizes the worst‑case average weighted squared residual over this set:

\