A Novel Segmentation, Mutual Information Network Framework for EEG Analysis of Motor Tasks
© Wang et al; licensee BioMed Central Ltd. 2009
Received: 29 October 2008
Accepted: 04 May 2009
Published: 04 May 2009
Monitoring the functional connectivity between brain regions is becoming increasingly important in elucidating brain functionality in normal and disease states. Current methods of detecting networks in the recorded electroencephalogram (EEG) such as correlation and coherence are limited by the fact that they assume stationarity of the relationship between channels, and rely on linear dependencies. In contrast to diseases of the brain cortex (e.g. Alzheimer's disease), with motor disorders such as Parkinson's disease (PD) the EEG abnormalities are most apparent during performance of dynamic motor tasks, but this makes the stationarity assumption untenable.
We therefore propose a novel EEG segmentation method based on the temporal dynamics of the cross-spectrogram of the computed Independent Components (ICs). We then utilize mutual information (MI) as the metric for determining also nonlinear statistical dependencies between EEG channels. Graphical theoretical analysis is then applied to the derived MI networks. The method was applied to EEG data recorded from six normal subjects and seven PD subjects off medication. One-way analysis of variance (ANOVA) tests demonstrated statistically significant difference in the connectivity patterns between groups.
The results suggested that PD subjects are unable to independently recruit different areas of the brain while performing simultaneous tasks compared to individual tasks, but instead they attempt to recruit disparate clusters of synchronous activity to maintain behavioral performance.
The proposed segmentation/MI network method appears to be a promising approach for analyzing the EEG recorded during dynamic behaviors.
Connectivity between brain regions is important for normal brain functioning, and may be impaired in many neurological diseases . The electroencephalogram (EEG), with its excellent temporal resolution (~1 msec), is the most widely available technology used for inferring transient synchronization between brain regions. Both linear and nonlinear measures have been applied to assess the interdependencies between EEG channels . For example, coherence and correlation methods [3, 4], which measure the dependencies between a pair of EEG signals in the frequency and time domains respectively, have been applied to the EEG to study the cortical synchrony that can be modulated as a function of task, and may systematically differ between normal and disease groups [5, 6]. Nevertheless, these measures consider only linear dependencies and may be particularly sensitive to outliers. Other methods may also be used to investigate both linear and non-linear relationships between multivariate time series in the EEG, such as the Synchronization Likelihood (SL), but this and related methods assume a fixed phase relationship between time series . However, in some diseases such as PD, transient phase-locked behavior between different parts of the motor system may be interrupted by "phase slips"  making the assumption of prolonged periods of phase synchrony potentially unsuitable.
An alternative to the linear methods of coherence and correlation and phase synchronization is to consider the mutual information (MI) between channels within a specified window. This enables estimation of both the linear and nonlinear statistical dependencies between time series and can be used to detect functional coupling. MI is a statistical technique that quantifies the information transmitted from one time series to another, with maximum value when two time series are the same and a value of zero if two time series are statistically independent. Previously, researchers have utilized MI as a suitable metric to investigate EEG coupling in various pathological conditions [9–11]. For example, by estimating the MI between the time series of multiple pairs of EEG channels, Jeong et al. demonstrated abnormalities in the information transmission between different cortical areas in Alzheimer's patients . Similar studies have used MI as a marker for cortico-cortical connections in schizophrenic patients  and odor stimulation .
Another disease where altered connectivity may be important is Parkinson's disease (PD), a movement disorder that is characterized by muscle rigidity, tremor, and bradykinesia (slowing of physical movement). These symptoms do not reflect a primary failure of the cortex (making resting EEG less likely to be abnormal), but rather the effects of failure of the basal ganglia to prime the cortex for preparation and execution of movement. As a result, PD patients have a difficult time performing simultaneous movements compared to normal subjects [12, 13]. In order to assess the indirect effects of basal ganglia dysfunction on the cortex in PD, it is necessary to have subjects perform a motor task. Furthermore, stressing the motor system by having PD subjects performing simultaneous movements is more likely to induce abnormalities in the recorded EEG.
However, as soon as a subject performs a dynamic motor task, the non-stationarity nature of the EEG must be considered . The non-stationarity likely reflects the switching of inherently metastable states of neural assemblies during task performance causing abrupt transitions. The non-stationary property of EEG suggests that techniques assuming stationarity may result in misleading interpretations. To address this concern, a conventional approach is to incorporate a sliding time window into the original signal models, and assume that the stationarity assumption is valid for the segment of data in the window. However, the selection of an appropriate (possibly time-varying) window length is non-trivial and could have a significant effect on the analysis results.
In order to obtain quasi-stationary segments in EEG signals and select task-related segments, we first propose a novel segmentation method of the EEG based on the temporal dynamics of the cross-spectrogram of the Independent Components (ICs), and then compute the MI between channels within the temporally-segmented regions. We then apply graph theoretical analysis to the network of each group defined by edges whose MI values exceed a suitable threshold, and compute the clustering coefficient (C) and shortest path length (L) [15, 16]. In order to accommodate the magnitude of MI values above a given threshold, the intra-group (ie. task) and inter-group (ie. subject groups) network differences are further analyzed by one-way analysis of variance (ANOVA).
For motor tasks, changes in the EEG are most likely related to event-related synchronization/desynchronization (ERS/ERD), particularly in the beta band . In self-paced movements, ERD corresponds to changes in coherence between brain regions . Thus we suggest the use of transient synchronization of between ICs as suitable markers for segmentation [19, 20]. We emphasize that the concept of "task-related sections" is flexible and may be dependent upon the underlying behavioral paradigm subjects are asked to perform. For example, if a subject was asked to push a button every 10 seconds, then transient synchronization of ICs occurring approximately every 10 seconds may be a suitable marker for segmentation. Here we demonstrate the proposed segmentation method in a less obvious situation: ongoing modulation of continual manual force production.
To our knowledge, this is the first application of using transient synchronization of ICs for temporal segmentation of time series. Also, it is the first application of joint MI and network analysis to assess information transmission abnormalities between different cortical areas in PD. The main contributions of this paper are as follows:
Propose a novel EEG segmentation approach to address the non-stationary nature of EEG data especially during performance of motor tasks, and to select task-related segments.
Present a coherent MI-based network analysis framework for modeling EEG to determine dependencies between EEG channels and infer statistically significant difference between groups.
Demonstrate how the proposed framework can be used for assessing the EEG during dynamic motor behaviors in pathological conditions, such as PD.
The paper is organized as follows: the detailed discussion of the proposed framework is presented in the Methods section. The Results and Discussion section describes the EEG experimental design and summarizes the results in a real case study of PD. Finally, we summarize and conclude the paper in Conclusions.
We recorded five minutes of EEG data while subjects performed simple hand movements in order to gain insight into the difficulty that PD subjects face when performing simultaneous movements.
EEG data contain a wide range of frequency components, many of which are not of clinical or physiological interest. The data are therefore initially bandpass-filtered by a 4th order Butterworth filter between 0.5–55 Hz .
ICA Noise Removal
Here the infomax-ICA algorithm  is applied to decompose the EEG signals. ICA finds a coordinate frame in which the data projections have minimal temporal overlap by minimizing the mutual information among the data projections or maximizing the joint entropy of a nonlinear function of s. It is most appropriate to perform ICA decomposition on sources that are linearly mixed in the recorded signals without time delays. After the artifactual sources are identified, the corresponding columns of the mixture matrix (i.e. calculated as the pseudo-inverse of W) that multiply the artifactual sources are set to zero to eliminate the artifacts and thus obtain the "corrected" EEG signal. In our study, we only remove well-known, obvious artifacts by identifying 1 to 2 components (e.g. representing eye-movements and/or electrocardiac signals) by visual inspection. Failure to remove these artifacts may result in correspondence between EEG channels being falsely attributed to synchronized brain activity.
EEG Segmentation based on the Cross-Spectrogram
We note that if the derived ICs were truly independent, then the cross-spectrum would not be significant. However, in real data many of the assumptions of ICA are violated. The data are not stationary, and the time courses are not temporally white. By using infomax-ICA, which does not incorporate time delays, the derived components will be maximally independent at zero lag. As such, it will deal with the problem of volume conduction – where a deep electrical source may impart common electrical activity to two or more channels. Even though ICs are maximally independent over the whole time range, they may exhibit partial synchronization within specific time/frequency window , through which the transient coupling of neural networks might be revealed. By examining the ICs within a short moving window, the non-stationary nature of the EEG will be explored, and significant dependencies between ICs become apparent. Recent studies such as [19, 20] have also explored the transient synchrony between ICs and suggested transient correlation between ICs.
with f being the frequency normalized by the sampling frequency. By using a short-term time shifting window, we are able to obtain localized frequency contents of the two signals and their relationship with respect to time and frequency. The cross-spectrum is computed based on a short (3 s) time window shifted by 0.5 s to obtain the localized time information. Power in the higher frequency ranges, such as the gamma band (> 20 Hz) are more likely to distributed over a broad frequency range. To avoid any potential confounds from the AC current at 60 Hz, we look for sharp increases in activity in the range 45 Hz–55 Hz as a good marker for transient broadband artifacts that were not eliminated by the ICA Noise Removal step. In contrast, by examining the cross-spectrogram of pairs of ICs within the lower frequency bands of physiologic interest, we can identify task-related segments.
Mutual Information based Network
where P X (x) is the probability that x is drawn from X and P XY (x, y) is the joint probability density function for the measurements of X and Y that produce the values x and y.
where I r (X, Y) is in the range [0, 1], and H(X) and H(Y) are the entropies. Entropy H(X), defined as-Σ x P X (x)log2P X (x), is regarded as a measure of uncertainty about a random variable X.
In order to graphically represent a large set of data, we have derived both a relevance network  and an MI network. A relevance network, originally devised for graphically depicting the relationship between genes , can be generalized to take large data sets of experimental data and graphically depict the result of pair-wise MI. It is obtained by applying a threshold and only the connections that are above the threshold are displayed in the network. The relevance network can then be used in graphical theoretical analysis (discussed in Methods – Network Analysis).
In addition, we have also taken into account the magnitude of MI values and obtained an MI network from the one-way ANOVA test (with details given in Methods – Network Analysis (Statistical Analysis)). The connections are established in the MI network if their MI values exceed a specified threshold and the ANOVA tests indicate significantly different values between groups.
Graphical Theoretical Analysis on Relevance Network
with n being the total number of vertices in the graph. Such C measures the local connectivity and ranges from 0 to 1. The higher the C, the greater the intensity of connections within a cluster.
The L of a graph is the mean of all shortest paths (shortest distance) connecting all pairs of vertices. It has a value greater than 1 and measures the global connectivity of the graph. A detail graphical explanation of a graph and graph theoretical measures can be found in .
In addition to the graphical theoretic analysis on the thresholded matrices, we created an "MI network" which also incorporated the magnitude of the MI values. As before, MI values were first thresholded by zeroing values less than the 95th percentile on a null distribution. A null distribution was obtained by repeatedly (n = 100) randomly permuting the order of the second signal and computing the pair-wise MI based on them . The MI differences between segments are analyzed by one-way ANOVA with subject number, groups, and tasks as factors . The connections between any two channels are established for the MI network if they are significantly different according to the ANOVA test and have magnitudes that are greater than the permutation threshold. The normality of the distribution of the MI values is verified by the Kolmogorov-Smirnov (KS) test .
In the one-way ANOVA test for each pair-wise MI I(X, Y), the effect of a factor (e.g. Group) is tested, by comparing with the F-test the variance of I(X, Y) explained by the factor against the variance of the residuals. Consequently, a p-value was calculated for each possible connection in the MI network. To account for the effect of testing multiple connections simultaneously, the p-values are corrected for multiple hypothesis testing using Storey's positive-false-discovery-rate (pFDR) procedure  which computes a q-value, the expected ratio of falsely rejected hypotheses among all those being rejected. Connections whose q-values were smaller than 5% are considered statistically significant.
Because we are more interested in the connection with greater MI values, the permutation test is used in conjunction with the ANOVA test to select the relevant features for the MI based network. We have chosen the largest observed value of the permutation test as our threshold. A connection is thus established in the MI network if the MI values are significantly different based on the q-value and are above the maximum observation from the permutation test.
Results and Discussion
Subjects and Experiment Design
All research was approved by the University of British Columbia Ethics Board. After giving informed consent, seven PD and six age-matched control subjects volunteered to participate in the study. All patients were diagnosed with mild to moderately severe PD (Hoehn and Yahr stage 1–3) . The control subjects were confirmed to be without active neurological disorders by a qualified neurologist. All patients were taken off L-Dopa medication after overnight withdrawal of at least 12 hours.
EEG Data Preprocessing
Noisy EEG Segment Removal
Task-Related EEG Segmentation
The amount that subjects were asked to squeeze was based on two sinusoids with period of 10 and 18 seconds. For the BO condition, the color change occurred every 20 seconds. Therefore, autocorrelations of the cross-spectrogram of the ICs over the three physiologically-relevant frequency bands  5–8 Hz (Theta), 8–12 Hz (Alpha), 12–30 Hz (Beta) that have a peak at 10 seconds or 18–20 seconds are selected for sections of EEG segmentation. Depending on the features of each dataset, approximately five pairs were chosen for each task. Only segments that are above the mean plus the mean absolute deviation are considered as task-related and obtained for further analysis.
Mutual Information based Network Analysis
Graphical Theoretical Analysis Applied to the Relevance Network
The MI matrix for each subject is converted to a graph separately, and the means of the cluster coefficient C and the shortest path length L of the graph within the group (ie. N-SQ, PD-SQ, N-BO, PD-BO) were computed as a function of the threshold T.
As T is varied from 0.01 to 0.3, the graphs started to break into subgraphs. In addition, at higher T, some subjects start to have empty graphs meaning the graphs contain no connection at all. Therefore, when we interpret the results, we need to make note of where those points are and they are summarized in Table 1. Again, we see that the points between N and PD do not differ much because their means are very close. The overall mean C of the graph for each group as a function of T was computed and compared. Because the means of the four groups (ie. N-SQ, PD-SQ, N-BO, PD-BO) as a function of T follow the same pattern and we are more interested in the differences between groups, the deviation from the overall mean of the four groups as a function of T is illustrated at the top panel of Fig. 11. The overall mean of the four group as a function of T is displayed at the bottom left corner of the top panel. The bottom panel shows the region that is significantly different between groups (at 2: N vs. PD during SQ; at 3: N vs. PD during BO; at 4: SQ vs. BO for N; at 5: SQ vs. BO for PD). The significance of the C between groups is tested with the pair-wise t-test (p < 0.05). In general, the intensity of connections within the cluster does not differ significantly between tasks (SQ and BO) within a group. The intensity of connections within the cluster is higher for PD compared to N for SQ task and BO task, especially at higher threshold T (T > 0.15). The consistently lower C seen in N compared to PD across all frequency bands probably reflects the widespread, excessive synchronization seen in PD [34, 35]. Unlike prior studies that have emphasized synchronization in the beta band, we have observed excessive synchronization in all bands, especially the theta band.
Squeeze Task (SQ) and Both Task (BO): Threshold T when graphs start to split into subgraphs or become empty graphs
Squeeze Task (SQ): Clustering coefficient C per vertex and the p-value of the pair-wise T-test between Normal and Parkinson's subjects.
Both Task (BO): Clustering coefficient C per vertex and the p-value of the pair-wise T-test between Normal and Parkinson's subjects. (T = 0.2 for 5–8 Hz; T = 0.18 for 8–12 Hz; T = 0.15 for 12–30 Hz)
Statistical Analysis on MI Network
This paper proposed a novel segmentation, mutual information network framework for EEG connectivity analysis for subjects performing a motor task. The greatest EEG changes during motor performance are typically event-related synchronization/desynchronization: EEG responses are not phase-locked to motor performance, but rather tend to be associated with augmentation or attenuation of specific frequency bands. This means that standard methods of averaging the EEG time-locked to the motor performance will tend to be inaccurate, necessitating the use of alternate methods. In addition, ERS/ERD is typically investigated in univariate fashion, where each EEG channel is examined independently for altered localized neuronal synchrony resulting in changes in the frequency spectra at that channel. Here we have used the cross-spectrum of ICs as a marker for segmentation. This has multiple benefits: first, consistent with ERS/ERD, it examines the data in the frequency domain, second, it allows the examination of multiple channels simultaneously, as each IC will consist of a linear weighting of all channels, and lastly, it will allow unmixing of the raw data so that task related activity can be extracted from ongoing background brain rhythms. After segmentation, we used mutual information to measure both linear and nonlinear dependencies, without assuming strict phase locking of signals. This allowed us to create a relevance network suitable for graph-theoretic analysis methods, and an MI network which further incorporated the magnitude of the MI at each channel pair to allow statistical analysis with ANOVA.
The proposed method provided several novel insights into abnormalities in PD subjects. The well known clinical observation of difficulty in performing simultaneous movements  appears to be related to an inability to recruit different brain areas independently. When normal subjects were asked to perform tasks with two hands compared to a task with one hand, there was no significant difference in the theta bands, and only mild changes in connectivity in the alpha and beta bands (Fig. 16). In contrast, when PD subjects attempted to perform simultaneous movements, they appeared unable to recruit different brain areas independently, resulting multiple areas of synchronization in the theta range, and also in the beta range. These results appear novel, as previous research has emphasized excessive synchronization in the beta range in PD .
Additionally, the widespread synchrony that is normally already seen during regular unimanual or bimanual performance was of a different form in PD (Figs. 11, 12, and 17). PD subjects tended to have higher cluster coefficients, C, and lower shortest path lengths, L suggesting focal clusters of synchronous activity. Taken together, these results suggest that normal subjects can synchronously activate broad areas or cortex independently in response to differing task demands. In contrast, PD subjects appear to have islands of hypersynchronicity that cannot be recruited independently. This is consistent with known biology of PD, where bradykinesia is most likely the result of "noisy" basal ganglia input to the frontal cortex and appears to critically depend upon dopamine depletion  as would be seen in the PD subjects off medication in this study. The higher MI values in the frontal region at lower and medium frequency bands and motor cortex at higher frequency in PD (Fig. 17) are also consistent with the cortical regions that receive output from the basal ganglia. Since the EEG is normally assumed to reflect cortical activity, isolation of clear abnormalities from PD subjects who have predominantly basal ganglia dysfunction and preserved cortical function is notable.
Although we suggest that our results demonstrate strong evidence to support the proposed framework as a tool to study EEG signals, there are nevertheless limitations of the proposed method. For example, currently only pair-wise MI, one particular case of calculating the MI between M random variables, is investigated. In a more general presentation, the corresponding M–dimensional MI can be defined. Since pairwise independence does not necessarily imply global independence, M–dimensional MI may reveal additional information from that of pair-wise MI, and thus may be a fruitful avenue to explore in the future for EEG analysis. Similarly, pair-wise MI may suffer from a high false discovery rate, i.e. nodes are erroneously associated while in truth they only indirectly interact through one or more other nodes. Therefore, to prune the reconstructed network of such false positives, we can extend the current work by exploring the concept of conditional mutual information (CMI) instead. Also, in the current approach, we did not consider the temporal information embedded in the time-series EEG data. As one future work, we intend to introduce temporality into the proposed MI network construction.
This work was supported by the Canadian Natural Sciences and Engineering Research Council (NSERC) under grant STPGP 365208-08, and by the CHRP grant.
- Rowe J, Stephan KE, Friston K, Frackowiak R, Lees A, Passingham R: Attention to action in Parkinson's disease: Impaired effective connectivity among frontal cortical regions. Brain 2002,125(2):276–289. 10.1093/brain/awf036View ArticleGoogle Scholar
- Perda E, Quiroga RQ, Bhattacharya J: Nonlinear multivariate analysis of neurophysiological signsl. Progress in Neurobiology 2005, 77: 1–37. 10.1016/j.pneurobio.2005.10.003View ArticleGoogle Scholar
- Shen B, Nadkarni M, Zappulla RA: Spectral modulation of cortical connections measured by EEG coherence in human. Clinical Neurophysiology 1999, 110: 115–125. 10.1016/S0013-4694(98)00104-7View ArticleGoogle Scholar
- Weiss S, Rappelsberger P: Long-range EEG synchronisation during word encoding correlates with successful memory performance. Cognitive Brain Ressearch 2000, 9: 299–312. 10.1016/S0926-6410(00)00011-2View ArticleGoogle Scholar
- Blandini F, Nappi G, Tassorelli C, Martignoni E: Functional changes of the basal ganglia circuitry in Parkinson's disease. Progress in Neurobiology 2000, 62: 63–88. 10.1016/S0301-0082(99)00067-2View ArticleGoogle Scholar
- Samuel M, Ceballos-Baumann AO, Blin J, Uema T, Boecker H, Passingham RE, Brooks DJ: Evidence for lateral premotor and parietal overactivity in Parkinson's disease during sequential and bimanual movements. A PET study. Brain 1997,120(6):963–976. 10.1093/brain/120.6.963View ArticleGoogle Scholar
- Stam CJ, van Dijk BW: Synchronization likelihood: an unbiased measure of generalized synchronization in multivariate data sets. Phys D 2002, 163: 236–251. 10.1016/S0167-2789(01)00386-4MathSciNetView ArticleGoogle Scholar
- Hurtado JM, Rubchinsky LL, Sigvardt KA, Wheelock VL, T PC: Temporal Evolution of Oscillations and Synchrony in GPi/Muscle Pairs in Parkinson's Disease. Journal of Neurophysiology 2005, 93: 1569–1584. 10.1152/jn.00829.2004View ArticleGoogle Scholar
- Jeong J, Gore JC, Peterson BS: Mutual information analysis of the EEG in patients with Alzheimer's disease. Clinical Neurophysiology 2001, 112: 827–835. 10.1016/S1388-2457(01)00513-2View ArticleGoogle Scholar
- Na SH, Jin SH, Kim SY, J HB: EEG in schizophrenic patients: mutual information analysis. Clinical Neurophysiology 2002, 113: 1954–1960. 10.1016/S1388-2457(02)00197-9View ArticleGoogle Scholar
- Min BC, Jin SH, Kang IH, Lee DH, Kang JK, Lee ST, Sakamoto K: Analysis of Mutual Information Content for EEG Responses to Odor Stimulation for Subjects Classified by Occupation. Chem Senses 2003, 28: 741–749. 10.1093/chemse/bjg066View ArticleGoogle Scholar
- Benecke R, Rothwell JC, Dick JP, Day BL, Marsden CD: Performance of simultaneous movements in patients with Parkinson's disease. Brain 1986, 109: 793–757. 10.1093/brain/109.4.739View ArticleGoogle Scholar
- Berardelli A, Rothwell JC, Thompson PD, Hallett M: Pathophysiology of bradykinesia in Parkinson's disease. Brain 2001, 124: 2131–2146. 10.1093/brain/124.11.2131View ArticleGoogle Scholar
- Kaplan AY, Fingelkurts AA, Fingelkurts AA, Borisov SV, Darkhovsky BS: Nonstationary nature of the brain activity as revealed by EEG/MEG: Metholdological, practical and conceptual challenges. Signal Processing 2005, 85: 2190–2212. 10.1016/j.sigpro.2005.07.010View ArticleGoogle Scholar
- Sporns O, Zwi JD: The small world of the cerebral cortex. Neuroinformatics 2004, 2: 145–162. 10.1385/NI:2:2:145View ArticleGoogle Scholar
- Stam CJ, Jones BF, Nolte G, Brekspear M, Scheltens p: Small-world Networks and functional connectivity in Alzheimer's Disease. Cerebral Cortex 2007, 17: 92–99. 10.1093/cercor/bhj127View ArticleGoogle Scholar
- Rappelsberger P, Pfurtscheller G, Filz O: Calculation of eventrelated coherence A new method to study short lasting coupling between brain areas. Brain Topography 1994,7(2):121–127. 10.1007/BF01186770View ArticleGoogle Scholar
- Leocani L, Toro C, Manganotti P, Zhuang P, Hallett M: Eventrelated coherence and eventrelated desynchronization/synchronization in the 10 Hz and 20 Hz EEG during self-paced movements. Electroencephalography and Clinical Neurophysiology/Evoked Potentials Section 2007,104(3):199–206. 10.1016/S0168-5597(96)96051-7View ArticleGoogle Scholar
- Gupta D, James CJ: Narrowband vs. Broadband Phase Synchronization Analysis Applied to Independent Components of Ictal and Interictal EEG. IEEE-EMBS Proceedings 2007.Google Scholar
- Hong B, Acharya S, Gao S, Thakor NV: Transient phase synchrony of independent cognitive components underlying scalp EEG. IEEE-EMBS Proceedings 2005.Google Scholar
- Pfurtscheller G, da Silva FHL: Event-related EEG/MEG synchronization and desynchronization: basic principles. Clinical Neurophysiology 1999, 110: 1842–1857. 10.1016/S1388-2457(99)00141-8View ArticleGoogle Scholar
- Jung TP, Makeig S, Humphries C, Lee TW, McKeown ML, Iragui V, Sejnowski TJ: Removing electroencephalographic artifacts by blind source separation. Psychophysiology 2000, 37: 163–178. 10.1017/S0048577200980259View ArticleGoogle Scholar
- Bell A, Sejnowski T: An information-maximization approach to blind separation and blind deconvolution. Neural Computation 1995, 7: 1129–1159. 10.1162/neco.19126.96.36.1999View ArticleGoogle Scholar
- Delorme A, Makeig S: EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. Journal of Neuroscience Methods 2004, 134: 9–21. 10.1016/j.jneumeth.2003.10.009View ArticleGoogle Scholar
- Oppenheim AV, Schafer RW: Discrete-Time Signal Processing. New Jersey: Prentice-Hall; 1999:737.Google Scholar
- Cassidy MJ, Brown P: Hidden Markov based autoregressive analysis of stationary and non-stationary electrophysiological signals for functional coupling studies. Journal of Neuroscience Methods 2002, 116: 35–53. 10.1016/S0165-0270(02)00026-2View ArticleGoogle Scholar
- Chiang J, Wang Z, McKeown MJ: A Hidden Markov, Multivariate Autoregressive (HMM-mAR) Network Framework for Analysis of Surface EMG (sEMG) Data. IEEE Trans on Signal Processing 2008,56(8):4069–4081. 10.1109/TSP.2008.925246MathSciNetView ArticleGoogle Scholar
- Cover TM, Thomas JA: Elements of information theory. New York: John Willey & Sons; 1991.View ArticleGoogle Scholar
- Butte AJ, Kohane IS: Mutual information relevance networks: functional genomic clustering using pairwise entropy measurements. Pacific Symposium on Biocomputing 2000, 5: 415–42.Google Scholar
- Francois D, Wertz V, Verleysen M: The permutation test for feature selection by mutual information. ESANN Proceedings 2006, 26–28.Google Scholar
- Hogg RV, Ledolter J: Engineering statistics. New York: MacMillan; 1987.Google Scholar
- D SJ: A direct approach to false discovery rates. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2002,64(3):479–498. 10.1111/1467-9868.00346MathSciNetView ArticleGoogle Scholar
- Hoehn M, Yahr M: Parkinsonism: onset, progression and mortality. Neurology 1967, 17: 427–442.View ArticleGoogle Scholar
- Salenius S, Avikainen S, Kaakkola S, Hari R, Brown P: Defective cortical drive to muscle in Parkinson's disease and its improvement with levodopa. Brain 2002, 125: 491–500. 10.1093/brain/awf042View ArticleGoogle Scholar
- Timmermann L, Gross J, Dirks M, Volkmann J, Freund HJ, Schnitzler A: The cerebral oscillatory network of parkinsonian resting tremor. Brain 2003, 126: 199–212. 10.1093/brain/awg022View ArticleGoogle Scholar
- Chen C, Litvak V, Gilbertson Te: Excessive synchronization of basal ganglia neurons at 20 Hz slows movement in Parkinson's disease. Experimental Neurology 2007, 205: 214–221. 10.1016/j.expneurol.2007.01.027View ArticleGoogle Scholar
- Rivlin-Etzion M, Marmor O, Heimer G, Raz A, Nini A, Bergman H: Basal ganglia oscillations and pathophysiology of movement disorders. Current Opinion in Neurobiology 2006,16(6):629–637. 10.1016/j.conb.2006.10.002View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.