A maximum likelihood based technique for validating detrended fluctuation analysis (ML-DFA)

Detrended Fluctuation Analysis (DFA) is widely used to assess the presence of long-range temporal correlations in time series. Signals with long-range temporal correlations are typically defined as having a power law decay in their autocorrelation function. The output of DFA is an exponent, which is the slope obtained by linear regression of a log-log fluctuation plot against window size. However, if this fluctuation plot is not linear, then the underlying signal is not self-similar, and the exponent has no meaning. There is currently no method for assessing the linearity of a DFA fluctuation plot. Here we present such a technique, called ML-DFA. We scale the DFA fluctuation plot to construct a likelihood function for a set of alternative models including polynomial, root, exponential, logarithmic and spline functions. We use this likelihood function to determine the maximum likelihood and thus to calculate values of the Akaike and Bayesian information criteria, which identify the best fit model when the number of parameters involved is taken into account and over-fitting is penalised. This ensures that, of the models that fit well, the least complicated is selected as the best fit. We apply ML-DFA to synthetic data from FARIMA processes and sine curves with DFA fluctuation plots whose form has been analytically determined, and to experimentally collected neurophysiological data. ML-DFA assesses whether the hypothesis of a linear fluctuation plot should be rejected, and thus whether the exponent can be considered meaningful. We argue that ML-DFA is essential to obtaining trustworthy results from DFA.Comment: 22 pages, 7 figure

Modelling and analysis of amplitude, phase and synchrony in human brain activity patterns

The critical brain hypothesis provides a framework for viewing the human brain as a critical system, which may transmit information, reorganise itself and react to external stimuli efficiently. A critical system incorporates structures at a range of spatial and temporal scales, and may be associated with power law distributions of neuronal avalanches and power law scaling functions. In the temporal domain, the critical brain hypothesis is supported by a power law decay of the autocorrelation function of neurophysiological signals, which indicates the presence of long-range temporal correlations (LRTCs). LRTCs have been found to exist in the amplitude envelope of neurophysiological signals such as EEG, EMG and MEG, which reveal patterns of local synchronisation within neuronal pools. Synchronisation is an important tool for communication in the nervous system and can also exist between disparate regions of the nervous system. In this thesis, inter-regional synchronisation is characterised by the rate of change of phase difference between neurophysiological time series at different neuronal regions and investigated using the novel phase synchrony analysis method. The phase synchrony analysis method is shown to recover the DFA exponents in time series where these are known. The method indicates that LRTCs are present in the rate of change of phase difference between time series derived from classical models of criticality at critical parameters, and in particular the Ising model of ferromagnetism and the Kuramoto model of coupled oscillators. The method is also applied to the Cabral model, in which Kuramoto oscillators with natural frequencies close to those of cortical rhythms are embedded in a network based on brain connectivity. It is shown that LRTCs in the rate of change of phase difference are disrupted when the network properties of the system are reorganised. The presence of LRTCs is assessed using detrended fluctuation analysis (DFA), which assumes the linearity of a log-log plot of detrended fluctuation magnitude. In this thesis it is demonstrated that this assumption does not always hold, and a novel heuristic technique, ML-DFA, is introduced for validating DFA results. Finally, the phase synchrony analysis method is applied to EEG, EMG and MEG time series. The presence of LRTCs in the rate of change of phase difference between time series recorded from the left and right motor cortices are shown to exist during resting state, but to be disrupted by a finger tapping task. The findings of this thesis are interpreted in the light of the critical brain hypothesis, and shown to provide motivation for future research in this area

A power-law distribution of phase-locking intervals does not imply critical interaction

Neural synchronisation plays a critical role in information processing, storage and transmission. Characterising the pattern of synchronisation is therefore of great interest. It has recently been suggested that the brain displays broadband criticality based on two measures of synchronisation - phase locking intervals and global lability of synchronisation - showing power law statistics at the critical threshold in a classical model of synchronisation. In this paper, we provide evidence that, within the limits of the model selection approach used to ascertain the presence of power law statistics, the pooling of pairwise phase-locking intervals from a non-critically interacting system can produce a distribution that is similarly assessed as being power law. In contrast, the global lability of synchronisation measure is shown to better discriminate critical from non critical interaction.Comment: (v3) Fixed error in Figure 1; (v2) Added references. Minor edits throughout. Clarified relationship between theoretical critical coupling for infinite size system and 'effective' critical coupling system for finite size system. Improved presentation and discussion of results; results unchanged. Revised Figure 1 to include error bars on r and N; results unchanged; (v1) 11 pages, 7 figure

Markers of criticality in phase synchronization

The concept of the brain as a critical dynamical system is very attractive because systems close to criticality are thought to maximize their dynamic range of information processing and communication. To date, there have been two key experimental observations in support of this hypothesis: (i) neuronal avalanches with power law distribution of size and (ii) long-range temporal correlations (LRTCs) in the amplitude of neural oscillations. The case for how these maximize dynamic range of information processing and communication is still being made and because a significant substrate for information coding and transmission is neural synchrony it is of interest to link synchronization measures with those of criticality. We propose a framework for characterizing criticality in synchronization based on an analysis of the moment-to-moment fluctuations of phase synchrony in terms of the presence of LRTCs. This framework relies on an estimation of the rate of change of phase difference and a set of methods we have developed to detect LRTCs. We test this framework against two classical models of criticality (Ising and Kuramoto) and recently described variants of these models aimed to more closely represent human brain dynamics. From these simulations we determine the parameters at which these systems show evidence of LRTCs in phase synchronization. We demonstrate proof of principle by analysing pairs of human simultaneous EEG and EMG time series, suggesting that LRTCs of corticomuscular phase synchronization can be detected in the resting state and experimentally manipulated. The existence of LRTCs in fluctuations of phase synchronization suggests that these fluctuations are governed by non-local behavior, with all scales contributing to system behavior. This has important implications regarding the conditions under which one should expect to see LRTCs in phase synchronization. Specifically, brain resting states may exhibit LRTCs reflecting a state of readiness facilitating rapid task-dependent shifts toward and away from synchronous states that abolish LRTCs

The Parkinsonian subthalamic network: measures of power, linear, and non-linear synchronization and their relationship to L-DOPA treatment and OFF state motor severity

In this paper we investigated the dopaminergic modulation of neuronal interactions occurring in the subthalamic nucleus (STN) during Parkinson's disease (PD). We utilized linear measures of local and long range synchrony such as power and coherence, as well as Detrended Fluctuation Analysis for Phase Synchrony (DFA-PS)- a recently developed non-linear method that computes the extent of long tailed autocorrelations present in the phase interactions between two coupled signals. Through analysis of local field potentials (LFPs) taken from the STN we seek to determine changes in the neurodynamics that may underpin the pathophysiology of PD in a group of 12 patients who had undergone surgery for deep brain stimulation. We demonstrate up modulation of alpha-theta (5–12 Hz) band power in response to L-DOPA treatment, whilst low beta band power (15–20 Hz) band-power is suppressed. We also find evidence for significant local connectivity within the region surrounding STN although there was evidence for its modulation via administration of L-DOPA. Further to this we present evidence for a positive correlation between the phase ordering of bilateral STN interactions and the severity of bradykinetic and rigidity symptoms in PD. Although, the ability of non-linear measures to predict clinical state did not exceed standard measures such as beta power, these measures may help identify the connections which play a role in pathological dynamics

Modeling the ongoing dynamics of short and long-range temporal correlations in broadband EEG during movement

Electroencephalogram (EEG) undergoes complex temporal and spectral changes during voluntary movement intention. Characterization of such changes has focused mostly on narrowband spectral processes such as Event-Related Desynchronization (ERD) in the sensorimotor rhythms because EEG is mostly considered as emerging from oscillations of the neuronal populations. However, the changes in the temporal dynamics, especially in the broadband arrhythmic EEG have not been investigated for movement intention detection. The Long-Range Temporal Correlations (LRTC) are ubiquitously present in several neuronal processes, typically requiring longer timescales to detect. In this paper, we study the ongoing changes in the dynamics of long- as well as short-range temporal dependencies in the single trial broadband EEG during movement intention. We obtained LRTC in 2 s windows of broadband EEG and modeled it using the Autoregressive Fractionally Integrated Moving Average (ARFIMA) model which allowed simultaneous modeling of short- and long-range temporal correlations. There were significant (p < 0.05) changes in both broadband long- and short-range temporal correlations during movement intention and execution. We discovered that the broadband LRTC and narrowband ERD are complementary processes providing distinct information about movement because eliminating LRTC from the signal did not affect the ERD and conversely, eliminating ERD from the signal did not affect LRTC. Exploring the possibility of applications in Brain Computer Interfaces (BCI), we used hybrid features with combinations of LRTC, ARFIMA, and ERD to detect movement intention. A significantly higher (p < 0.05) classification accuracy of 88.3 ± 4.2% was obtained using the combination of ARFIMA and ERD features together, which also predicted the earliest movement at 1 s before its onset. The ongoing changes in the long- and short-range temporal correlations in broadband EEG contribute to effectively capturing the motor command generation and can be used to detect movement successfully. These temporal dependencies provide different and additional information about the movement

Complexity of multi-dimensional spontaneous EEG decreases during propofol induced general anaesthesia

Emerging neural theories of consciousness suggest a correlation between a specific type of neural dynamical complexity and the level of consciousness: When awake and aware, causal interactions between brain regions are both integrated (all regions are to a certain extent connected) and differentiated (there is inhomogeneity and variety in the interactions). In support of this, recent work by Casali et al (2013) has shown that Lempel-Ziv complexity correlates strongly with conscious level, when computed on the EEG response to transcranial magnetic stimulation. Here we investigated complexity of spontaneous high-density EEG data during propofol-induced general anaesthesia. We consider three distinct measures: (i) Lempel-Ziv complexity, which is derived from how compressible the data are; (ii) amplitude coalition entropy, which measures the variability in the constitution of the set of active channels; and (iii) the novel synchrony coalition entropy (SCE), which measures the variability in the constitution of the set of synchronous channels. After some simulations on Kuramoto oscillator models which demonstrate that these measures capture distinct ‘flavours’ of complexity, we show that there is a robustly measurable decrease in the complexity of spontaneous EEG during general anaesthesia