ξ-π: a nonparametric model for neural power spectra decomposition

Abstract

AbstractThe power spectra estimated from the brain recordings are the mixed representation of aperiodic transient activity and periodic oscillations, i.e., aperiodic component (AC) and periodic component (PC). Quantitative neurophysiology requires precise decomposition preceding parameterizing each component. However, the shape, statistical distribution, scale, and mixing mechanism of AC and PCs are unclear, challenging the effectiveness of current popular parametric models such as FOOOF, IRASA, BOSC, etc. Here,ξ-πwas proposed to decompose the neural spectra by embedding the nonparametric spectra estimation with penalized Whittle likelihood and the shape language modeling into the expectation maximization frame-work.ξ-πwas validated on the synthesized spectra with loss statistics and on the sleep EEG and the large sample iEEG with evaluation metrics and neurophysiological evidence. Compared to FOOOF, both the simulation presenting shape irregularities and the batch simulation with multiple isolated peaks indicated thatξ-πimproved the fit of AC and PCs with less loss and higher F1-score in recognizing the centering frequencies and the number of peaks; the sleep EEG revealed thatξ-πproduced more distinguishable AC exponents and improved the sleep state classification accuracy; the iEEG showed thatξ-πapproached the clinical findings in peak discovery. Overall,ξ-πoffered good performance in the spectra decomposition, which allows flexible parameterization using descriptive statistics or kernel functions.ξ-πmay be a promising tool for brain signal decoding in fields such as cognitive neuroscience, brain-computer interface, neurofeedback, and brain diseases.