Fast Machine-Precision Spectral Likelihoods for Stationary Time Series

We provide in this work an algorithm for approximating a very broad class of symmetric Toeplitz matrices to machine precision in time with applications to fitting time series models. In particular, for a symmetric Toeplitz matrix with values where is piecewise smooth, we give an approximation , where is the DFT matrix, is diagonal, and the matrices and are in with . Studying these matrices in the context of time series, we offer a theoretical explanation of this structure and connect it to existing spectral-domain approximation frameworks. We then give a complete discussion of the numerical method for assembling the approximation and demonstrate its efficiency for improving Whittle-type likelihood approximations, including dramatic examples where a correction of rank to the standard Whittle approximation increases the accuracy from to digits for a matrix . The method and analysis of this work applies well beyond time series analysis, providing an algorithm for extremely accurate direct solves with a wide variety of symmetric Toeplitz matrices. The analysis employed here largely depends on asymptotic expansions of oscillatory integrals, and also provides a new perspective on when existing spectral-domain approximation methods for Gaussian log-likelihoods can be particularly problematic.
View on arXiv