Detecting coherence between oscillators in heavy-noise environments with a moving block bootstrap

Kazimieras Pukenas1

1Lithuanian Sports University, Kaunas, Lithuania

1Corresponding author

Mathematical Models in Engineering, Vol. 1, Issue 1, 2015, p. 20-25.
Received 25 February 2015; received in revised form 20 April 2015; accepted 10 May 2015; published 30 June 2015

Copyright © 2015 JVE International Ltd.

Table of Contents Abstract Full-text  Download PDF References
Views 0
Reads 0
Downloads 21
since May, 2018
Abstract.

In this paper, a novel algorithm based on the moving block-bootstrap (MBB) technique for the detection of coherence between oscillators in heavy-noise environments is proposed. To verify the null hypothesis of absent significant coherence against the alternative, the average coherence of MBB resampled data is compared with a certain threshold over all frequencies. The threshold is defined as the upper bound of a 95 % confidence interval for bootstrapping coherence performed with independent resampling indexes for both time series, so the dependence between the time series under consideration is destroyed. The benefits of the proposed method are illustrated by simulations of the phase synchronization effects of two linearly coupled chaotic Rössler oscillators.

Keywords: coherence, moving block bootstrap, null hypothesis, confidence interval.

1. Introduction

Identifying the couplings between processes is an important task in many areas of science and engineering. These interdependences have been quantified in various ways. Among the most robust and therefore the most frequently applied techniques are synchronization measures and cross-spectral analysis. This resulted from the fact that both approaches are robust with respect to noise, one of the key requirements in most applications. The cross-spectral analysis approach additionally provides information on the frequency at which an interaction is strongest, a feature that is strongly desirable in applications to noisy nonlinear oscillatory systems [1]. Moreover, it is argued [2, 3], that for noisy data in some cases the results for Fourier coherence are (slightly) more robust than for phase synchrony; if the signal is weak it is more likely that noise will destroy the phase structure. Detecting relationships between experimental data becomes a challenge when the data are very noisy. When the signals are nonlinear, and especially chaotic, applying the noise reduction techniques, including chaos-based approaches for such a purpose, can lead to erroneous estimation of the actual degree of interdependence between oscillators with narrowing bandwidth [4, 5]. For high noise levels significant coherence detection can be formulated as hypothesis testing. Moving block bootstrap (MBB) method [6-10] was introduced as an adaption of the ordinary bootstrap for serially correlated data and significantly improves statistical inference for noisy data. The general idea of MBB is to provide the distribution of a statistics based on the measured signals alone, when the analytic derivation of the statistics is not known. Once this empirical distribution is obtained from the bootstrap, confidence intervals can be derived and hypothesis tests can be performed based on the empirical quantiles. Based on this technique, the significance of the coupling between two time series is commonly assessed by setting a threshold level in the coherence function [6-8]. When noise is strong, however, the accuracy of the MBB approach to calculate confidence intervals from the asymptotic distribution function of the sample coherence function is degraded. On the other hand, it has been shown that some approaches based on surrogate data and the corresponding null hypothesis that the oscillators are unsynchronized or expected coherence is equal to zero at all frequencies can lead to false conclusions, because surrogate techniques like phase randomization suffer from the problem that not only is the interaction between the systems neutralized but the systems themselves are also linearized [11, 12]. This problem calls for a threshold level to be set in the estimated genuine coherence in order to assess whether the two series are significantly coupled or not.

The aim of this paper is to propose a simple coherence detection algorithm in heavy-noise environments based on the two-way MBB technique with a novel defined significance threshold. The performances of the proposed method are illustrated by simulations of the phase synchronization effects between two coupled chaotic Rössler oscillators.

2. Coherence and moving block bootstrap MBB

The coherence between two processes is a measure of the linear relationship of the two at a specific frequency. For two real stationary zero-mean processes Xt and Yt, the cross spectrum SXYω is defined as the Fourier transform of the cross-covariance function FXYτ=XtYt-τ [13]. The coherence function, also called the magnitude squared coherence function, CXYω2 at frequency ω is given by the modulus of the cross-spectrum SXYω normalized by the respective auto-spectra SXω and SYω [13]:

(1)
C X Y ω 2 = S X Y ω 2 S X ω S Y ω   ,         - < ω < .

For simplicity we shall refer to CXYω2 as the coherence. It measures the extent to which Yt is determinable from Xt at frequency ω by linear time-invariant operations and is bounded between zero and one. An estimate of the coherence C^XYω2 is obtained by replacing the cross- and auto-spectra in Eq. (1) with their respective estimated quantities obtained by averaging or smoothing cross- and auto-periodograms.

Briefly, MBB is a nonparametric bootstrap method that can be applied to dependent observations [6-9]. The concept of MBB consists of calculating the estimator of interest over replicated pseudo-data obtained by random drawing with replacement from the blocks of consecutive data [8]. In more detail, from the n available samples of the time series Xt random selection of blocks of length b from the original data takes place. Then these blocks are concatenated together to obtain a replicated bootstrap pseudo-data. As a result of the MBB procedure, approximations of unknown distributions of root statistics are obtained. These approximations allow the construction of confidence intervals, and then building tests for different problems. Typically, hundreds and thousands of resamples are recommended for reliable testing.

3. Description of the algorithm

To obtain a mean coherence measure with MBB, we proceed as follows. Let X=X1, Y1, , Xn,Yn be a bivariate time series of length n. Let the integer b be the block length, define Bt,b=Xt,Yt,, Xt+b-1,Yt+b-1, and k=n/b is assumed to be an integer. Using random number generator to generate k integer values i1,…, ik from the uniform distribution on the set 1, 2, …, n-b+1 and joining the k blocks each of length b, such as Bi1,b, …, Bik,b, we draw the MBB sample X*=X1*, Y1*, , Xn*,Yn* of the same size n. The MBB version of the coherence estimator C^XY*ω2 is calculated by replacing in Eq. (1) the estimates S^Xω, S^Yω, and S^XYω with their bootstrap counterparts, i. e.:

(2)
C ^ X Y * ω 2 = S ^ X Y * ω 2 S ^ X * ω S ^ Y * ω   .

Repeating the MBB procedure B times we obtain a total of N bootstrap statistics C^XY*1ω2,…, C^XY*Nω2 and a mean coherence measure is calculated:

(3)
C ^ X Y * ω 2 = 1 B i = 1 B C ^ X Y * i ω 2 .

The mean coherence measure radically reduces spurious maxima that may lead to misleading interpretations. However, for high noisy time series, the variance at some frequencies outside the dominant frequency range:

(4)
σ 2 * ω = B - 1 i = 1 B C ^ X Y * i ω 2 - B - 1 j = 1 B C ^ X Y * j 2 2 ,

remains high, and non-zero mean coherence values are likely to occur at some frequencies outside the dominant frequency range, as illustrated in Fig. 1.

Moreover, the upper bound of the 95 % confidence interval of spurious peaks can exceed the estimated mean coherence at dominant frequency. This problem calls for appropriate threshold level to be set in the estimated coherence so we can assess whether the two series are significantly coupled or not. We propose a novel but simple method for defining the threshold, which is based on two MBB procedures: a) bootstrap resampling in the context of bivariate time series model and b) bootstrap resampling for each variable separately. Namely, the threshold is defined as the upper bound of a 95 % confidence interval for bootstrapping coherence performed with independent bootstrap indexes for time series Xt and Yt, which results in destroying of the dependence between the time series under consideration. The bootstrap using independent bootstrap indexes for time series Xt and Yt can be treated as a two-dimensional surrogate method in respect of coherence. As a result, the structural properties of the original time series are maintained but the bootstrapped series are structurally uncoupled. Because no dominant coherence between randomly resampled time series Xt and Yt should occur, it is reasonable to define the constant threshold as the new upper bound of a 95 % confidence interval over all frequencies of the spectrum. The null hypothesis that the two series are not coherent against the alternative that the two series are significantly coupled, i. e., the observed mean coherence measure derived from the bootstrap statistics is greater than the defined threshold, is examined by performing t-tests at each frequency of the spectrum.

Fig. 1. The mean coherence measure with 95 % confidence bounds between two coupled Rössler oscillators derived after 480 bootstrap realizations at SNR = –15 dB

4. Simulation results

We consider the following system of two coupled chaotic Rössler oscillators [14, 15]:

(5)
x ˙ 1,2 = - ω 1,2 y 1,2 - z 1,2 + c x 2,1 - x 1,2 ,
y ˙ 1,2 = ω 1,2 x 1,2 + 0.15 y 1,2 ,
z ˙ 1,2 = 0.2 + z 1,2 x 1,2 - 8.5 ,

where ω1= 0.98 and ω2= 1.02. With coupling strength c= 0.035, the coherent systems are synchronized. The spectral peak of the measured data is located around the normalized frequency 0.05 (×π rad/sample). We generate a time series using the ODE45 integrator of Matlab. The solution was sampled at time intervals t= 0.1 and 2500 samples, after the first 200 samples have been discarded, are adopted for analysis. The magnitude squared coherence of the input signals Xt and Yt were estimated by Welch’s averaged, modified periodogram method and periodic Hamming windows of length to obtain eight equal sections Xt and Yt with a 50 % overlap were used. All data processing and analyses were performed with Matlab software (The MathWorks, Natick, MA).

For moving block-bootstrapping, the appropriate choice of the block length b is a key parameter, and the choice of b for a small sample size is not straightforward. For the bootstrap to be effective, b should be large enough for as much of the dependence structure as possible to be retained in the blocks, but not so large that the number of blocks becomes small, resulting in poor estimation of underlying distribution [6]. The method based on the decay rate of the autocorrelation function [16] and a similar approach with respect to discontinuity of the autocorrelation function at time lag zero for noisy time series [7] were proposed for selection of the optimal block length. At strength noise, however, the decay rate of the autocorrelation function cannot be reliably estimated. Therefore, block length, estimated from the decay rate of the autocorrelation function (by fitting an exponential function to the envelope of the empirical autocorrelation function and excluding the zero time lag) was approximated to the floor integer of periods (six) of the autocorrelation function. A similar solution is used in seasonal block bootstrap, which was designed for time series with periodic structure and in which the block length is an integer multiple of the period length [17].

Fig. 2. The mean coherence values (black) between two coupled Rössler oscillators and threshold values (red) derived from 100 independent Monte Carlo trials at SNR= – 15 dB

The performance of the algorithm was evaluated by means of a 100-trial Monte Carlo simulation whereby the signals were contaminated by various realizations of additive zero mean white Gaussian noise at a defined signal-to-noise ratio (SNR). Furthermore, these simulations were performed for two versions of random selection of blocks: (i) the blocks of length b started at the time indexes i1,…, ik, where i1,…, ik, are iid from a discrete uniform distribution 1, 2, …, n-b+1 as shown above, and (ii) the blocks were selected randomly but with respect to the pseudo periodicity of the processes under consideration, i.e., the blocks started at the time indexes, the distance between which is an integer multiple of the approximated period length. Number of bootstrap realizations B= 480 was chosen at every trial. Fig. 2 shows the mean coherence values from 100 independent Monte Carlo trials as well as defined threshold values at SNR = –15 dB. As can be seen, the algorithm reliably detects the coherence peak at the oscillation frequency in the vast majority of trials and a small number of false positive errors can occur at some frequencies. The false positive and false negative error rates across trials at different SNR are shown in Table 1.

Table 1. False positive and false negative error rates at different SNR

Start points of blocks
Error rate, %
SNR
–5 dB
–10 dB
–15dB
Full random
False positive
2
2
1
False negative
0
0
2
With respect to the periodicity of the processes
False positive
2
1
1
False negative
0
0
0

5. Conclusions

All simulations results performed on coupled chaotic Rössler oscillators, verify the applicable of the proposed algorithm for detecting coherence at the oscillation frequency between oscillators in a strong white Gaussian noise environment. In accordance with this approach, the null hypothesis of absence significant coherence (i.e., at all frequencies) is rejected if the average coherence of MBB resampled data exceeds the 1001α % quantile for bootstrapping coherence obtained from the surrogate time series, i. e., derived from the MBB procedure using independent resampling indexes for the time series Xt and Yt. Differently from other surrogate techniques like phase randomization, when not only is the interaction between the systems neutralized but the systems themselves are also linearized [12], the MBB procedure using independent resampling indexes for both time series preserves the structural properties of the original time series. The simulations show that the sensitivity of the algorithm increases (without loss of the specificity) if the distance between blocks is defined in terms of the oscillation period. It is worth noting that for practical applications, significant coherence between oscillators under strong white Gaussian noise background can be detected subjectively on the basis of the fairly pronounced spike at the oscillation frequency of the lower 5 % confidence interval of bootstrapped coherence.

References

  1. Sommerlade L., Mader M., Mader W., Timmer J., Thiel M., Grebogi C., Schelter B. Optimized spectral estimation for nonlinear synchronizing systems. Physical Review E, Vol. 89, 2014, p. 032912.
  2. Nolte G., Bai O., Wheaton L., Mari Z., Vorbach Sh., Hal M. Identifying true brain interaction from EEG data using the imaginary part of coherency. Clinical Neurophysiology – Journal, Vol. 115, 2004, p. 2292-2307.
  3. Zhang H., Benz H. L., Thakor N. V., Bezerianos A. Connectivity mapping of human brain by phase based evolution map approach. International Journal of Bifurcation and Chaos, Vol. 22, Issue 9, 2012, p. 1250225.
  4. Tung W., Gao J., Hu J., Yang L. Detecting chaos in heavy-noise environments. Physical Review E, Vol. 83, 2011, p. 046210.
  5. Xu L., Chen Z., Hu K., Stanley H. E., Ivanov P. Ch. Spurious detection of phase synchronization in coupled nonlinear oscillators. Physical Review E, Vol. 73, 2006, p. 065201.
  6. Zoubir A., Iskander R. Bootstrap Techniques for Signal Processing. Cambridge University Press, New York, 2004.
  7. Mader M., Mader W, Sommerlade L., Timmer J., Schelter B. Block-bootstrapping for noisy data. Journal of Neuroscience Methods, Vol. 219, 2013, p. 285-291.
  8. Maiz S., Elbadaoui M., Bonnardot F., Serviere Ch. New second order cyclostationary analysis and application to the detection and characterization of a runner's fatigue. Signal Processing, Vol. 102, 2014, p. 188-200.
  9. Synowiecki R. Consistency and application of moving block bootstrap for nonstationary time series with periodic and almost periodic structure. Bernoulli, Vol. 13, Issue 4, 2007, p. 1151-1178.
  10. Allefeld C., Kurths J. Testing for phase synchronization. International Journal of Bifurcation and Chaos, Vol. 14, Issue 2, 2004, p. 405-416.
  11. Faes L., Pinna G. D., Porta A., Maestri R., Nollo G. Surrogate data analysis for assessing the significance of the coherence function. IEEE Transactions on Biomedical Engineering, Vol. 51, Issue 7, 2004, p. 1156-1166.
  12. Schelter B., Winterhalder M., Timmer J., Peifer M. Testing for phase synchronization. Physics Letters A, Vol. 366, 2007, p. 382-390.
  13. Bloomfeld P. Fourier Analysis of Time Series: an Introduction. Wiley, New York, 1976.
  14. Lai Y.-Ch., Frei M. G., Osorio I. Detecting and characterizing phase synchronization in nonstationary dynamical systems. Physical Review E, Vol. 73, 2006, p. 026214.
  15. Sun J., Small M. Unified framework for detecting phase synchronization in coupled time series. Physical Review E, Vol. 80, 2009, p. 046219.
  16. Peifer M., Schelter B., Guschlbauer B., Hellwig B., Lücking C., et al. On studentising and block length selection for the bootstrap on time series. Biometrical Journal, Vol. 47, Issue 3, 2005, p. 346-357.
  17. Dudek A. E., Maiz S., Elbadaoui M. Generalized seasonal block bootstrap in frequency analysis of cyclostationary signals. Signal Processing, Vol. 104, 2014, p. 358-368.