Stefanos Georgiadis

Single-trial dynamical estimation of event-related potentials: a Kalman filter-based approach

A method for single-trial dynamical estimation of event-related potentials (ERPs) is presented. The method is based on recursive Bayesian mean square estimation and the estimators are obtained with a Kalman filtering procedure. We especially focus on the case that previous trials contain prior information of relevance to the trial being analyzed. The potentials are estimated sequentially using the previous estimates as prior information. The performance of the method is evaluated with simulations and with real P300 responses measured using auditory stimuli. Our approach is shown to have excellent capability of estimating dynamic changes form stimulus to stimulus present in the parameters of the ERPs, even in poor signal-to-noise ratio (SNR) conditions.

Time-varying analysis of heart rate variability signals with a Kalman smoother algorithm

A time-varying parametric spectrum estimation method for analysing non-stationary heart rate variability signals is presented. As a case study, the dynamics of heart rate variability during an orthostatic test is examined. In this method, the non-stationary signal is first modelled with a time-varying autoregressive model and the model parameters are estimated recursively with a Kalman smoother algorithm. The benefit of using the Kalman smoother is that the lag error present in a Kalman filter, as well as in all other adaptive filters, can be avoided. The spectrum estimates for each time instant are then obtained from the estimated model parameters. Statistics of the obtained spectrum estimates are derived using the error propagation principle. The obtained spectrum estimates can further be decomposed into separate components and, thus, the time variation of low- and high-frequency components of heart rate variability can be examined separately. By using the presented method, high resolution time-varying spectrum estimates with no lag error can be produced. Other benefits of the method are the straightforward procedure for evaluating the statistics of the spectrum estimates and the option of spectral decomposition.

Linear and nonlinear tremor acceleration characteristics in patients with Parkinson's disease

The purpose of the study was to evaluate linear and nonlinear tremor characteristics of the hand in patients with Parkinsons disease (PD) and to compare the results with those of healthy old and young control subjects. Furthermore, the aim was to study correlation between tremor characteristics and clinical signs. A variety of nonlinear (sample entropy, cross-sample entropy, recurrence rate, determinism and correlation dimension) and linear (amplitude, spectral peak frequency and total power, and coherence) hand tremor parameters were computed from acceleration measurements for PD patients (n = 30, 68.3 ± 7.8 years), and old (n = 20, 64.2 ± 7.0 years) and young (n = 20, 18.4 ± 1.1 years) control subjects. Nonlinear tremor parameters such as determinism, sample entropy and cross-sample entropy were significantly different between the PD patients and healthy controls. These parameters correlated with the Unified Parkinsons disease rating scale (UPDRS), tremor and finger tapping scores, but not with the rigidity scores. Linear tremor parameters such as the amplitude and the maximum power (power corresponding to peak frequency) also correlated with the clinical findings. No major difference was detected in the tremor characteristics between old and young control subjects. The study revealed that tremor in PD patients is more deterministic and regular when compared to old or young healthy controls. The nonlinear tremor parameters can differentiate patients with PD from healthy control subjects and these parameters may have potential in the assessment of the severity of PD (UPDRS).

Time-Varying Analysis of Heart Rate Variability with Kalman Smoother Algorithm

A time-varying parametric spectrum estimation method for analyzing nonstationary heart rate variability signals is presented. In the method, the nonstationary signal is first modeled with time-varying autoregressive model and the model parameters are estimated recursively with a Kalman smoother algorithm. The spectrum estimates for each time are then obtained from the estimated model parameters. Statistics of the obtained spectrum estimates are derived using the error propagation principle. The obtained spectrum estimates can further be decomposed into separate components and, thus, the time-variation of low and high frequency components of heart rate variability can be examined separately

Directional Connectivity between Frontal and Posterior Brain Regions Is Altered with Increasing Concentrations of Propofol

Recent studies using electroencephalography (EEG) suggest that alteration of coherent activity between the anterior and posterior brain regions might be used as a neurophysiologic correlate of anesthetic-induced unconsciousness. One way to assess causal relationships between brain regions is given by renormalized partial directed coherence (rPDC). Importantly, directional connectivity is evaluated in the frequency domain by taking into account the whole multichannel EEG, as opposed to time domain or two channel approaches. rPDC was applied here in order to investigate propofol induced changes in causal connectivity between four states of consciousness: awake (AWA), deep sedation (SED), loss (LOC) and return of consciousness (ROC) by gathering full 10/20 system human EEG data in ten healthy male subjects. The target-controlled drug infusion was started at low rate with subsequent gradual stepwise increases at 10 min intervals in order to carefully approach LOC (defined as loss of motor responsiveness to a verbal stimulus). The direction of the causal EEG-network connections clearly changed from AWA to SED and LOC. Propofol induced a decrease (p = 0.002–0.004) in occipital-to-frontal rPDC of 8-16 Hz EEG activity and an increase (p = 0.001–0.040) in frontal-to-occipital rPDC of 10–20 Hz activity on both sides of the brain during SED and LOC. In addition, frontal-to-parietal rPDC within 1–12 Hz increased in the left hemisphere at LOC compared to AWA (p = 0.003). However, no significant changes were detected between the SED and the LOC states. The observed decrease in back-to-front EEG connectivity appears compatible with impaired information flow from the posterior sensory and association cortices to the executive prefrontal areas, possibly related to decreased ability to perceive the surrounding world during sedation. The observed increase in the opposite (front-to-back) connectivity suggests a propofol concentration dependent association and is not directly related to the level of consciousness per se.

Electroencephalogram reactivity to verbal command after dexmedetomidine, propofol and sevoflurane‐induced unresponsiveness

Although electroencephalogram reactivity (i.e. transient changes in electrical brain activity following external stimulus) might be useful in depth‐of‐anaesthesia monitoring, it has not been systematically examined with different anaesthetics at doses titrated to unresponsiveness. Three 10‐subject groups of healthy volunteers received dexmedetomidine, propofol or sevoflurane in escalating pseudo‐steady‐state concentrations at 10‐min intervals until they did not open their eyes to command. The electroencephalogram was continuously recorded and spectral variables were calculated with short‐time Fourier transform and time‐varying autoregressive modelling. Electroencephalogram reactivity was most prominent in the midfrontal derivations (termed F3 and F4). During drug‐induced unresponsiveness, electroencephalogram reactivity was still present in all drug groups. Dexmedetomidine, propofol and sevoflurane induced distinct suppression patterns on the electroencephalogram reactivity at the same clinical endpoint (unresponsiveness). Reactivity was best maintained with propofol, while only minimally preserved with dexmedetomidine and sevoflurane. Thus, it may be difficult to harness reactivity for depth‐of‐anaesthesia monitoring.

Time-varying spectrum estimation of heart rate variability signals with Kalman smoother algorithm

A time-varying parametric spectrum estimation method for analyzing dynamics of heart rate variability (HRV) signals is presented. In the method, HRV signal is first modeled with a time-varying autoregressive model and the model parameters are solved recursively with a Kalman smoother algorithm. Time-varying spectrum estimates are then obtained from the estimated model parameters. The obtained spectrum can be further decomposed into separate components, which is especially advantageous in HRV applications where low frequency (LF) and high frequency (HF) components are generally aimed to be distinguished. As case studies, the dynamics of HRV signals recorded during 1) orthostatic test, 2) exercise test and 3) simulated driving task are analyzed.

A subspace method for dynamical estimation of evoked potentials

It is a challenge in evoked potential (EP) analysis to incorporate prior physiological knowledge for estimation. In this paper, we address the problem of single-channel trial-to-trial EP characteristics estimation. Prior information about phase-locked properties of the EPs is assesed by means of estimated signal subspace and eigenvalue decomposition. Then for those situations that dynamic fluctuations from stimulus-to-stimulus could be expected, prior information can be exploited by means of state-space modeling and recursive Bayesian mean square estimation methods (Kalman filtering and smoothing). We demonstrate that a few dominant eigenvectors of the data correlation matrix are able to model trend-like changes of some component of the EPs, and that Kalman smoother algorithm is to be preferred in terms of better tracking capabilities and mean square error reduction. We also demonstrate the effect of strong artifacts, particularly eye blinks, on the quality of the signal subspace and EP estimates by means of independent component analysis applied as a prepossessing step on the multichannel measurements.

Analysis of heart rate variability dynamics during propofol and dexmedetomidine anesthesia

It has been observed that heart rate variability (HRV) diminishes during anesthesia, but the exact mechanisms causing it are not completely understood. The aim of this paper was to study the dynamics of HRV during low dose propofol (N=9) and dexmedetomidine (N=8) anesthesia by using state-of-the-art time-varying methods, and thereby ultimately try to improve the safety of anesthesia. The time-varying spectrum is estimated by using a Kalman smoother approach. The results show that there is an overall increase in HRV and decrease in heart rate prior to loss of consciousness. For dexmedetomidine these changes are more considerable than for propofol. For dexmedetomidine the variability also seems to start decreasing right after loss of consciousness, whereas for propofol HRV continues increasing.

Nonlinear parameters of surface EMG in schizophrenia patients depend on kind of antipsychotic therapy

We compared a set of surface EMG (sEMG) parameters in several groups of schizophrenia (SZ, n = 74) patients and healthy controls (n = 11) and coupled them with the clinical data. sEMG records were quantified with spectral, mutual information (MI) based and recurrence quantification analysis (RQA) parameters, and with approximate and sample entropies (ApEn and SampEn). Psychotic deterioration was estimated with Positive and Negative Syndrome Scale (PANSS) and with the positive subscale of PANSS. Neuroleptic-induced parkinsonism (NIP) motor symptoms were estimated with Simpson-Angus Scale (SAS). Dyskinesia was measured with Abnormal Involuntary Movement Scale (AIMS). We found that there was no difference in values of sEMG parameters between healthy controls and drug-naïve SZ patients. The most specific group was formed of SZ patients who were administered both typical and atypical antipsychotics (AP). Their sEMG parameters were significantly different from those of SZ patients taking either typical or atypical AP or taking no AP. This may represent a kind of synergistic effect of these two classes of AP. For the clinical data we found that PANSS, SAS, and AIMS were not correlated to any of the sEMG parameters. Conclusion: with nonlinear parameters of sEMG it is possible to reveal NIP in SZ patients, and it may help to discriminate between different clinical groups of SZ patients. Combined typical and atypical AP therapy has stronger effect on sEMG than a therapy with AP of only one class.