Differentiating resting brain states using ordinal symbolic analysis

Differentiating resting brain states using ordinal symbolic analysis

C. Quintero-Quiroz Universitat Politècnica de Catalunya, Departament de Física, Colom 11, 08222 Terrassa, Barcelona, Spain.    Luis Montesano Bitbrain, Zaragoza, Spain.    A. J. Pons    M. C. Torrent Universitat Politècnica de Catalunya, Departament de Física, Colom 11, 08222 Terrassa, Barcelona, Spain.    J. García-Ojalvo Department of Experimental and Health Sciences, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain.    C. Masoller Universitat Politècnica de Catalunya, Departament de Física, Colom 11, 08222 Terrassa, Barcelona, Spain.
July 20, 2019
Abstract

Symbolic methods of analysis are valuable tools for investigating complex time-dependent signals. In particular, the ordinal method defines sequences of symbols according to the ordering in which values appear in a time series. This method has been shown to yield useful information, even when applied to signals with large noise contamination. Here we use ordinal analysis to investigate the transition between eyes closed (EC) and eyes open (EO) resting states. We analyze two EEG datasets (with 71 and 109 healthy subjects) with different recording conditions (sampling rates and number of electrodes in the scalp). Using as diagnostic tools the permutation entropy, the entropy computed from symbolic transition probabilities, and an asymmetry coefficient (that measures the asymmetry of the likelihood of the transitions between symbols) we show that ordinal analysis applied to the raw data distinguishes the two brain states. In both datasets we find that the EO state is characterized by higher entropies and lower asymmetry coefficient, as compared to the EC state. Our results thus show that these diagnostic tools have potential for detecting and characterizing changes in time-evolving brain states.

time series analysis, ordinal analysis, EEG, brain dynamics

In the “big data” era, many efforts are being devoted to extracting useful information from complex signals. The human brain is one of the most complex systems that one can try to understand. In the last decades, the development and popularization of recording techniques such as electroencephalography (EEG), magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI), have provided the scientific community with a huge amount of data: different types of brain signals, recorded with different spatio-temporal resolution, under different behavioral or cognitive states, from healthy or from dysfunctioning subjects. The underlying brain states are, in spite of many efforts, still poorly understood. Here we use a symbolic analysis tool to investigate EEG signals recorded from healthy subjects during a simple behavioral task: the subjects remain in resting state with eyes closed (EC state) during an interval of time, and then open their eyes (EO state). We show that symbolic analysis applied to the raw EEG signals detects the transition and identifies subtle differences between the EC and EO brain states.

I Introduction

Changes in brain states detected through the analysis of electroencephalography (EEG) signals can be used for translating brain signals into operational commands, and in fact, EEG analysis is one of the techniques used for brain-computer interfaces.

Several methods have been used to detect underlying changes in the behavior of dynamical systems from observed data, and one of these, ordinal analysisBandt and Pompe (2002); Cao et al. (2004); Masoller et al. (2015), has been demonstrated to be computationally efficient and to perform well even with very noisy data Zanin et al. (2012); Quintero-Quiroz et al. (2015). Due to these advantages, ordinal analysis has been used in the field of neuroscience, specifically in the area of epilepsy, for detecting, anticipating and characterizing seizures Li et al. (2014); Nicolaou and Georgiou (2012); Olofsen, Sleigh, and Dahan (2008); Bruzzo et al. (2008); Li, Ouyang, and Richards (2007); Veisi, Pariz, and Karimpour (2007); Ren et al. (2015).

Since the early 1930’s it is well known that alpha waves dominate the EEGs of healthy individuals when they are resting with their eyes closed, and that this activity diminishes when their eyes are opened Smith (1938); Jasper (1936); Adrian and Matthews (1934); Berger (1933). Therefore, a simple method to detect the Eyes-Closed (EC) to Eyes-Open (EO) transition is by using the Fourier spectrum to estimate the difference of the power of the alpha frequency components Barry et al. (2007); Barry and De Blasio (2017). However, this approach has the drawback of requiring a certain time-window for computing the power spectra. Another approach to study the EC-EO transition is to use the synchronization likelihoodJin et al. (2014) or the mutual informationTan et al. (2013) to find changes in the functional brain networks that characterize the two brain states. However, constructing functional brain networks is computationally demanding, and comparing them is a challenging task because it is not always possible to discriminate reliably between differences that are due to constrains imposed by method of network construction, or due to genuine changes in brain statesVan Wijk, Stam, and Daffertshofer (2010); Schieber et al. (2017).

The aim of this paper is to investigate if the ordinal approach can accurately discriminate between EC and EO brain states. In Sec. II we describe the datasets analyzed, in Sec. III we describe the ordinal-pattern methodology and the quantifiers used to characterize the EC and EO states. Sec. IV presents the results obtained and Sec. V summarizes our conclusions.

Ii Datasets

We use two EEG datasets with different number of subjects and recording conditions, which are summarized in the Table 1. Dataset one (DTS1) was collected by the Bitbrain companyBitbrain Technologies (). The EEG signals were recorded from 71 healthy subjects that remained with eyes closed and eyes open during a period of two minutes each. Dataset two (DTS2), which is freely availableSchalk et al. (2004); Goldberger et al. (2000), consists of EEG recordings of 109 subjects performing the same task, in this case for a period of one minute in each of the two states.

DTS1 DTS2
Sampling rate(Hz) 256 160
Time task(seg) 120 60
Total points 30720 9600
Number of electrodes 16 64
Number of subjects 71 109
Table 1: Description of the datasets used.

We removed the artifacts related to eye blinking following the standard procedure: we applied the Independent Component Analysis (ICA) using the function ICA from the MNE library on PythonGramfort et al. (2014, 2013) and filtered out the component related to the blinks (see Fig. 1).

Figure 1: Example of Independent Component Analysis (ICA) related to eye blinks in the EEGs, for a given subject of DTS1. Left: spatial contribution of the selected ICA component, ploted in red in the left panel. Right: all the individual components obtained from the ICA function, in red the component related to the eye blinks, shown in the left panel.

It is well known that alpha waves are a dominant component in EEG signals during eyes closed conditions, and are reduced when the eyes are openBerger (1933); Barry et al. (2007). Therefore, in order to determine whether changes detected through ordinal analysis are only due to the change of the strength of alpha waves, we analyze and compare the results obtained from the raw time series, and from filtered time series where both eye blinking artifacts and the frequency component of the alpha band were removed (by using a band pass filter between 14 and 31 Hz, see Fig. 2).

Figure 2: Filtering of the alpha band from the ICA data. The left panel displays the power spectrum of the EEG of a subject before (blue line) and after (green line) filtering. The right top panel displays the post-processed EEG’s time-series (after the filtering), the vertical line indicates the time of the eyes closed – eyes open transition. The subject is the same as in Fig. 1.

Iii Methods

We apply ordinal analysis in non-overlapping windows of 1 secondNicolaou and Georgiou (2012), and thus, the number of data points in the window, , is equal to for DTS1 and to for DTS2. Then, for each electrode , the time-series, is transformed into a sequence of symbols, , by using the ordinal ruleBandt and Pompe (2002); Zanin et al. (2012), explained in what follows.

To define the ordinal patterns we consider vectors of dimension formed by consecutive data points, i.e. , and then assign a symbol according to the ordinal relationship (from the largest to the smallest value) of the entries in the vectors. For example, with there are 2 ordinal patterns (): corresponding to the ordinal pattern ‘01’ and corresponding to the ordinal pattern ‘10’. Then we computed the frequency of occurrence of the D! different patterns in the signal of electrode , and by averaging over all the electrodes, computed the probability of each pattern.

Then, the permutation entropy (PE) is calculated as:

(1)

where is the probability of pattern along all the electrodes. In this way, the PE is a measure of the entropy of the brain EEG signals, in the given time window. If the EEG signals are generated by fully random processes, all symbols are equally probable and the PE is maximum, .

Additional diagnostic tools were proposed by Masoller et al.Masoller et al. (2015), which are based in the transition probabilities (TPs) between consecutive symbols defined from non-overlapping data values. The transition probability from pattern to pattern is the relative number of times pattern is followed by pattern , in the sequence :

(2)

With this definition, the transition probabilities are normalized such that . Then, exploiting this normalization, an entropy can be associated to the transition probabilities of each pattern as , and its average

(3)

is another measure of the entropy of the EEG signal. If a signal is generated by fully random processes, all transition probabilities will be equal and thus, for all , and . In addition, we calculate the transition asymmetry coefficientMasoller et al. (2015),

(4)

which is equal to zero if transition probabilities are fully symmetric ( for all ), and equal to one if they are fully asymmetric (either or , for all ). If the EEG signals are generated by fully random processes, then the transition probabilities will be all equal and .

In the following section the analysis is performed with non overlapping patterns of length (similar results were found with ). There are possible patterns and 2424=576 possible transitions. For the dataset DTS1 (DTS2), taking together the 16 (71) electrodes, in each time window the symbolic sequence contains 4048 (10048) patterns and 4032 (9984) transitions. While the number of patterns is clearly sufficient to compute the probabilities of the 24 patterns with good accuracy, longer sequences are needed to compute the 576 TPs with similar accuracy. Nevertheless, we will show that the TP-based diagnostic tools, and , can also detect changes in datasets. As a measure of statistical significance we calculate the -value using Welch’s t-test and consider, as null hypothesis, that the signals represent the same state.

Iv Results

We begin by calculating the PE for the raw, unfiltered time series. Figure 3 displays the results obtained from DTS1 and DTS2 and we can see that there is a significant difference between the PE values of the eyes-closed and eyes-open states. The entropy is computed for each subject, and then is averaged over all the subjects (71 or 109, depending on the dataset). The shaded area represents one standard deviation of the PE values of all subjects, and we note that there is large variability, thus, the PE value does not allow a full discrimination between the two states. We note that the average PE value is slightly different for the two datasets, which is attributed to the fact that they have different spatial and temporal resolution. We also note that the average value of the PE is significantly different from the maximum possible value (which occurs when the patterns are equally probable, and for , PE). This reveals the presence of patterns with high and low frequency of occurrence in the symbolic sequence. An inspection of the ordinal probabilities reveals that the “trend” patterns (generated by consecutively increasing or decreasing data values) are more expressed, in all the channels and for both, EO and EC states. For future work, it would be interesting to investigate if the probabilities of certain patterns (which could be defined by using a lag, as inParlitz et al. (2012)), allow for a better discrimination of the two states.

Figure 3: Permutation entropy, Eq. 1, from raw time series of DTS1 (top) and DTS2 (bottom). In DTS1 the subjects open their eyes at 120 s; in DTS2, at 60 s. The blue line indicates the mean value of the PE for all the subjects, and the shaded area indicates one standard deviation of the PE values.

Comparing with the results obtained from the filtered time series, displayed in Fig. 4, we note that the PE values remain almost unchanged, which suggests that the PE captures changes in brain dynamics which are not due to the change in the strength of alpha oscillations during the EC-EO transition.

Figure 4: Permutation entropy, Eq. 1, computed from filtered time series of DTS1 (top) and DTS2 (bottom).

Figure 5 displays the TP-based measures, and , using the DST2 (similar results were found in DST1), although there is a clear transition around 60s in the mean values, the dispersion in the values of the different subjects is higher than in the PE analysis (likely due to the limited length of the time series, which does not allow a precise estimation of the TPs).

Figure 5: The entropy defined from the transition probabilities, Eq. 3, and the asymmetry coefficient, Eq. 4, computed for the DTS2.

In Fig. 6 (for DTS1) and Fig. 7 (for DTS2) we present the topographic visualization for the different electrodes, of the PE value averaged over all the subjects, for the EC and EO conditions. We also present the difference of PE values (PE-open - PE-closed), and the -value. The results are consistent for the two datasets, the discrepancies are due to their different spatial resolution. We also note the low p-values obtained, which confirm the significance of the uncovered differences.

Figure 6: Topographic visualization of the analysis of the raw (top row) and filtered (lower row) EEG signals of DTS1, average over the subjects. Panels a) and e) display the permutation entropy for EC conditions; b) and f) for EO conditions; c) and g) display the difference of the PE values; d) and h) display the value.
Figure 7: As in Fig. 6, but for the dataset DTS2.

V Conclusions

We have used ordinal time series analysis to investigate EEG signals recorded under eyes-closed (EC) and eyes-open (EO) resting conditions. We have analyzed two datasets with different spatial and temporal resolutions, and contrasted the results of the analysis of raw time series and filtered time series (where eye blinking artifacts were removed and the alpha frequency band was filtered out). We used three diagnostic measures, the permutation entropy, PE, which is computed from the probabilities of the ordinal patterns, and two measures, the transition entropy and the asymmetry coefficient, which are computed from the transition probabilities between patterns.

We have found, in both datasets, that the EO state is characterized by higher entropy values, accompanied by a lower asymmetry coefficient, with respect to the EC state. We have also identified which brain regions are more important for distinguishing the two states. No significant difference was detected between the raw data and the pre-processed data, which suggests that the ordinal method can be directly applied to EEG signals, avoiding the need of data pre-processing. Thus, ordinal analysis can be a computationally efficient tool, which could provide extra valuable information for new brain-computer interface protocols.

Acknowledgments

This work was supported in part by ITN NETT (FP7 289146), the Spanish MINECO (FIS2015-66503 and FIS2015-66503-C3-2-P) and the program ICREA ACADEMIA of Generalitat de Catalunya.

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
338823
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description