Optimal Sparse Estimation of High Dimensional Heavy-tailed Time Series

TL;DR

This paper develops optimal estimation techniques for high-dimensional heavy-tailed VAR models, achieving rates comparable to iid data and allowing for more general sparsity-inducing penalties without relying on mixing conditions.

Contribution

It establishes optimal consistency rates for heavy-tailed VARs with temporal dependence and introduces a new concentration bound for heavy-tailed processes, enabling more flexible penalties.

Findings

Optimal finite-sample bounds matching iid rates

New concentration bounds for heavy-tailed processes

Enhanced estimation methods for complex sparsity patterns

Abstract

Recently, high dimensional vector auto-regressive models (VAR), have attracted a lot of interest, due to novel applications in the health, engineering and social sciences. The presence of temporal dependence poses additional challenges to the theory of penalized estimation techniques widely used in the analysis of their iid counterparts. However, recent work (e.g., [Basu and Michailidis, 2015, Kock and Callot, 2015]) has established optimal consistency of $\ell_1$-LASSO regularized estimates applied to models involving high dimensional stable, Gaussian processes. The only price paid for temporal dependence is an extra multiplicative factor that equals 1 for independent and identically distributed (iid) data. Further, [Wong et al., 2020] extended these results to heavy tailed VARs that exhibit "$\beta$-mixing" dependence, but the rates rates are sub-optimal, while the extra factor is intractable. This paper improves these results in two important directions: (i) We establish optimal consistency rates and corresponding finite sample bounds for the underlying model parameters that match those for iid data, modulo a price for temporal dependence, that is easy to interpret and equals 1 for iid data. (ii) We incorporate more general penalties in estimation (which are not decomposable unlike the $\ell_1$ norm) to induce general sparsity patterns. The key technical tool employed is a novel, easy-to-use concentration bound for heavy tailed linear processes, that do not rely on "mixing" notions and give tighter bounds.