Robust calculation of slopes in detrended fluctuation analysis and its application to envelopes of human alpha rhythms

Detrended fluctuation analysis (DFA) is a popular method to analyze long-range temporal correlations in time series of many different research areas but in particular also for electrophysiological recordings. Using the classical DFA method, the cumulative sum of data are divided into segments, and the variance of these sums is studied as a function of segment length after linearly detrending them in each segment. The starting point of the proposed new method is the observation that the classical method is inherently non-stationary without justification by a corresponding non-stationarity of the data. This leads to unstable estimates of fluctuations to the extent that it is impossible to estimate slopes of the fluctuations other than by fitting a line over a wide range of temporal scales. We here use a modification of the classical method by formulating the detrending as a strictly stationary operation. With this modification the detrended fluctuations can be expressed as a weighted average across the power spectrum of a signal. Most importantly, we can also express the slopes, calculated as analytic derivatives of the fluctuations with respect to the scales, as statistically robust weighted averages across the power spectra. The method is applied to amplitudes of brain oscillations measured with magnetoencephalography in resting state condition. We found for envelopes of the the alpha rhythm that fluctuations as a function of time scales in a double-logarithmic plot differ substantially from a linear relation for time scales below 10 seconds. In particular we will show that model selections fail to determine accurate scaling laws, and that standard parameter settings are likely to yield results depending on signal to noise ratios than on true long range temporal correlations.