Abstract:Classical nonparametric spectral analysis uses sliding windows to capture the dynamic nature of most real-world time series. This universally accepted approach fails to exploit the temporal continuity in the data and is not well-suited for signals with highly structured time-frequency representations. For a time series whose time-varying mean is the superposition of a small number of oscillatory components, we formulate nonparametric batch spectral analysis as a Bayesian estimation problem. We introduce prior distributions on the time-frequency plane that yield maximum a posteriori (MAP) spectral estimates that are continuous in time yet sparse in frequency. Our spectral decomposition procedure, termed spectrotemporal pursuit, can be efficiently computed using an iteratively reweighted least-squares algorithm and scales well with typical data lengths. We show that spectrotemporal pursuit works by applying to the time series a set of data-derived filters. Using a link between Gaussian mixture models, l1 minimization, and the expectation-maximization algorithm, we prove that spectrotemporal pursuit converges to the global MAP estimate. We illustrate our technique on simulated and real human EEG data as well as on human neural spiking activity recorded during loss of consciousness induced by the anesthetic propofol. For the EEG data, our technique yields significantly denoised spectral estimates that have significantly higher time and frequency resolution than multitaper spectral estimates. For the neural spiking data, we obtain a new spectral representation of neuronal firing rates. Spectrotemporal pursuit offers a robust spectral decomposition framework that is a principled alternative to existing methods for decomposing time series into a small number of smooth oscillatory components.

A new approach to nonparametric estimation of multivariate spectral density function using basis expansion

A nonparametrically corrected likelihood for Bayesian spectral analysis of multivariate time series

Multidimensional Variational Line Spectra Estimation.

Adaptive Bayesian Spectral Analysis of High-Dimensional Nonstationary Time Series

Bayesian Spectral Modeling for Multiple Time Series

A New Approach to Spectral Estimation: a Tunable High-Resolution Spectral Estimator

Nonparametric collective spectral density estimation with an application to clustering the brain signals

Multivariate Nonparametric Volatility Density Estimation

Spectral Initialization for High-Dimensional Phase Retrieval with Biased Spatial Directions

A Nonlinear Subspace Approach for Parametric Estimation of PDFs from Short Data Records with Application to Rayleigh Fading

Spectral Analysis of High-dimensional Time Series

Spectral Density Estimation for a Class of Spectrally Correlated Processes

On Generalized Expectation-Based Estimation of a Population Spectral Distribution from High-Dimensional Data

Nonparametric Estimation of Band-limited Probability Density Functions

Multivariate, Heteroscedastic Empirical Bayes via Nonparametric Maximum Likelihood

Two new algorithms for maximum likelihood estimation of sparse covariance matrices with applications to graphical modeling

Nonparametric methods for volatility density estimation

Posterior consistency for the spectral density of non-Gaussian stationary time series

Approximate Message Passing for orthogonally invariant ensembles: Multivariate non-linearities and spectral initialization

Robust spectrotemporal decomposition by iteratively reweighted least squares

Sublinear Time Spectral Density Estimation