Measure of the electroencephalographic effects of sevoflurane using recurrence dynamics

Network Analysis of Time Series: Novel Approaches to Network Neuroscience

In the last two decades, there has been an explosion of interest in modeling the brain as a network, where nodes correspond variously to brain regions or neurons, and edges correspond to structural or statistical dependencies between them. This kind of network construction, which preserves spatial, or structural, information while collapsing across time, has become broadly known as "network neuroscience." In this work, we provide an alternative application of network science to neural data: network-based analysis of non-linear time series and review applications of these methods to neural data. Instead of preserving spatial information and collapsing across time, network analysis of time series does the reverse: it collapses spatial information, instead preserving temporally extended dynamics, typically corresponding to evolution through some kind of phase/state-space. This allows researchers to infer a, possibly low-dimensional, "intrinsic manifold" from empirical brain data. We will discuss three methods of constructing networks from nonlinear time series, and how to interpret them in the context of neural data: recurrence networks, visibility networks, and ordinal partition networks. By capturing typically continuous, non-linear dynamics in the form of discrete networks, we show how techniques from network science, non-linear dynamics, and information theory can extract meaningful information distinct from what is normally accessible in standard network neuroscience approaches.

Information-Theoretical Quantifier of Brain Rhythm Based on Data-Driven Multiscale Representation

This paper presents a data-driven multiscale entropy measure to reveal the scale dependent information quantity of electroencephalogram (EEG) recordings. This work is motivated by the previous observations on the nonlinear and nonstationary nature of EEG over multiple time scales. Here, a new framework of entropy measures considering changing dynamics over multiple oscillatory scales is presented. First, to deal with nonstationarity over multiple scales, EEG recording is decomposed by applying the empirical mode decomposition (EMD) which is known to be effective for extracting the constituent narrowband components without a predetermined basis. Following calculation of Renyi entropy of the probability distributions of the intrinsic mode functions extracted by EMD leads to a data-driven multiscale Renyi entropy. To validate the performance of the proposed entropy measure, actual EEG recordings from rats (n = 9) experiencing 7 min cardiac arrest followed by resuscitation were analyzed. Simulation and experimental results demonstrate that the use of the multiscale Renyi entropy leads to better discriminative capability of the injury levels and improved correlations with the neurological deficit evaluation after 72 hours after cardiac arrest, thus suggesting an effective diagnostic and prognostic tool.

Parameter selection in permutation entropy for an electroencephalographic measure of isoflurane anesthetic drug effect

The permutation entropy (PE) of the electroencephalographic (EEG) signals has been proposed as a robust measure of anesthetic drug effect. The calculation of PE involves the somewhat arbitrary selection of embedding dimension (m) and lag (τ) parameters. Previous studies of PE include the analysis of EEG signals under sevoflurane or propofol anesthesia, where different parameter settings were determined using a number of different criteria. In this study we choose parameter values based on the quantitative performance, to quantify the effect of a wide range of concentrations of isoflurane on the EEG. We analyzed a set of previously published EEG data, obtained from 29 patients who underwent elective abdominal surgery under isoflurane general anesthesia combined with epidural anesthesia. PE indices using a range of different parameter settings (m = 3-7, τ = 1-5) were calculated. These indices were evaluated as regards: the correlation coefficient (r) with isoflurane end-tidal concentration, the relationship with isoflurane effect-site concentration assessed by the coefficient of determination (R (2)) of the pharmacokinetic-pharmacodynamic models, and the prediction probability (PK). The embedding dimension (m) and lag (τ) have significant effect on the r values (Two-way repeated-measures ANOVA, p < 0.001). The proposed new permutation entropy index (NPEI) [a combination of PE(m = 3, τ = 2) and PE(m = 3, τ = 3)] performed best among all the parameter combinations, with r = 0.89(0.83-0.94), R (2) = 0.82(0.76-0.87), and PK = 0.80 (0.76-0.85). Further comparison with previously suggested PE measures, as well as other unrelated EEG measures, indicates the superiority of the NPEI. The PE can be utilized to indicate the dynamical changes of EEG signals under isoflurane anesthesia. In this study, the NPEI measure that combines the PE with m = 3, τ = 2 and that with m = 3, τ = 3 is optimal.

Measuring the effects of sevoflurane on electroencephalogram using sample entropy

BACKGROUND

Monitoring the effect of anesthetic drugs on the neural system is a major ongoing challenge for anesthetists. During the past few years, several electroencephalogram (EEG)-based methods such as the response entropy (RE) as implemented in the Datex-Ohmeda M-Entropy Module have been proposed. In this paper, sample entropy is used to quantify the predictability of EEG series, which could provide an index to show the effect of sevoflurane anesthesia. The dose-response relation of sample entropy is compared with that of RE.

METHODS

EEG data from 21 subjects is collected during the induction of general anesthesia with sevoflurane. The sample entropy is applied to the EEG recording. Pharmacokinetic-pharmacodynamic modeling and prediction probability statistic are used to evaluate the efficiency of sample entropy in comparison with RE.

RESULTS

Both methods track the gross changes in EEG, especially the occurrence of burst-suppression pattern at high doses of anesthetics. However, our method produces faster reaction to transients in EEG during the induction of anesthesia as indicated from the pharmacokinetic and pharmacodynamic modeled parameters and analysis around the point of loss of consciousness. Also, sample entropy correlated more closely with effect-site sevoflurane concentration than the RE. In addition, our proposed method exhibits greater resistance to noise in the EEG signals.

CONCLUSION

The results demonstrate that sample entropy can estimate the sevoflurane drug effect on the EEG more effectively than the commercial RE index with a stronger noise resistance.

GPGPU-aided ensemble empirical-mode decomposition for EEG analysis during anesthesia

Ensemble empirical-mode decomposition (EEMD) is a novel adaptive time-frequency analysis method, which is particularly suitable for extracting useful information from noisy nonlinear or nonstationary data. Unfortunately, since the EEMD is highly compute-intensive, the method does not apply in real-time applications on top of commercial-off-the-shelf computers. Aiming at this problem, a parallelized EEMD method has been developed using general-purpose computing on the graphics processing unit (GPGPU), namely, G-EEMD. A spectral entropy facilitated by G-EEMD was, therefore, proposed to analyze the EEG data for estimating the depth of anesthesia (DoA) in a real-time manner. In terms of EEG data analysis, G-EEMD has dramatically improved the run-time performance by more than 140 times compared to the original serial EEMD implementation. G-EEMD also performs far better than another parallelized implementation of EEMD bases on conventional CPU-based distributed computing technology despite the latter utilizes 16 high-end computing nodes for the same computing task. Furthermore, the results obtained from a pharmacokinetics/pharmacodynamic (PK/PD) model analysis indicate that the EEMD method is slightly more effective than its precedent alternative method (EMD) in estimating DoA, the coefficient of determination R(2) by EEMD is significantly higher than that by EMD (p < 0.05, paired t-test) and the prediction probability P(k) by EEMD is also slighter higher than that by EMD (p < 0.2, paired t-test).

Anaesthesia monitoring by recurrence quantification analysis of EEG data

Appropriate monitoring of the depth of anaesthesia is crucial to prevent deleterious effects of insufficient anaesthesia on surgical patients. Since cardiovascular parameters and motor response testing may fail to display awareness during surgery, attempts are made to utilise alterations in brain activity as reliable markers of the anaesthetic state. Here we present a novel, promising approach for anaesthesia monitoring, basing on recurrence quantification analysis (RQA) of EEG recordings. This nonlinear time series analysis technique separates consciousness from unconsciousness during both remifentanil/sevoflurane and remifentanil/propofol anaesthesia with an overall prediction probability of more than 85%, when applied to spontaneous one-channel EEG activity in surgical patients.

Hilbert-Huang versus Morlet wavelet transformation on mismatch negativity of children in uninterrupted sound paradigm

BACKGROUND

Compared to the waveform or spectrum analysis of event-related potentials (ERPs), time-frequency representation (TFR) has the advantage of revealing the ERPs time and frequency domain information simultaneously. As the human brain could be modeled as a complicated nonlinear system, it is interesting from the view of psychological knowledge to study the performance of the nonlinear and linear time-frequency representation methods for ERP research. In this study Hilbert-Huang transformation (HHT) and Morlet wavelet transformation (MWT) were performed on mismatch negativity (MMN) of children. Participants were 102 children aged 8-16 years. MMN was elicited in a passive oddball paradigm with duration deviants. The stimuli consisted of an uninterrupted sound including two alternating 100 ms tones (600 and 800 Hz) with infrequent 50 ms or 30 ms 600 Hz deviant tones. In theory larger deviant should elicit larger MMN. This theoretical expectation is used as a criterion to test two TFR methods in this study. For statistical analysis MMN support to absence ratio (SAR) could be utilized to qualify TFR of MMN.

RESULTS

Compared to MWT, the TFR of MMN with HHT was much sharper, sparser, and clearer. Statistically, SAR showed significant difference between the MMNs elicited by two deviants with HHT but not with MWT, and the larger deviant elicited MMN with larger SAR.

CONCLUSION

Support to absence ratio of Hilbert-Huang Transformation on mismatch negativity meets the theoretical expectations, i.e., the more deviant stimulus elicits larger MMN. However, Morlet wavelet transformation does not reveal that. Thus, HHT seems more appropriate in analyzing event-related potentials in the time-frequency domain. HHT appears to evaluate ERPs more accurately and provide theoretically valid information of the brain responses.

Using recurrence plot for determinism analysis of EEG recordings in genetic absence epilepsy rats

OBJECTIVE

Understanding the transition of brain activity towards an absence seizure is a challenging task. In this paper, we use recurrence quantification analysis to indicate the deterministic dynamics of EEG series at the seizure-free, pre-seizure and seizure states in genetic absence epilepsy rats.

METHODS

The determinism measure, DET, based on recurrence plot, was applied to analyse these three EEG datasets, each dataset containing 300 single-channel EEG epochs of 5-s duration. Then, statistical analysis of the DET values in each dataset was carried out to determine whether their distributions over the three groups were significantly different. Furthermore, a surrogate technique was applied to calculate the significance level of determinism measures in EEG recordings.

RESULTS

The mean (+/-SD) DET of EEG was 0.177+/-0.045 in pre-seizure intervals. The DET values of pre-seizure EEG data are significantly higher than those of seizure-free intervals, 0.123+/-0.023, (P<0.01), but lower than those of seizure intervals, 0.392+/-0.110, (P<0.01). Using surrogate data methods, the significance of determinism in EEG epochs was present in 25 of 300 (8.3%), 181 of 300 (60.3%) and 289 of 300 (96.3%) in seizure-free, pre-seizure and seizure intervals, respectively.

CONCLUSIONS

Results provide some first indications that EEG epochs during pre-seizure intervals exhibit a higher degree of determinism than seizure-free EEG epochs, but lower than those in seizure EEG epochs in absence epilepsy.

SIGNIFICANCE

The proposed methods have the potential of detecting the transition between normal brain activity and the absence seizure state, thus opening up the possibility of intervention, whether electrical or pharmacological, to prevent the oncoming seizure.

Synchronization Measurement of Multiple Neuronal Populations

The purpose of the present paper is to develop a method, based on equal-time correlation, correlation matrix analysis and surrogate resampling, that is able to quantify and describe properties of synchronization of population neuronal activity recorded simultaneously from multiple sites. Initially, Lorenz-type oscillators were used to model multiple time series with different patterns of synchronization. Eigenvalue and eigenvector decomposition was then applied to identify “clusters” of locally synchronized activity and to calculate a “global synchronization index.” This method was then applied to multichannel data recorded from an in vitro model of epileptic seizures. The results demonstrate that this novel method can be successfully used to analyze synchronization between multiple neuronal population series.