Robust Estimation of 1/f Activity Improves Oscillatory Burst Detection

Abstract

AbstractNeural oscillations often occur as transient bursts with variable amplitude and frequency dynamics. Quantifying these effects is important for understanding brain-behaviour relationships, especially in continuous datasets. To robustly measure bursts, rhythmical periods of oscillatory activity must be separated from arrhythmical background 1/f activity, which is ubiquitous in electrophysiological recordings. The Better OSCillation (BOSC) framework achieves this by defining a power threshold above the estimated background 1/f activity, combined with a duration threshold. Here we introduce a modification to this approach called fBOSC which uses a spectral parametrisation tool to accurately model background 1/f activity in neural data. fBOSC (which is openly available as a MATLAB toolbox) is robust to power spectra with oscillatory peaks and can also model non-linear spectra. Through a series of simulations, we show that fBOSC more accurately models the 1/f power spectrum compared with existing methods. fBOSC was especially beneficial where power spectra contained a “knee” below ∼0.5-10 Hz, which is typical in neural data. We also found that, unlike other methods, fBOSC was unaffected by oscillatory peaks in the neural power spectrum. Moreover, by robustly modelling background 1/f activity, the sensitivity for detecting oscillatory bursts was standardised across frequencies (e.g. theta- and alpha-bands). Finally, using openly available resting state magnetoencephalography and intracranial electrophysiology datasets, we demonstrate the application of fBOSC for oscillatory burst detection in the theta-band. These simulations and empirical analyses highlight the value of fBOSC in detecting oscillatory bursts, including in datasets that are long and continuous with no distinct experimental trials.GRAPHICAL ABSRACTTo determine a power threshold for burst detection, the Better OSCillation framework (BOSC) estimates background 1/f activity by modelling neural power spectra. Here we introduce a modification, termed fBOSC, to more robustly estimate 1/f activity in situations with prominent oscillatory peaks and/or the presence of a non-linear “knee” in the power spectrum. This was shown to standardise burst detection across frequency bands in both simulated and empirical data.