Frontiers in Neurology (Dec 2013)
Non-Linear Dynamical Analysis of EEG Time Series Distinguishes Patients with Parkinson’s Disease from Healthy Individuals
Abstract
The pathophysiology of Parkinson’s disease (PD) is known to involve altered patterns of neuronal firing and synchronization in cortical-basal ganglia circuits. We used delay differential equations (DDE) as non-linear time-domain classification tools to analyze electroencephalographic (EEG) recordings from PD patients on and off dopaminergic therapy and healthy individuals. Two sets of 50 1-s segments of 64-channel EEG activity were recorded from nine PD patients on and off medication and nine age-matched controls. The 64 EEG channels were grouped into 10 clusters covering different brain regions for analysis. DDE models were fitted to individual trials, and model coefficients and error were used as features for classification. The best models were selected using repeated random sub-sampling and classification performance was measured using the area under the ROC curve A'. In a companion paper, we show that DDEs can uncover hidden dynamical structure from short segments of simulated time series of known dynamical systems in high noise regimes. Using the same method for finding the best models, we found here that even short segments of resting state EEG in PD patients and controls contained dynamical structure, and, moreover, that PD patients exhibited a greater dynamic range than controls. DDE model output on the means from one set of 50 trials provided nearly complete separation of PD patients off medication from controls across brain regions A' varied from 0.95-1.0. For distinguishing PD patients on vs. off medication, classification performance A' ranged from 0.86-1.0 across brain regions. Moreover, the generalizability of the model to the second set of 50 trials was excellent, with A' ranging from 0.81-0.94 across brain regions for controls vs. PD off medication, and from 0.62-0.82 for PD on medication vs. off. Finally, model features significantly predicted individual patients' motor severity, as assessed with standard clinical rating scales.
Keywords