An Orchestra Where Every Instrument Enters on Its Own Schedule ✎

Original: "Adaptive Bayesian Power Spectrum Analysis of Multivariate Nonstationary Time Series"
· Zeda Li, Robert T. Krafty
arXiv:1706.05661 · 2017-06-18 · arXiv non-exclusive · independent analysis · 2 min
A new method splits a complex signal into segments, allowing different parts to change independently, and pinpoints moments of change with high accuracy.

An orchestra plays a complex piece, and each group of instruments enters on its own schedule. The violins might slow down while the woodwinds keep the same tempo, and the percussion might suddenly burst into a roll. The conductor hears everything at once, but pinpointing the exact moment when the flutes, not the cellos, changed rhythm is nearly impossible. That is what signals from the brain during sleep or from the ocean surface during El Niño years look like: many voices, each living its own life.

Standard analysis of such a signal either assumes all instruments change simultaneously or requires changes to be smooth. The new method allows each part of the signal to have its own schedule of changes.

Previously, researchers could either fix common breakpoints for all channels or smooth over changes, losing sharp jumps. The new approach works like an attentive listener who divides the entire recording into short segments where the music sounds roughly the same and builds a frequency picture for each segment. The number of segments is not known in advance, and the method decides where to draw the boundaries by exploring options through random jumps and smart search.

On a sleep brain-wave recording, the method found the transition from slow-wave to REM sleep with 94 percent probability, and the boundary fell within a quarter of a second of the one marked manually by a physician.

For El Niño data, the method suggested that the character of oscillations changed in the 1980s, assigning this event a probability of 65 percent. This is not a final verdict but rather a hint for climatologists about where to look for causes of change. In tests on synthetic signals with known breakpoints, the method made fewer errors than its competitors and ran faster than one of the most accurate previous algorithms.

Such a tool will be useful wherever the dynamics of connections between multiple signals matter: in sleep medicine for diagnosing disorders, in neuroscience for analyzing activity across different brain regions, and in climatology for monitoring ocean and atmosphere indices. The orchestra is complex, but now we have a way to hear when each group of instruments changes its tempo.

In EEG analysis, the method automatically identified the transition from slow-wave to REM sleep within a quarter of a second, comparable to the performance of an experienced clinician.

The idea of adaptively partitioning a signal into segments with different dynamics echoes the analysis of brain waves in the movie Inception, where characters track changes in states of consciousness through activity patterns.

f^{-1}(u,\omega) = \Theta(u,\omega) \Psi(u,\omega)^{-1} \Theta(u,\omega)^*
The inverse spectral matrix is expressed through a lower triangular matrix Θ and a positive diagonal matrix Ψ, ensuring positive definiteness.
Concepts
nonlinear laplacian spectral analysis wavelet analysis Variational Method bayesian hierarchical model probabilistic model wavelet transform noninvasive brain monitoring gaussian process metropolis hastings algorithm parameter posterior distribution power spectrum time reversal permutation test nested sampling Boltzmann equation Bayesian inference stochastic processes spectral analysis ai empirical orthogonal function