Autoregressive with Slack Time Series Model for Forecasting a Partially-Observed Dynamical Time Series
Akifumi Okuno, Yuya Morishita, Yoh-ichi Mototake
TL;DR
The paper tackles forecasting partially observed dynamical time series defined by an evolution function, where some state variables are unobserved. It proposes the autoregressive with slack time Series (ARS) model, which jointly estimates the evolution matrix $\boldsymbol B$ and imputes missing variables as slack time series $\{\boldsymbol z^{\dagger}_j\}$, leveraging a time-invariant linear dynamical assumption. A theoretical guarantee in the 2D case shows ARS can reconstruct the underlying true dynamics from a single observed dimension, and numerical experiments on circular motion and Lorenz-like dynamics demonstrate improved forecasting over conventional AR. The work clarifies connections to state-space formulations and higher-order AR models, discusses limitations such as overparameterization and optimization efficiency, and outlines future directions including extensions to interaction terms and practical applications in plasma physics.
Abstract
This study delves into the domain of dynamical systems, specifically the forecasting of dynamical time series defined through an evolution function. Traditional approaches in this area predict the future behavior of dynamical systems by inferring the evolution function. However, these methods may confront obstacles due to the presence of missing variables, which are usually attributed to challenges in measurement and a partial understanding of the system of interest. To overcome this obstacle, we introduce the autoregressive with slack time series (ARS) model, that simultaneously estimates the evolution function and imputes missing variables as a slack time series. Assuming time-invariance and linearity in the (underlying) entire dynamical time series, our experiments demonstrate the ARS model's capability to forecast future time series. From a theoretical perspective, we prove that a 2-dimensional time-invariant and linear system can be reconstructed by utilizing observations from a single, partially observed dimension of the system.
