- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2012 (2012), Article ID 239210, 13 pages
Inferring Functional Neural Connectivity with Phase Synchronization Analysis: A Review of Methodology
1Med-X Research Institute, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China
2Department of Automation, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
Received 25 November 2011; Accepted 31 January 2012
Academic Editor: Dimitris Kugiumtzis
Copyright © 2012 Junfeng Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Functional neural connectivity is drawing increasing attention in neuroscience research. To infer functional connectivity from observed neural signals, various methods have been proposed. Among them, phase synchronization analysis is an important and effective one which examines the relationship of instantaneous phase between neural signals but neglecting the influence of their amplitudes. In this paper, we review the advances in methodologies of phase synchronization analysis. In particular, we discuss the definitions of instantaneous phase, the indexes of phase synchronization and their significance test, the issues that may affect the detection of phase synchronization and the extensions of phase synchronization analysis. In practice, phase synchronization analysis may be affected by observational noise, insufficient samples of the signals, volume conduction, and reference in recording neural signals. We make comments and suggestions on these issues so as to better apply phase synchronization analysis to inferring functional connectivity from neural signals.
Segregation and integration are two fundamental organization principles of neural systems in brain . The neural organization can be investigated with functional neural connectivity in both local and global brain regions. The connectivity in local regions reflects specialized functions of local cortex, while the connectivity among spatially separated brain regions plays important roles in advanced cognitive function [2–7]. Various measures, such as phase synchronization index (PSI), mutual information, correlation coefficient, coherence, and partial directed coherence, have been applied to inferring neural connectivity among brain units at multiple temporal and spatial scales with neural signals including electroencephalography (EEG), magnetoencephalography (MEG), functional magnetic resonance imaging, and local field potential [8–13]. These measures can be classified into several families such as correlation/coherence family, phase synchronization family, and Granger causality family. The measures in the same family usually yield results with strong correlation to each other and thus could provide little additional information on neural connectivity, while some measures belonging to different families may have weak correlation to each other in inferring neural connectivity, which implies that they each could characterize interdependence of signals in different aspects .
Among these measures, PSI quantifies the relationship between instantaneous phases (IP, represents the rhythm of oscillation or signal wave) of coupled systems/brain units but neglects the effect of their amplitude. PSI has been demonstrated to be effective in inferring neural connectivity, especially when the connectivity is too weak to be detected by other measures [8, 14]. Note that there are other types of synchronization, such as complete synchronization and generalized synchronization [14, 15]. But most of them are defined for theoretical models of coupled oscillators and difficult to be applied to neural signal analysis. The coupled systems/units are claimed to be in phase synchronization (PS) when the difference of IPs is bounded with respect to time. PS is ubiquitous in both natural and man-made systems, such as neural oscillations [4, 16], coupled chaotic oscillators , chaotic laser arrays , and electrochemical oscillations . In addition, two metrics of PS, that is, the phase-lock interval and the lability of global synchronization, hold power law probability distributions for both simulation systems (i.e., the Ising model and the Kuramoto model) and brain networks at a broadband frequency. These results imply that the IP association among multiple units or brain regions hold the property of criticality .
PSI, also called phase-locking value in literature, is proposed to quantify the level of PS and has been applied to a wide range of neuroscience research [3, 5]. For example, PSI has been used to examine the alternation of cortical connectivity prior to seizures  and neuronal synchrony after ischemic stroke , to identify mental states in brain-computer interface [23, 24], to investigate cognitive dysfunctions of mental disorder , and to gain new strategies for clinical treatment [22, 26]. However, to get a reliable estimation of PSI from observed neural signals is not so easy, especially when the signals are with a small number of samples and/or contaminated by noise [10, 11, 27, 28]. Note that, in literature, other measures, such as phase clustering index, used different definition of phase, which is not within the framework of IP [29, 30]. In this paper, we review the methodology of PS analysis in inferring functional neural connectivity from the following aspects.IP definition. To detect PS in a pair of signals, various IP definitions have been proposed. Most of these IP definitions are based on particular transforms, such as the Hilbert transform  and the wavelet transform . But actually, these IP definitions can be unified into one framework which defines IP with specific filter applied to signals .PSI and its significance test. We will introduce two commonly used PSIs, which are based on entropy [2, 33] and circular statistics [34, 35], respectively. We will further introduce several strategies to provide significance test for the estimated PSI.Practical problems in PS detection. The estimation of PSI is affected by observational noise, volume conduction, the number of samples in observed signals, reference, and other factors. We will summarize the advances on these issues and give suggestions for a better PS detection.Extensions of PS analysis. Lots of extensions/variations of PS analysis have been proposed for multitrial signals and multivariate signals. We will give a brief overview of these methods and some comments to their applications in neural signal analysis.
2. Definition of Instantaneous Phase
Before estimating the level of PS, the IP of each signal should be defined (Figure 1). The most basic definition of IP is based on the Hilbert transform [14, 17, 36], which can be directly applied to coherent signals (i.e., narrow band signals with one prominent spectral component). But, for noncoherent data such as raw neural signals, this IP definition may yield negative instantaneous frequency (IF, defined as the derivative of IP with respect to time), which is physically meaningless [37–39]. One way is to process noncoherent signal with a specific narrow bandpass filter. An alternative way is to define IP based on other transforms such as the wavelet transform [18, 19, 31, 40]. PS detection based on these IP definitions has been compared numerically with both simulation data and experimental signals, yielding similar connectivity [10, 11]. In addition, an analytical study on these IP definitions has unified them into one framework which defines IP with a specific bandpass filter applied . In this section, we will briefly introduce issues on IP definition based on results in .
2.1. Definitions of Instantaneous Phase
2.1.1. IP Definition Based on the Hilbert Transform
For a given signal , its analytic signal is defined as where is the Hilbert transform of (here P.V. means that the integral is taken in the sense of the Cauchy principal value). Then the IP of is defined as In the frequency domain, appears as , where is the Fourier transform of and is the Fourier transform of . Then in numerical implementation, the estimation of analytic signal can be obtained by the inverse Fourier transform of . In the time domain, can be expressed as the convolution of with the complex-response filter ; that is,
2.1.2. IP Definition Based on Gaussian Filter
Another method defines IP as where is the convolution of with a narrow-band Gaussian filter , which is shifted by the nominal frequency . In the frequency domain, can be expressed as , where . PS analysis based on this IP definition has successfully detected PS in coupled laser arrays, which was not revealed by IP defined with the Hilbert transform . The possible reason is that the Gaussian filter extracts the components of the laser data in a particular frequency band which are in PS, while PS analysis based on the Hilbert transform does not use filter to extract the components, and thus the underlying PS level is submerged by noise and components in other bands.
2.1.3. IP Definition Based on the Wavelet Transform
With the Gabor wavelet , IP is defined as where is the convolution of with , and is the envelope . The difference between IP definitions based on the Gaussian filter and on the wavelet transform lies in that the filters and are with Gaussian windows of different amplitudes and width; that is, is scaled by . In the frequency domain, appears as , and is .
2.1.4. Unified Framework for IP Definitions
The IP definitions discussed above all can be expressed as the argument of the signal after a specific filter. In other words, these IP definitions can be unified into one common framework ; that is, where Here is a filter with envelope , nominal frequency , and response duration . The filter assures that the extracted signal wave is coherent. is analytic in an asymptotic sense as the bandwidth of is smaller than [36, 41, 42]. This method is exactly the one based on a Gaussian filter when the envelope is a Gaussian function. Filter can be expressed as in the frequency domain, where is the Fourier transform of .
The analytic signal can be interpreted as a combination of the Hilbert transform and a real bandpass filter. Let , where is the real part of . Then appears as in the frequency domain. With (4), the analytic signal of can be obtained by performing the inverse Fourier transform to ; that is, where denotes the inverse Fourier transform operator.
Beside IP definitions within this framework, there are IP definitions proposed in other aspects [43, 44]. For example, the IP and PSI based on Rihaczek distribution are more robust to noise than the method based on the wavelet transform .
2.2. Constraints for IP Definition
For a real signal , its imaginary counterpart is usually unobservable, where and denote the real and the imaginary parts of the complex variable , respectively. Generally, the imaginary part is estimated from by a certain operation; that is, . Among various operators proposed, the Hilbert transform is the only one that satisfies the following three conditions .The associated amplitude is continuous and differentiable.The IPs of signals and is the same; that is, the operator should possess the property , where is a constant.For any constant amplitude , frequency , and phase , the operator satisfies .
In addition, the Bedrosian theorem gives more constraints on IP definition based on the Hilbert transform [36, 42]. This theorem states that the relation holds for a low-frequency term and a high-frequency term whose spectra do not overlap in the frequency domain. For real signal , is its quadrature model. It seems that in this model is a good candidate for definition of IP. However, it is difficult to estimate and from only the observable signal with this model. There is a difference between and the analytic signal (1). This difference approaches zero in an asymptotic sense as and fulfill the Bedrosian theorem . Note that the Bedrosian theorem sets a constraint to the signal, which is similar to the second condition; that is, , that set on the operator. For noncoherent signal, bandpass filter is usually applied, so that the extracted signal wave is coherent and fulfills the Bedrosian theorem. The effective bandwidth of filter is , which should be less than . As the filter becomes narrower close to delta function (i.e., ) in the frequency domain, the IF of the components in the pass band approaches the nominal frequency in an asymptotic sense .
Generally, to infer functional connectivity with PS analysis, we would recommend to define IP by combining the Hilbert transform with specific bandpass filter; that is, a bandpass filter is first applied to extracting the neural signal waves in specific frequency band, and then IP is defined based on the Hilbert transform.
3. Phase Synchronization Analysis
3.1. Phase Synchronization Indexes
Let and denote the cumulative IP of signals observed from two coupled units, respectively. Then the signal pair is said to be in : PS when the inequality . holds, where and are positive integers. In this paper, we focus on the case of 1 : 1 PS, which is defined based on samples of one trial in the time domain. Most conclusions for 1 : 1 PS can be easily extended to the case of : PS [47, 48]. In neuroscience research, the estimated PSI is usually taken as one kind of functional connectivity in neural signals. Here we introduce two PSIs which have been commonly used to quantify the functional connectivity of neural signals. They are based on entropy [2, 33] and circular statistics [34, 35], respectively. Let denote , that is, the observation of at time , and let denote the estimation of at time , where is the sampling interval.The entropy-based PSI is estimated by where is the entropy of distribution , , , and is the number of bins of distribution [2, 33].The mean phase coherence (MPC) of IP difference is another commonly used PSI. It is defined as and can be estimated via where is the estimation of and is the number of samples.
These two PSIs quantify how concentrated the distribution of IP difference is. Their values are in [0, 1], with PSI = 1 implying that the signal pair is with exact rhythm locking and PSI = 0 indicating no PS at all. Beside these two indexes, other measures, such as the index based on conditional probability  and the index based on mutual information , have also been applied to the quantification of relationship between IPs. More discussions on PSI can be found in .
3.2. Significance Tests for PSI
In inferring neural connectivity, a highly concerned question is whether the estimated PSI could indicate that the signal pair is with significant connective strength or not. This question is also concerned in brain network research, as brain network is constructed by setting edges between nodes (i.e., brain regions) with connective strength over a predefined threshold. Accumulated studies show that both functional and structural brain networks is of small-worldness with sparse edge number [28, 50–55]. The topology of brain networks are dependent on the number of edges, that is, the number of connectivity with significantly big strength. One way to decide the value of this common threshold is based on significance test for estimated PSI [8, 34, 49, 56]. Here we introduce three different strategies in providing significance threshold for PS analysis.
3.2.1. Artificial Surrogate Tests
The first strategy is based on artificial surrogate data. Surrogate methods usually produce artificial data by randomizing the concerned property but keeping as much as possible other properties of the original signal [49, 57]. Then whether the original signal possesses the concerned property could be tests with the artificial surrogate data. Various surrogate methods have been proposed [8, 57–61]. For example, surrogate method based on autoregressive model generates surrogate realizations by fitting an autoregressive model of the original signals and using independent white noise as the inputs to the model. The surrogate data generated by this method are linear stochastic processes and have the same power spectra of the original signals [57, 60]. In another study, a method named twin surrogate is proposed based on the recurrence properties. It is demonstrated to be suitable to provide significance test for PS in Rössler oscillators .
In , four different surrogate methods, which generate artificial surrogate data by shuffling the rank order, the phase spectra, the IP of original EEG signals, are compared in significance test for PS analysis. Results show that the phase-shuffled surrogate method is workable for PS analysis. This method generates surrogate data by shuffling the phase spectra of original signal but keeping the amplitude spectra unchanged. Let denote the discrete Fourier transform of original signal ; that is, Then a surrogate realization of is generated by the inverse discrete Fourier transform of ; that is, where with a uniform random sequence [57, 58].
A level of significance corresponds to a probability of a false rejection. A larger number of surrogate realizations could offer a greater power in discrimination . In , 99 surrogate pairs are generated for each original EEG signal pair. Then if the PSI of one original signal pair is larger than the 5th biggest in all the 100 PSIs (i.e., the PSIs of the original signal pair and its 99 surrogate pairs), the original signal pair is claimed to be in PS with a 95% level of significance. Figure 2 shows the histogram of the MPC-based PSIs for 435 original EEG signal pairs and their phase-shuffled surrogate pairs for one subject, where one common significance threshold is estimated for all 435 original signal pairs in each case.
3.2.2. Analytical Significance Test
The second strategy is based on the distribution test of IP series [34, 56]. When the distribution of PSI is well approximated by an analytical model, a significance threshold then could be offered analytically. Empirical distributions of IP difference of coupled Rössler systems have been tested under the assumption that IP obeys specific distributions and the IPs of different samples are independent . But this is applicable only for special cases, as the distribution of IP difference may vary with respect to different dynamical systems. In another study, the IP difference is formulated as an increment process under the assumption that the increments could be represented by an -mixing process. A theoretical significance level is proposed for the MPC-based PSI. Simulation results show that the significance level is workable when the time series is with enough samples and two dependent parameters could be reliably estimated .
3.2.3. Surrogate Test Based on Intersubject Signal Pairs
The third strategy is based on the assumption that the neural signals of different subjects are independent. In , intersubject EEG signal pair denotes two EEG signals from two different subjects, and intra-subject EEG signal pair denotes two EEG signals from the same subject. Then under this assumption, the intersubject EEG signal pairs can be used as surrogate pairs for intra-subject EEG signal pairs. Compared with above artificial surrogate data, intersubject surrogate pairs seem to possess more inherent features of EEG signals but with no association between each other. Then significance test based on intersubject EEG signal pairs may be used as a standard criterion in evaluating the performance of artificial surrogate methods. Significance test with this strategy shows that the histogram of the MPC-based PSIs for intersubject surrogate pairs and intra-subject EEG signal pairs is similar to that based on the phase-shuffled surrogate method (Figure 2) , and the significance threshold suggested by intersubject surrogate pairs is close to that by the phase-shuffled surrogate method. This implies that the phase-shuffled surrogate method is workable in providing significance test for PS analysis. In our opinion, more studies on significance test for PS analysis are needed, and the studies combined the three different strategies introduced above would be promising when they may reach consistent results.
4. Problems in Phase Synchronization Detection
In real applications, the observed neural signals are more or less contaminated by noise, and the samples collected are usually limited. Obviously, the observational noise will degrade the estimation of PSI and even submerge non-trivial PS in neural signals. In addition, the MPC-based PSI is a biased estimator, which implies that the reliability of functional connectivity inferred by it will decrease as the samples collected in signals are insufficient. To get a reliable PS detection in real applications, we must take the issues into consideration. In this section, we would discuss challenges of PS analysis in practice. In particular, we would introduce the advances on the effect of noise in PS detection, the influence of signal duration and estimation bias of PSI, the effect of volume conduction, and the influence of reference in PS analysis.
4.1. Effect of Noise in PS Detection
The effect of noise in PS detection has been examined by both numerical computation [10, 11, 27, 33, 62, 63] and analytical study . A bandpass prefiltering can suppress the effect of noise but may introduce spurious connectivity as well . In this section, we mainly review a theoretical study on the effect of noise in IP estimation and PS detection based on the unified framework of IP definition .
Let denote the noisy signal, where is the clean signal and is the additive noise. The analytic signal of can be defined with a bandpass filter ; that is, where . Let denote the estimate error of IP due to the noise term , where denotes the estimate of from the noisy signal . It has been proved that the distribution of is a normal distribution; that is, , when the instantaneous signal-to-noise ratio, defined as in the pass band, is larger than 5. Here is the variance of , and . The wrapped , that is, , obeys the wrapped normal distribution, [35, 64],
For the MPC-based PSI, , the effect of noise appears as that is, the noise introduces a degrading factor, , into the true index . This factor is determined by only the level of in-band noise. If in (22) is not a constant, actually obeys a conditional distribution, and (22) turns out to be For observed signal , the distribution of IP error is a scale mixture of normal distributions (SMN) with different variances. If the probability density function of is known, the empirical distribution of IP error can be approximated as where are the respective empirical probabilities which are estimated from . Note that is the instantaneous amplitude of the clean signal, and thus it is difficult to obtain its distribution with only the observed noisy signal . Simulations are performed by considering the SMN of the IP error as a normal distribution with constant standard deviation ; that is, , and computational outputs support the theoretical results (24). When the in-band noise level is not so high (signal-to-noise ratio greater than 10 dB), the estimated PSI is not much affected by observational noise .
4.2. Influence of Signal Duration and Estimation Bias
Another question in PS analysis is how long a epoch of signals is enough to get a reliable quantification of neural connectivity? One microstate of EEG signals usually lasts 50 to 200 ms [3, 9, 66]. However, various durations (from 100 ms to 10 ms) of EEG signals have been used in examining functional connectivity [12, 24, 67, 68]. A comprehensive study has been performed regarding this issue based on both surrogate tests and intersubject EEG surrogate test . Results show that a duration of EEG waves of wave cycles is suitable for PS analysis. The value of PSI is not only determined by the relationship between two IP sequences but also highly affected by the duration of signals used for PSI estimation. Too short duration of signals will result in too large PSI, while too long duration of signals will yield too small PSI, both of which may not offer a good indication on the level of PS (Figure 2) .
In another study, the IP difference of a signal pair with no synchrony is assumed to obey uniform distribution. Under this assumption, the expected value of the MPC-based PSI is where is the number of samples used for PSI estimation [43, 69]. This means that the MPC-based PSI is a positively biased estimators for finite sample sizes . More samples used will result in a less biased PSI. Then even for the same duration of signals, signals observed with higher sampling rate will yield a better estimation of the MPC-based PSI than those with lower sampling rate. The effect of sampling rate on both the MPC-based PSI and the entropy-based PSI has been evaluated. Results demonstrate that the MPC-based PSI is less affected by sampling rate than the entropy-based PSI and thus is recommended for PS analysis of EEG signals measured with low sampling rate .
In many experimental studies, the number of samples or trials of measured data is limited. Bias will be unavoidable in PSI estimation due to insufficient samples. In addition, even when the observed signals are with sufficient samples, they may be nonstationary. In this case, we still want to estimate PSI with samples in short time window, so as to obtain functional connectivity with a certain degree of temporal resolution. Regarding this problem, a new measure called pairwise phase consistency has been proposed . This measure is defined as the mean of the cosine of IP difference across all given signal pairs; that is, Pairwise phase consistency quantifies the similarity of relative IP among trials or samples and is bias-free. In addition, it is demonstrated that pairwise phase consistency is equivalent to the squared MPC-based PSI in population statistic .
4.3. Influence of Volume Conduction
The activities of a single source within the brain can be observed by many sensors on the scalp . This is usually referred as volume conduction. Then the PSI between signals measured by different sensors, especially by spatially adjacent sensors, may be a trivial artefact due to volume conduction, but not a true interaction of the underlying brain activities . There are two ways to tackle this problem. One way is to estimate the true sources of underlying brain activities with observed EEG/MEG signals using inverse method and then quantify the relationship of the estimated sources rather than the relationship of the observed signals. However, the true sources are usually not easy to be obtained by inverse methods.
An alternative way is based on the imaginary part of coherency of observed signals. We assume that two signals and collected by sensors and are from a linear superposition of independent sources , and the mapping of each source to sensors is instantaneous with no distortion . Then in the frequency domain, can be expressed as where and are the Fourier transform of and , respectively, and are the contribution coefficient of to . The cross-spectrum of and is where denotes the Kronecker-delta function. Here is real, which implies that volume conduction of multiple sources strongly affects the real part of the cross-spectrum between and but does not affect the imaginary part [75, 76].
The complex coherency of and is defined as In , the imaginary part of is defined as a synchronization measure for and in the frequency domain. The mean of over all frequencies is equal to the mean of the cross-correlation of the corresponding analytic signals over time ; that is, where and denote the average, respect to frequency and time respectively. Then the imaginary part of is However, it has been demonstrated that the imaginary part of is not a good index of PS as it depends on the amplitudes of signals [76, 77]. Further, a measure called phase lag index, was proposed based on the consideration that “the existence of a consistent, nonzero phase lag between two times series cannot be explained by volume conduction from a single strong source” . Recently, it is demonstrated that the performance of phase lag index could be degraded by small perturbations . To deal with this problem, a weighted phase lag index is defined as Compared with phase lag index, the weighted phase lag index is demonstrated with “reduced sensitivity to additional, uncorrelated noise sources and increased statistical power to detect changes in phasesynchronization” .
4.4. Influence of Reference
The influence of reference on EEG signals is a long-lasting problem in quantifying functional connectivity [78–80]. Various reference strategies, such as bipolar EEG , average common reference EEG , and Laplacian EEG , have been proposed. However, the cautions on these reference strategies have been extensively reported as well [81, 83].
The influence of reference has been examined with both analytical analysis and computational simulation . Let denote the reference signal, where is coefficient and is time-dependent. Let denote the scalp signal measured by the th sensor, and, denote the signal rereferenced to . For analytic signal defined with the Hilbert transform, we have Then we have as . This result indicates that the coefficient has a great influence on the MPC-based PSI . A sufficiently large will lead to larger PSI for referential signals than that for nonreferential signals . Simulation study demonstrates that the PSI of two referential signals may monotonically increase to 1 or decrease first and then increase to 1 as the coefficient increases from 0 to . In addition, a method based on independent component analysis has been demonstrated to be an appropriate method to generate reference for EEG signals in quantifying connectivity .
5. Extensions of Phase Synchronization Analysis
Various extensions of PS analysis have been proposed to infer the relationship in multivariate or multitrial signals based on the concept of IP [85–90]. For example, a method named frequency flow analysis has been used to examine the global synchronization of multivariate signals . If the IF, derived from IP, of each variable almost equals to the IFs of other variables in a certain frequency band for a duration, the multivariate signals are called in PS accordingly. In this section, we will further introduce two extensions/variations of PS analysis, that is, the trial-based PSI [85, 86] and the partial PSI .
A trial-based PSI is proposed to examine the variation of IP difference between signal channels across trials under repeated stimulus [85, 86]. For a data set of trials and each trial with samples, the trial-based PSI is defined as where and are the IPs of the channels and in the th trial. The PSI measures the intertrial variability of IP difference between channels and at instant . The PSI with value close to 1 implies that the IP difference varies little across trials at instant , while the PSI with value close to zero means that the IP difference varies approximately uniformly across trials.
The pairwise PSI discussed above can infer the strength of connectivity in two signals but cannot indicate whether the connectivity is induced by the direct coupling between them or due to the indirect interaction mediated by other units. To deal with this problem, a measure called partial PSI is generalized from PS analysis following the idea of partial coherence [87, 91]. Combined with the pairwise PSI, the partial PSI can be used to distinguish the direct and indirect interdependencies among interacted systems/units . For a set of time series , a matrix is defined with PSIs of bivariate signals as Let denote the inverse matrix of . Then a measure called partial PSI is defined as for and , conditioning on the remaining processes . A partial PSI would imply that the association is induced by indirect coupling if the pair-wise PSI between and is significant.
6. Discussions and Conclusions
We give a technical review on PS analysis in this paper. In particular, we discuss IP definitions, PSI estimation and its significance test, the issues that may affect PS detection, and extensions of PSI. PS analysis is a method to quantify the mutual rhythmic interaction of coupled systems/units and has been used to infer functional connectivity from observed neural signals such as EEG and MEG. The main advantage of PS analysis is that it could detect weak interaction between signal pairs by only taking the IPs of signals into consideration but neglecting the influence of instantaneous amplitudes of signals. In addition, PS analysis could work for nonstationary signals. These merits imply that PS analysis is suitable for neuroscience research, as we are usually interested in the relationship between neural oscillations in particular frequency bands such as beta waves ([12, 30] Hz) and gamma waves ([30, 80] Hz) rather than the interaction between broadband raw signals.
While inferring functional neural connectivity with PS analysis, several cautions and limitations should be taken into consideration. First, the observed neural signal is usually with broadband spectra and unavoidably contaminated by noise. In this case, bandpass filter should be used to extract neural oscillations in the raw signals. Second, for spatially adjacent neural recordings, PS would be affected by volume conduction and so did correlation coefficient and mutual information. In this case, the phase lag index  or the weighted phase lag index  is recommended. Third, PS analysis does not work in a black-box way. The PSIs estimated with different durations of neural signals are not recommended to be compared. Fourth, users should be aware that PSI quantifies the variation instead of the mean of IP difference within a period. Therefore, PS analysis might not be suitable for analyzing those relatively stable components in neural signals, such as event-related potentials estimated from multiple EEG trials, as PS would ignore the amplitudes of components in event-related potential and miss the latencies of these components as well.
Various measures, such as cross-correlation, coherence, nonlinear interdependence , mutual information , partial directed coherence , correlation-entropy coefficient , and coherence entropy coefficient , have been proposed in functional connectivity analysis based on EEG, MEG signals and simulated data from various aspects [10–13, 93–95]. Results show that they can reveal a similar tendency of global connectivity . The MPC-based PSI has been argued to be slightly better than coherence for both estimation and detection purposes . In another study, total 34 different measures are classified into several families such as correlation/coherence family, mutual information family, PS family, and the Granger causality family. These measures are comprehensively compared, and results show that PS measures, Granger’s causality measures and stochastic event synchrony measures are only weakly correlated with correlation coefficient, and are mutually uncorrelated as well . These results imply that these measures each could characterize interdependence of signals from different aspects. With this consideration, we could have a good characterization of functional connectivity between neural signals with only a few representative measures of different families.
This work was supported by the National Natural Science Foundation of China (nos. 60901025, 61174045, and 60935001), National Basic Research Program of China (973 Program) (no. 2011CB013304), and the International Science and Technology Cooperation Program of China (no. 2011DFA10950).
- G. Tononi, O. Sporns, and G. M. Edelman, “A measure for brain complexity: relating functional segregation and integration in the nervous system,” Proceedings of the National Academy of Sciences of the United States of America, vol. 91, no. 11, pp. 5033–5037, 1994.
- P. Tass, M. G. Rosenblum, J. Weule et al., “Detection of n:m phase locking from noisy data: application to magnetoencephalography,” Physical Review Letters, vol. 81, no. 15, pp. 3291–3294, 1998.
- F. Varela, J. P. Lachaux, E. Rodriguez, and J. Martinerie, “The brainweb: phase synchronization and large-scale integration,” Nature Reviews Neuroscience, vol. 2, no. 4, pp. 229–239, 2001.
- R. T. Canolty, E. Edwards, S. S. Dalal et al., “High gamma power is phase-locked to theta oscillations in human neocortex,” Science, vol. 313, no. 5793, pp. 1626–1628, 2006.
- J. Fell and N. Axmacher, “The role of phase synchronization in memory processes,” Nature Reviews Neuroscience, vol. 12, no. 2, pp. 105–118, 2011.
- B. He, L. Yang, C. Wilke, and H. Yuan, “Electrophysiological imaging of brain activity and connectivity—challenges and opportunities,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 7, pp. 1918–1931, 2011.
- P. J. Uhlhaas and W. Singer, “Abnormal neural oscillations and synchrony in schizophrenia,” Nature Reviews Neuroscience, vol. 11, no. 2, pp. 100–113, 2010.
- J. M. Hurtado, L. L. Rubchinsky, and K. A. Sigvardt, “Statistical method for detection of phase-locking episodes in neural oscillations,” Journal of Neurophysiology, vol. 91, no. 4, pp. 1883–1898, 2004.
- M. Ding, S. L. Bressler, W. Yang, and H. Liang, “Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment,” Biological Cybernetics, vol. 83, no. 1, pp. 35–45, 2000.
- R. Q. Quiroga, A. Kraskov, T. Kreuz, and P. Grassberger, “Performance of different synchronization measures in real data: a case study on electroencephalographic signals,” Physical Review E, vol. 65, no. 4, Article ID 041903, 14 pages, 2002.
- T. Kreuz, F. Mormann, R. G. Andrzejak, A. Kraskov, K. Lehnertz, and P. Grassberger, “Measuring synchronization in coupled model systems: a comparison of different approaches,” Physica D, vol. 225, no. 1, pp. 29–42, 2007.
- L. Astolfi, F. De Vico Fallani, F. Cincotti et al., “Estimation of effective and functional cortical connectivity from neuroelectric and hemodynamic recordings,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 17, no. 3, pp. 224–233, 2009.
- J. Dauwels, F. Vialatte, T. Musha, and A. Cichocki, “A comparative study of synchrony measures for the early diagnosis of alzheimer's disease based on EEG,” Neuroimage, vol. 49, no. 1, pp. 668–693, 2010.
- A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in Nonlinear Sciences, Cambridge University Press, Cambridge, UK, 2001.
- R. Brown and L. Kocarev, “A unifying definition of synchronization for dynamical systems,” Chaos, vol. 10, no. 2, pp. 344–349, 2000.
- M. Le Van Quyen and A. Bragin, “Analysis of dynamic brain oscillations: methodological advances,” Trends in Neurosciences, vol. 30, no. 7, pp. 365–373, 2007.
- M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “Phase synchronization of chaotic oscillators,” Physical Review Letters, vol. 76, no. 11, pp. 1804–1807, 1996.
- D. J. DeShazer, R. Breban, E. Ott, and R. Roy, “Detecting phase synchronization in a chaotic laser array,” Physical Review Letters, vol. 87, no. 4, Article ID 044101, 4 pages, 2001.
- I. Z. Kiss, Q. Lv, and J. L. Hudson, “Synchronization of non-phase-coherent chaotic electrochemical oscillations,” Physical Review E, vol. 71, no. 3, Article ID 035201, 4 pages, 2005.
- M. G. Kitzbichler, M. L. Smith, S. R. Christensen, and E. Bullmore, “Broadband criticality of human brain network synchronization,” Plos Computational Biology, vol. 5, no. 3, article e1000314, 2009.
- M. Chávez, M. Le Van Quyen, V. Navarro, M. Baulac, and J. Martinerie, “Spatio-temporal dynamics prior to neocortical seizures: amplitude versus phase couplings,” IEEE Transactions on Bio-medical Engineering, vol. 50, no. 5, pp. 571–583, 2003.
- W. Wu, J. Sun, Z. Jin et al., “Impaired neuronal synchrony after focal ischemic stroke in elderly patients,” Clinical Neurophysiology, vol. 122, no. 1, pp. 21–26, 2011.
- E. Gysels and P. Celka, “Phase synchronization for the recognition of mental tasks in a brain-computer interface,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 12, no. 4, pp. 406–415, 2004.
- C. Brunner, R. Scherer, B. Graimann, G. Supp, and G. Pfurtscheller, “Online control of a brain-computer interface using phase synchronization,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 12, pp. 2501–2506, 2006.
- P. J. Uhlhaas and W. Singer, “Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology,” Neuron, vol. 52, no. 1, pp. 155–168, 2006.
- C. Hammond, H. Bergman, and P. Brown, “Pathological synchronization in parkinson's disease: networks, models and treatments,” Trends in Neurosciences, vol. 30, no. 7, pp. 357–364, 2007.
- M. Winterhalder, B. Schelter, J. Kurths, A. Schulze-Bonhage, and J. Timmer, “Sensitivity and specificity of coherence and phase synchronization analysis,” Physics Letters Section A, vol. 356, no. 1, pp. 26–34, 2006.
- J. Sun, X. Hong, and S. Tong, “Phase synchronization analysis of EEG signals: an evaluation based on surrogate tests,” IEEE Transactions on Biomedical Engineering. Submitted.
- R. A. Dobie and M. J. Wilson, “Analysis of auditory evoked potentials by magnitude-squared coherence,” Ear and Hearing, vol. 10, no. 1, pp. 2–13, 1989.
- S. Kalitzin, J. Parra, D. N. Velis, and F. H. Lopes Da Silva, “Enhancement of phase clustering in the EEG/meg gamma frequency band anticipates transitions to paroxysmal epileptiform activity in epileptic patients with known visual sensitivity,” IEEE Transactions on Biomedical Engineering, vol. 49, no. 11, pp. 1279–1286, 2002.
- A. E. Hramov and A. A. Koronovskii, “An approach to chaotic synchronization,” Chaos, vol. 14, no. 3, pp. 603–610, 2004.
- J. Sun and M. Small, “Unified framework for detecting phase synchronization in coupled time series,” Physical Review E, vol. 80, no. 4, Article ID 046219, 2009.
- L. Xu, Z. Chen, K. Hu, H. E. Stanley, and P. C. Ivanov, “Spurious detection of phase synchronization in coupled nonlinear oscillators,” Physical Review E, vol. 73, no. 6, Article ID 065201, 2006.
- B. Schelter, M. Winterhalder, J. Timmer, and M. Peifer, “Testing for phase synchronization,” Physics Letters Section A, vol. 366, no. 4-5, pp. 382–390, 2007.
- S. R. Jammalamadaka and A. SenGupta, Topics in Circular Statistics, World Scientific Publishing, Singapore, 2001.
- E. Bedrosian, “A product theorem for Hilbert transform,” Proceedings of the IEEE, vol. 51, no. 5, pp. 868–869, 1963.
- P. J. Loughlin, “Spectrographic measurement of instantaneous frequency and the time-dependent weighted average instantaneous frequency,” Journal of the Acoustical Society of America, vol. 105, no. 1, pp. 264–274, 1999.
- P. M. Oliveira and V. Barroso, “Definitions of instantaneous frequency under physical constraints,” Journal of the Franklin Institute, vol. 337, no. 4, pp. 303–316, 2000.
- T. Pereira, M. S. Baptista, and J. Kurths, “Phase and average period of chaotic oscillators,” Physics Letters Section A, vol. 362, no. 2-3, pp. 159–165, 2007.
- G. V. Osipov, B. Hu, C. Zhou, M. V. Ivanchenko, and J. Kurths, “Three types of transitions to phase synchronization in coupled chaotic oscillators,” Physical Review Letters, vol. 91, no. 2, Article ID 024101, 4 pages, 2003.
- D. Vakman, “Computer measuring of frequency stability and the analytic signal,” IEEE Transactions on Instrumentation and Measurement, vol. 43, no. 4, pp. 668–671, 1994.
- A. H. Nuttall and E. Bedrosian, “On the quadrature approximation to the hilbert transform of modulated signals,” Proceedings of the IEEE, vol. 54, no. 10, pp. 1458–1459, 1966.
- S. Aviyente and A. Y. Mutlu, “A time-frequency-based approach to phase and phase synchrony estimation,” IEEE Transactions on Signal Processing, vol. 59, no. 7, pp. 3086–3098, 2011.
- M. Wacker and H. Witte, “Adaptive phase extraction: incorporating the Gabor transform in the matching pursuit algorithm,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 10, pp. 2844–2851, 2011.
- D. Vakman, “On the analytic signal the teager-kaiser energy algorithm and other methods for defining amplitude and frequency,” IEEE Transactions on Signal Processing, vol. 44, no. 4, pp. 791–797, 1996.
- L. Cohen, Time-Frequency Analysis, Prentice-Hall, Englewood Cliffs, NJ, USA, 1995.
- D. Hoyer, O. Hoyer, and U. Zwiener, “A new approach to uncover dynamic phase coordination and synchronization,” IEEE Transactions on Biomedical Engineering, vol. 47, no. 1, pp. 68–74, 2000.
- M. Wacker and H. Witte, “On the stability of the n:m phase synchronization index,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 2, pp. 332–338, 2011.
- M. Paluš, “Detecting phase synchronization in noisy systems,” Physics Letters Section A, vol. 235, no. 4, pp. 341–351, 1997.
- D. J. Watts and S. H. Strogatz, “Collective dynamics of “small-world” networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998.
- E. T. Bullmore and D. S. Bassett, “Brain graphs: graphical models of the human brain connectome,” Annual Review of Clinical Psychology, vol. 7, pp. 113–140, 2011.
- D. S. Bassett, E. Bullmore, B. A. Verchinski, V. S. Mattay, D. R. Weinberger, and A. Meyer-Lindenberg, “Hierarchical organization of human cortical networks in health and schizophrenia,” Journal of Neuroscience, vol. 28, no. 37, pp. 9239–9248, 2008.
- Y. Liu, M. Liang, Y. Zhou et al., “Disrupted small-world networks in schizophrenia,” Brain, vol. 131, no. 4, pp. 945–961, 2008.
- K. Wu, Y. Taki, K. Sato et al., “Age-related changes in topological organization of structural brain networks in healthy individuals,” Human Brain Mapping, vol. 33, no. 3, pp. 552–568, 2012.
- J. Sun, S. Tong, and G.-Y. Yang, “Reorganization of brain networks in aging and age-related diseases,” Aging and Disease. In press.
- C. Allefeld and J. Kurths, “Testing for phase synchronization,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 14, no. 2, pp. 405–416, 2004.
- T. Schreiber and A. Schmitz, “Surrogate time series,” Physica D, vol. 142, no. 3-4, pp. 346–382, 2000.
- J. Theiler, S. Eubank, A. Longtin, B. Galdrikian, and J. Doyne Farmer, “Testing for nonlinearity in time series: the method of surrogate data,” Physica D, vol. 58, no. 1–4, pp. 77–94, 1992.
- R. G. Andrzejak, A. Kraskov, H. Stögbauer, F. Mormann, and T. Kreuz, “Bivariate surrogate techniques: necessity, strengths, and caveats,” Physical Review E, vol. 68, no. 6, Article ID 066202, 5 pages, 2003.
- L. Faes, G. D. Pinna, A. Porta, R. Maestri, and G. Nollo, “Surrogate data analysis for assessing the significance of the coherence function,” IEEE Transactions on Biomedical Engineering, vol. 51, no. 7, pp. 1156–1166, 2004.
- M. Thiel, M. C. Romano, J. Kurths, M. Rolfs, and R. Kliegl, “Twin surrogates to test for complex synchronisation,” Europhysics Letters, vol. 75, no. 4, pp. 535–541, 2006.
- J. Kurths, M. C. Romano, M. Thiel et al., “Synchronization analysis of coupled noncoherent oscillators,” Nonlinear Dynamics, vol. 44, no. 1–4, pp. 135–149, 2006.
- J. Sun, J. Zhang, J. Zhou, X. Xu, and M. Small, “Detecting phase synchronization in noisy data from coupled chaotic oscillators,” Physical Review E, vol. 77, no. 4, Article ID 046213, 2008.
- B. C. Lovell and R. C. Williamson, “The statistical performance of some instantaneous frequency estimators,” IEEE Transactions on Signal Processing, vol. 40, no. 7, pp. 1708–1723, 1992.
- H. Hamdan, “Characterizing and approximating infinite scale mixtures of normals,” Communications in Statistics, vol. 35, no. 3, pp. 407–413, 2006.
- D. Van De Ville, J. Britz, and C. M. Michel, “EEG microstate sequences in healthy humans at rest reveal scale-free dynamics,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 42, pp. 18179–18184, 2010.
- P. J. Uhlhaas, D. E. J. Linden, W. Singer et al., “Dysfunctional long-range coordination of neural activity during gestalt perception in schizophrenia,” Journal of Neuroscience, vol. 26, no. 31, pp. 8168–8175, 2006.
- S. Hu, M. Stead, Q. Dai, and G. A. Worrell, “On the recording reference contribution to EEG correlation, phase synchorony, and coherence,” IEEE Transactions on Systems, Man, and Cybernetics, Part B, vol. 40, no. 5, pp. 1294–1304, 2010.
- G. V. Tcheslavski, Coherence and phase synchrony analysis of electroencephalogram, Ph.D. dissertation, the Virginia Polytechnic Institute, State University, Blacksburg, Va, USA, 2005.
- M. Vinck, R. Oostenveld, M. Van Wingerden, F. Battaglia, and C. M. A. Pennartz, “An improved index of phase-synchronization for electrophysiological data in the presence of volume-conduction, noise and sample-size bias,” Neuroimage, vol. 55, no. 4, pp. 1548–1565, 2011.
- M. Vinck, M. van Wingerden, T. Womelsdorf, P. Fries, and C. M. A. Pennartz, “The pairwise phase consistency: a bias-free measure of rhythmic neuronal synchronization,” Neuroimage, vol. 51, no. 1, pp. 112–122, 2010.
- J. Sarvas, “Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem,” Physics in Medicine and Biology, vol. 32, no. 1, article 004, pp. 11–22, 1987.
- P. L. Nunez, R. Srinivasan, A. F. Westdorp et al., “EEG coherency i: statistics, reference electrode, volume conduction, laplacians, cortical imaging, and interpretation at multiple scales,” Electroencephalography and Clinical Neurophysiology, vol. 103, no. 5, pp. 499–515, 1997.
- J. G. Stinstra and M. J. Peters, “The volume conductor may act as a temporal filter on the ecg and EEG,” Medical and Biological Engineering and Computing, vol. 36, no. 6, pp. 711–716, 1998.
- G. Nolte, O. Bai, L. Wheaton, Z. Mari, S. Vorbach, and M. Hallett, “Identifying true brain interaction from EEG data using the imaginary part of coherency,” Clinical Neurophysiology, vol. 115, no. 10, pp. 2292–2307, 2004.
- C. J. Stam, G. Nolte, and A. Daffertshofer, “Phase lag index: assessment of functional connectivity from multi channel EEG and meg with diminished bias from common sources,” Human Brain Mapping, vol. 28, no. 11, pp. 1178–1193, 2007.
- L. A. Wheaton, G. Nolte, S. Bohlhalter, E. Fridman, and M. Hallett, “Synchronization of parietal and premotor areas during preparation and execution of praxis hand movements,” Clinical Neurophysiology, vol. 116, no. 6, pp. 1382–1390, 2005.
- L. T. Trujillo, M. A. Peterson, A. W. Kaszniak, and J. J. B. Allen, “EEG phase synchrony differences across visual perception conditions may depend on recording and analysis methods,” Clinical Neurophysiology, vol. 116, no. 1, pp. 172–189, 2005.
- S. J. Schiff, “Dangerous phase,” Neuroinformatics, vol. 3, no. 4, pp. 315–317, 2005.
- R. Guevara, J. L. P. Velazquez, V. Nenadovic, R. Wennberg, G. Senjanović, and L. G. Dominguez, “Phase synchronization measurements using electroencephalographic recordings: what can we really say about neuronal synchrony?” Neuroinformatics, vol. 3, no. 4, pp. 301–313, 2005.
- H. P. Zaveri, R. B. Duckrow, and S. S. Spencer, “On the use of bipolar montages for time-series analysis of intracranial electroencephalograms,” Clinical Neurophysiology, vol. 117, no. 9, pp. 2102–2108, 2006.
- O. Bertrand, F. Perrin, and J. Pernier, “A theoretical justification of the average reference in topographic evoked potential studies,” Electroencephalography and Clinical Neurophysiology, vol. 62, no. 6, pp. 462–464, 1985.
- H. P. Zaveri, R. B. Duckrow, and S. S. Spencer, “The effect of a scalp reference signal on coherence measurements of intracranial electroencephalograms,” Clinical Neurophysiology, vol. 111, no. 7, pp. 1293–1299, 2000.
- S. Hu, M. Stead, and G. A. Worrell, “Automatic identification and removal of scalp reference signal for intracranial EEGs based on independent component analysis,” IEEE Transactions on Biomedical Engineering, vol. 54, no. 9, pp. 1560–1572, 2007.
- J.-P. Lachaux, E. Rodriguez, J. Martinerie, and F. J. Varela, “Measuring phase synchrony in brain signals,” Human Brain Mapping, vol. 8, no. 4, pp. 194–208, 1999.
- L. Li, J. X. Zhang, and T. Jiang, “Visual working memory load-related changes in neural activity and functional connectivity,” Plos One, vol. 6, no. 7, Article ID e22357, 2011.
- B. Schelter, M. Winterhalder, R. Dahlhaus, J. Kurths, and J. Timmer, “Partial phase synchronization for multivariate synchronizing systems,” Physical Review Letters, vol. 96, no. 20, Article ID 208103, 2006.
- D. Rudrauf, A. Douiri, C. Kovach et al., “Frequency flows and the time-frequency dynamics of multivariate phase synchronization in brain signals,” Neuroimage, vol. 31, no. 1, pp. 209–227, 2006.
- H. Osterhage, F. Mormann, T. Wagner, and K. Lehnertz, “Detecting directional coupling in the human epileptic brain: limitations and potential pitfalls,” Physical Review E, vol. 77, no. 1, Article ID 011914, 2008.
- C. Allefeld and S. Bialonski, “Detecting synchronization clusters in multivariate time series via coarse-graining of markov chains,” Physical Review E, vol. 76, no. 6, Article ID 066207, 2007.
- R. Dahlhaus, “Graphical interaction models for multivariate time series,” Metrika, vol. 51, no. 2, pp. 157–172, 2000.
- A. Gunduz and J. C. Principe, “Correntropy as a novel measure for nonlinearity tests,” Signal Processing, vol. 89, no. 1, pp. 14–23, 2009.
- A. V. Sazonov, C. K. Ho, J. W. M. Bergmans et al., “An investigation of the phase locking index for measuring of interdependency of cortical source signals recorded in the EEG,” Biological Cybernetics, vol. 100, no. 2, pp. 129–146, 2009.
- A. Bruns, “Fourier-, Hilbert- and wavelet-based signal analysis: Are they really different approaches?” Journal of Neuroscience Methods, vol. 137, no. 2, pp. 321–332, 2004.
- O. David, D. Cosmelli, and K. J. Friston, “Evaluation of different measures of functional connectivity using a neural mass model,” Neuroimage, vol. 21, no. 2, pp. 659–673, 2004.