The aim of the present study was to analyse the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD), one of the most frequent disorders among elderly population. For this pilot study, we recorded the MEGs with a 148-channel whole-head magnetometer in 20 patients with probable AD and 21 age-matched control subjects. Artefact-free epochs of 3392 samples were analysed with auto-mutual information (AMI). Average AMI decline rates were lower for the AD patients' recordings than for control subjects' ones. Statistically significant differences were found using a Student's t-test (p<0.01) in 144 channels. Mean AMI values were analysed with a receiver operating characteristic curve. Sensitivity, specificity and accuracy values of 75%, 90.5% and 82.9% were obtained. Our results show that AMI estimations of the magnetic brain activity are different in both groups, hence indicating an abnormal type of dynamics associated with AD. This study suggests that AMI might help medical doctors in the diagnosis of the disease.
We appreciate the interest of Dr Tang in our article. Certainly, our previous results should be taken with caution due to the small database size. Nevertheless, it must be noted that this limitation was clearly recognized in our article. Furthermore, our hypothesis is completely justified from the current state of the art in the analysis of electroencephalogram (EEG) signals from Alzheimer's disease (AD) patients. We evaluated whether the multiscale entropy (MSE) analysis of the EEG background activity was useful to distinguish AD patients and controls. We do believe that further discussions about risk factors or related clinicophysiological protein aspects are clearly beyond the scope of our article. For the sake of completeness, we now detail some results that complement our previous analysis. They suggest that the MSE analysis can provide relevant information about the dynamics of AD patients' EEG data. Thus, we must reaffirm our conclusions, although we again acknowledge that further studies are needed.
Santamarta D, Abásolo D, Martínez-Madrigal M, Hornero R (2012) Characterisation of the intracranial pressure waveform during infusion studies by means of central tendency measure., Acta Neurochir (Wien)154(9)pp. 1595-1602
In the present study an attempt was made to quantify and characterise the changes in the intracranial pressure (ICP) waveform over the wide pressure range covered during infusion studies by means of the central tendency measure (CTM). CTM is a non-linear approach using continuous chaotic modelling that summarises the degree of variability in a signal.
Hydrocephalus includes a number of disorders characterised by clinical symptoms, enlarged ventricles (observable using neuroimaging techniques) and altered cerebrospinal fluid (CSF) dynamics. Infusion tests are one of the available procedures to study CSF circulation in patients with clinical and radiological features of hydrocephalus. In them, intracranial pressure (ICP) is deliberately raised and CSF circulation disorders evaluated through measurements of the resulting ICP. In this study, we analysed seventy-seven ICP signals recorded during infusion tests using four spectral-based parameters: median frequency (MF) and relative power (RP) in three frequency bands. These measures provide a novel perspective for the analysis of ICP signals in the frequency domain. Each signal was divided into four artefact-free epochs (corresponding to the basal, early infusion, plateau and recovery phases of the infusion study). The four spectral parameters were calculated for each epoch. We analysed differences between epochs of the infusion test and correlations between these epochs and patient data. Statistically significant differences (p<1.7×10-3, Bonferroni-corrected Wilcoxon signed-rank tests) were found between epochs of the infusion test using MF and RP. Furthermore, some spectral parameters (MF in the basal phase, RP for the first frequency band and in the early infusion phase, RP for the second frequency band and in all phases of the infusion study and RP in the third frequency band and in the basal phase) revealed significant correlations (p<0.01) between epochs of the infusion test and signal amplitude in the basal and plateau phases. Our results suggest that spectral analysis of ICP signals could be useful for understanding CSF dynamics in hydrocephalus. © 2013 IPEM.
The aim of this pilot study was to assess the usefulness of independent component analysis (ICA) to detect cardiac artifacts and power line interferences in magnetoencephalogram (MEG) recordings. We recorded MEG signals from six subjects with, a 148-channel whole-head magnetometer (MAGNES 2500 WH, 4D Neuroimaging). Epochs of 50 s with power line noise, cardiac, and ocular artifacts were selected for analysis. We applied a statistical criterion to determine the number of sources, and a robust ICA algorithm to decompose the MEG epochs. Skewness, kurtosis, and a spectral metric were used to mark the studied artifacts. We found that the power fine interference could be easily detected by its frequency characteristics. Moreover, skewness outperformed kurtosis when identifying the cardiac artifact.
The aim of this study was to improve the diagnosis of Alzheimer's disease (AD) patients applying a blind source separation (BSS) and component selection procedure to their magnetoencephalogram (MEG) recordings. MEGs from 18 AD patients and 18 control subjects were decomposed with the algorithm for multiple unknown signals extraction. MEG channels and components were characterized by their mean frequency, spectral entropy, approximate entropy, and Lempel-Ziv complexity. Using Student's t-test, the components which accounted for the most significant differences between groups were selected. Then, these relevant components were used to partially reconstruct the MEG channels. By means of a linear discriminant analysis, we found that the BSS-preprocessed MEGs classified the subjects with an accuracy of 80.6%, whereas 72.2% accuracy was obtained without the BSS and component selection procedure.
García Gadañón M, Poza J, Hornero R, Santamarta D, Abásolo D (2014) Analysis of Intracranial Pressure Signals Using the Spectral Turbulence, IFMBE Proceedings41pp. 795-798
Hydrocephalus includes a range of conditions characterized by clinical symptoms, enlarged ventricles and disorders in cerebrospinal fluid (CSF) circulation. Infusion tests can be used to analyze CSF dynamics in patients with hydrocephalus. In infusion tests, intracranial pressure (ICP) is artificially raised, while the resulting ICP is measured in order to detect CSF circulation alterations. In this study, we analyzed 77 ICP signals recorded during infusion tests using the spectral turbulence (ST). Each signal was divided into four artifact-free epochs. The mean ST, , and the standard deviation of ST, SD[ST], were calculated for each epoch. Statistically significant differences were found between the basal phase of the infusion test and the remaining phases using and SD[ST] (p<1.7·10, Bonferroni-corrected Wilcoxon tests). Furthermore, we found significantly higher and significantly lower SD[ST] values in the plateau phase than in the basal phase. These findings suggest that the increase in ICP induced by infusion studies is associated with a significant loss of irregularity and variability of the spectral content of ICP signals. In conclusion, spectral analysis of ICP signals could be useful for understanding CSF dynamics in hydrocephalus. © Springer International Publishing Switzerland 2014.
Objective: Considering the increasing use of complexity estimates in neuropsychiatric populations, a normative study is critical to define the 'normal' behaviour of brain oscillatory complexity across the life span. Method: This study examines changes in resting-state magnetoencephalogram (MEG) complexity - quantified with the Lempel-Ziv complexity (LZC) algorithm - due to age and gender in a large sample of 222 (100 males/122 females) healthy participants with ages ranging from 7 to 84 years. Results: A significant quadratic (curvilinear) relationship (p < 0.05) between age and complexity was found, with LZC maxima being reached by the sixth decade of life. Once that peak was crossed, complexity values slowly decreased until late senescence. Females exhibited higher LZC values than males, with significant differences in the anterior, central and posterior regions (p < 0.05). Conclusions: These results suggest that the evolution of brain oscillatory complexity across the life span might be considered a new illustration of a 'normal' physiological rhythm. Significance: Previous and forthcoming clinical studies using complexity estimates might be interpreted from a more complete and dynamical perspective. Pathologies not only cause an 'abnormal' increase or decrease of complexity values but they actually 'break' the 'normal' pattern of oscillatory complexity evolution as a function of age. © 2012 International Federation of Clinical Neurophysiology.
Hydrocephalus is a condition characterized by altered cerebrospinal fluid (CSF) dynamics and chronic rises in intracranial pressure (ICP). However, the reason why hydrocephalic physiologies fail to inhibit dangerously high ICP levels is not known. Infusion studies are used to raise ICP and evaluate CSF circulation disorders. In this pilot study, ICP signals recorded during infusion tests from 33 patients with normal pressure hydrocephalus and 36 patients having developed a secondary form of normal pressure hydrocephalus were characterized using Permutation Entropy (PE), a symbolic non-linear method to quantify complexity. Each ICP signal was divided into four epochs - baseline (before infusion begins), infusion, plateau, and recovery (after infusion has stopped) - and the mean PE was calculated for each epoch. Statistically significant differences were found between PE for most epochs (p<;0.00833, Bonferroni-corrected Wilcoxon tests), with a significant decrease in the plateau phase. However, differences between PE for normal pressure and secondary hydrocephalus were not significant. Results suggest that the increase in ICP during infusion studies is associated with a significant decrease in PE. PE analysis of ICP signals could be useful for increasing our understanding of CSF dynamics in normal pressure hydrocephalus.
Zamarrón C, Hornero R, del Campo F, Abásolo D, Alvarez D (2006) Heart rate regularity analysis obtained from pulse oximetric recordings in the diagnosis of obstructive sleep apnea., Sleep Breath10(2)pp. 83-89
Approximate entropy (ApEn) is a technique that can be used to quantify the irregularity or variability of time series. We prospectively evaluated the validity of ApEn of heart rate data obtained from pulse oximetric recordings as a diagnostic test for obstructive sleep apnea (OSA) in patients clinically suspected of suffering this disease. A sample of 187 referred outpatients (147 men and 40 women), with a mean age of 57.9+/-12.8 years and a body mass index of 29.5+/-5.5 kg/m(2), clinically suspected of having OSA were studied using nocturnal pulse oximetric recording performed simultaneously with complete polysomnography. A diagnosis of OSA was confirmed in 111 (59.3%). Patients with OSA presented significantly higher ApEn levels than those without OSA (1.334+/-0.189 vs 1.167+/-0.182). Chronic obstructive pulmonary disease (COPD) was diagnosed for 42 patients. Among these patients, 22 (52.4%) were diagnosed with OSA. COPD patients with OSA showed significantly higher ApEn levels than COPD patients without OSA (1.337+/-0.193 vs 1.184+/-0.173; p=0.01). ApEn correlated significantly with apnea-hypopnea index (r=0.38; p=0.000). There was no significant correlation between ApEn and either age or body mass index. No significant changes were observed in ApEn throughout the night in OSA patients. Using receiver operating characteristic curve analysis, we obtained a diagnostic sensitivity of 71.2%, specificity of 78.9%, positive predictive value of 81.3%, and negative predictive value of 66% at a threshold of 1.272. We conclude that ApEn analysis of heart rate data obtained from pulse oximetric recordings could be a useful tool in the study of OSA.
Hornero R, Aboy M, Abásolo D, McNames J, Goldstein B (2005) Interpretation of approximate entropy: analysis of intracranial pressure approximate entropy during acute intracranial hypertension., IEEE Trans Biomed Eng52(10)pp. 1671-1680
We studied changes in intracranial pressure (ICP) complexity, estimated by the approximate entropy (ApEn) of the ICP signal, as subjects progressed from a state of normal ICP (< 20-25 mmHg) to acutely elevated ICP (an ICP "spike" defined as ICP > 25 mmHg for < or = 5 min). We hypothesized that the measures of intracranial pressure (ICP) complexity and irregularity would decrease during acute elevations in ICP. To test this hypothesis we studied ICP spikes in pediatric subjects with severe traumatic brain injury (TBI). We conclude that decreased complexity of ICP coincides with episodes of intracranial hypertension (ICH) in TBI. This suggests that the complex regulatory mechanisms that govern intracranial pressure are disrupted during acute rises in ICP. Furthermore, we carried out a series of experiments where ApEn was used to analyze synthetic signals of different characteristics with the objective of gaining a better understanding of ApEn itself, especially its interpretation in biomedical signal analysis.
The magnetoencephalogram (MEG) is contaminated with undesired signals, which are called artifacts. Some of the most important ones are the cardiac and the ocular artifacts (CA and OA, respectively), and the power line noise (PLN). Blind source separation (BSS) has been used to reduce the influence of the artifacts in the data. There is a plethora of BSS-based artifact removal approaches, but few comparative analyses. In this study, MEG background activity from 26 subjects was processed with five widespread BSS (AMUSE, SOBI, JADE, extended Infomax, and FastICA) and one constrained BSS (cBSS) techniques. Then, the ability of several combinations of BSS algorithm, epoch length, and artifact detection metric to automatically reduce the CA, OA, and PLN were quantified with objective criteria. The results pinpointed to cBSS as a very suitable approach to remove the CA. Additionally, a combination of AMUSE or SOBI and artifact detection metrics based on entropy or power criteria decreased the OA. Finally, the PLN was reduced by means of a spectral metric. These findings confirm the utility of BSS to help in the artifact removal for MEG background activity.
Understanding the dynamics of brain activity manifested in the EEG, local field potentials (LFP), and neuronal spiking is essential for explaining their underlying mechanisms and physiological significance. Much has been learned about sleep regulation using conventional EEG power spectrum, coherence, and period-amplitude analyses, which focus primarily on frequency and amplitude characteristics of the signals and on their spatio-temporal synchronicity. However, little is known about the effects of ongoing brain state or preceding sleep-wake history on the nonlinear dynamics of brain activity. Recent advances in developing novel mathematical approaches for investigating temporal structure of brain activity based on such measures, as Lempel-Ziv complexity (LZC) can provide insights that go beyond those obtained with conventional techniques of signal analysis. Here, we used extensive data sets obtained in spontaneously awake and sleeping adult male laboratory rats, as well as during and after sleep deprivation, to perform a detailed analysis of cortical LFP and neuronal activity with LZC approach. We found that activated brain states-waking and rapid eye movement (REM) sleep are characterized by higher LZC compared with non-rapid eye movement (NREM) sleep. Notably, LZC values derived from the LFP were especially low during early NREM sleep after sleep deprivation and toward the middle of individual NREM sleep episodes. We conclude that LZC is an important and yet largely unexplored measure with a high potential for investigating neurophysiological mechanisms of brain activity in health and disease.
The spurious increase in coherence of electroencephalogram (EEG) signals between distant electrode points has long been understood to be due to volume conduction of the EEG signal. Reducing the volume conduction components of EEG recordings in pre-processing attenuates this. However, the effect of volume conduction on non-linear signal processing of EEG signals is yet to be fully described. This pilot study aimed to investigate the impact of volume conduction on results calculated with a distance based, bivariate form of Lempel-Ziv Complexity (dLZC) by analyzing EEG signals from Alzheimer's disease (AD) patients and healthy age-matched controls with and without pre-processing with Current Source Density (CSD) transformation. Spurious statistically significant differences between AD patients and control's EEG signals seen without CSD pre-processing were not seen with CSD volume conduction mitigation. There was, however, overlap in the region of electrodes which were seen to hold this statistically significant information. These results show that, while previously published findings are still valid, volume conduction mitigation is required to ensure non-linear signal processing methods identify changes in signals only due to the purely local signal alone.
Fernández A, Gómez C, Hornero R, Abásolo D, López M (2006) Complexity analysis of the magnetoencephalogram background activity in Alzheimer disease patients, J PSYCHOPHYSIOL20(3)pp. 229-229 HOGREFE & HUBER PUBLISHERS
In this study we have investigated the electroencephalogram (EEG) background activity in patients with Alzheimer's disease (AD) using non-linear analysis methods. We calculated the Lempel-Ziv (LZ) complexity - applying two different sequence conversion methods - and the central tendency measure (CTM) of the EEG in 11 AD patients and 11 age-matched control subjects. CTM quantifies the degree of variability, while LZ complexity reflects the arising rate of new patterns along with the EEG time series. We did not find significant differences between AD patients and control subjects' EEGs with CTM. On the other hand, AD patients had significantly lower LZ complexity values (p<0.01) at electrodes P3 and O1 with a two-symbol sequence conversion, and P3, P4, O1 and T5 using three symbols. Our results show a decreased complexity of EEG patterns in AD patients. In addition, we obtained 90.9% sensitivity and 72.7% specificity at O1, and 72.7% sensitivity and 90.9% specificity at P3 and P4. These findings suggest that LZ complexity may contribute to increase the insight into brain dysfunction in AD in ways which are not possible with more classical and conventional statistical methods.
Simons S, Abásolo D, Escudero J (2012) Fuzzy Entropy and Multiscale Fuzzy Entropy of the Electroencephalogram in Alzheimer's Disease, Book of abstracts of the Royal Academy of Engineering Young Researchers Futures Meeting ? Neural Engineering
Al-Amri M, Abasolo D, Ghoussayni S, Ewins D (2013) Virtual reality based lower extremity exercises for children with cerebral palsy: Preliminary work, 2013 International Conference on Virtual Rehabilitation, ICVR 2013pp. 198-199
The work outlined in this paper aimed to evaluate the practicality of the Surrey Virtual Rehabilitation System (SVRS) during simple hip and knee flexion-extension exercises. Two clinical engineering trainees, a physiotherapist and two children with CP performed two tests using the SVRS whilst parents and a paediatric physiotherapist observed the tests. Feedback collected using both questionnaires and an open-discussion suggests that the SVRS is practical for rehabilitation purposes during both physical exercises. © 2013 IEEE.
This study assesses the connectivity alterations caused by Alzheimer's disease (AD) and mild cognitive impairment (MCI) in magnetoencephalogram (MEG) background activity. Moreover, a novel methodology to adaptively extract brain rhythms from the MEG is introduced. This methodology relies on the ability of empirical mode decomposition to isolate local signal oscillations and constrained blind source separation to extract the activity that jointly represents a subset of channels. Inter-regional MEG connectivity was analysed for 36 AD, 18 MCI and 26 control subjects in ´, ¸, ± and ² bands over left and right central, anterior, lateral and posterior regions with magnitude squared coherence?c(f). For the sake of comparison, c(f) was calculated from the original MEG channels and from the adaptively extracted rhythms. The results indicated that AD and MCI cause slight alterations in the MEG connectivity. Computed from the extracted rhythms, c(f) distinguished AD and MCI subjects from controls with 69.4% and 77.3% accuracies, respectively, in a full leave-one-out cross-validation evaluation. These values were higher than those obtained without the proposed extraction methodology.
Alzheimer's disease (AD) is one of the most prominent neurodegenerative disorders. The aim of this study is to analyze the magnetoencephalogram (MEG) background activity in AD patients using sample entropy (SampEn) and multiscale entropy (MSE). The former quantifies the signal regularity, while the latter is a complexity measure. These concepts, irregularity and complexity, are linked although the relationship is not straightforward. Five minutes of recording were acquired with a 148-channel whole-head magnetometer in 20 patients with probable AD and 21 control subjects. Our results show that MEG recordings are less complex and more regular in AD patients than in control subjects. Significant differences between both groups were found in some MEG channels with both methods (p<0.01, Student's t-test with Bonferroni's correction). Using receiver operating characteristic curves, accuracies of 75.6% with SampEn and of 87.8% with MSE were reached. Our findings show the usefulness of these entropy measures to increase our insight into AD.
In this study, we explored the ability of several spectral based measures to summarize the information of the power spectral density (PSD) function from spontaneous magnetoencephalographic (MEG) activity in Alzheimer's disease (AD). The MEGs of 20 AD patients and 21 elderly controls were recorded with eyes closed at rest during 5 min from 148 channels. Five spectral parameters were estimated from PSD: mean frequency (MF), individual alpha frequency (IAF), transition frequency (TF), 95% spectral edge frequency (SEF95) and spectral entropy (SE). To reduce the dimensionality of the problem, we applied a principal component analysis. According to our results, MF was the best discriminating index between both groups (85.00% sensitivity, 85.71% specificity) indicating a shift to the left of the power spectrum in AD. A significant MEG slowing was also observed using both IAF (p < 0.001) and TF (p < 0.01). The lowest classification statistics (65% sensitivity, 66.67% specificity) were obtained with SEF95. However, these results were also significant (p < 0.05). This fact points out that there is a variation in the spectral content at high frequencies of AD patients and controls. Finally, a significant decrease of irregularity in the AD group was observed with SE, with results close to those obtained with MF (90.00% sensitivity, 76.19% specificity). In conclusion, a complete description of PSD can help to increase our insight into brain dysfunction in AD and to extract spectral patterns specific to the disease.
The aim of the present research is to study the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD) using the Lempel-Ziv (LZ) complexity. We recorded the MEG with a 148-channel whole-head magnetometer (MAGNES 2500 WH, 4D Neuroimaging) in 10 patients with probable AD and 10 age-matched control subjects, during five minutes. Artefact-free epochs were selected for the non-linear analysis. In all MEG channels, the AD patients had lower complexity than control subjects. In 77 of them the differences were statistically significant (p < 0.01). These preliminary results suggest that cognitive dysfunction in AD is associated with a decreased complexity in certain regions of the brain.
The auto mutual information function (AMIF) evaluates the signal predictability by assessing linear and non-linear dependencies between two measurements taken from a single time series. Furthermore, the AMIF rate of decrease (AMIFRD) is correlated with signal entropy. This metric has been used to analyze biomedical data, including cardiac and brain activity recordings. Hence, the AMIFRD can be a relevant parameter in the context of biomedical signal analysis. Thus, in this pilot study, we have analyzed a synthetic sequence (a Lorenz system) and real biosignals (electroencephalograms recorded with eyes open and closed) with the AMIFRD. We aimed at illustrating the application of this parameter to biomedical time series. Our results show that the AMIFRD can detect changes in the non-linear dynamics of a sequence and that it can distinguish different physiological conditions.
Simons S, Abásolo D (2014) Permutation entropy of the electroencephalogram background activity in Alzheimer's disease: Investigation into the incidence of repeated values, BIOSIGNALS 2014 - 7th Int. Conference on Bio-Inspired Systems and Signal Processing, Proceedings; Part of 7th Int. Joint Conference on Biomedical Engineering Systems and Technologies, BIOSTEC 2014pp. 97-103
This pilot study applied Permutation Entropy (PE), a non-linear symbolic measure, and a novel modification (modPE), to investigate the regularity of electroencephalogram (EEG) signals from 11 Alzheimer's disease (AD) patients and 11 age-matched controls given input parameters n (embedding vector), Ä (coarse graining) and slide (difference between the start of two concurrent embedding vectors). PE discriminated better than modPE with controls showing reduced regularity over AD patients. Increasing Ä identified the greatest differences between EEG signals. Longer embedding vectors were also more able to identify differences. The greatest difference between groups was at Fp1 with n,Ä,slide = 3,10,1 (p=0.0112 Kruskal Wallis with Bonferroni). Subject and epoch based leave-one-out cross validation was carried out with thresholding from Receiver Operating Characteristic Curves. The greatest ability to correctly identify AD patients and controls were 81.82% (Fp2 n,Ä,slide = 7,4,4, PE and modPE, F7 n,Ä,slide = 3,10,1, PE and modPE) and 90.91% (Fp1 n,Ä,slide = 3,10,1, PE and modPE), respectively. The maximum accuracy (both groups correctly identified) was 81.82% seen at many electrode and input combinations. All are with subject based analysis. This suggests that PE can identify changes in EEG signals in AD, given appropriate variables. However, modPE makes little improvement over PE. Copyright © 2014 SCITEPRESS - Science and Technology Publications. All rights reserved.
Hornero R, Abásolo D, Jimeno N, Sánchez CI, Poza J, Aboy M (2006) Variability, regularity, and complexity of time series generated by schizophrenic patients and control subjects., IEEE Trans Biomed Eng53(2)pp. 210-218
We analyzed time series generated by 20 schizophrenic patients and 20 sex- and age-matched control subjects using three nonlinear methods of time series analysis as test statistics: central tendency measure (CTM) from the scatter plots of first differences of data, approximate entropy (ApEn), and Lempel-Ziv (LZ) complexity. We divided our data into a training set (10 patients and 10 control subjects) and a test set (10 patients and 10 control subjects). The training set was used for algorithm development and optimum threshold selection. Each method was assessed prospectively using the test dataset. We obtained 80% sensitivity and 90% specificity with LZ complexity, 90% sensitivity, and 60% specificity with ApEn, and 70% sensitivity and 70% specificity with CTM. Our results indicate that there exist differences in the ability to generate random time series between schizophrenic subjects and controls, as estimated by the CTM, ApEn, and LZ. This finding agrees with most previous results showing that schizophrenic patients are characterized by less complex neurobehavioral and neuropsychologic measurements.
The Lempel-Ziv (LZ) complexity is a non-linear time series analysis metric that reflects the arising rate of new patterns along with the sequence. Thus, it captures its temporal sequence and, quite conveniently, it can be computed with short data segments. In the present work, a detailed analysis on LZ complexity is presented within the context of atrial fibrillation (AF) organization estimation. As the analysed time series depend on the original sampling rate (fs), we evaluated the relationship between LZ complexity and fs. Furthermore, different implementations of LZ complexity were tested. Our results show the usefulness of LZ complexity to estimate AF organization and suggest that the signals from a terminating paroxysmal AF group are more organized (i.e. less complex) than those from the non-terminating paroxysmal AF group. However, the diagnostic accuracy was not as high as that obtained with sample entropy (SampEn), another non-linear metric, with the same database in a previous study (92% vs. 96%). Nevertheless, the LZ complexity analysis of AF organization with sampling frequencies higher than 2048 Hz, or even its combination with SampEn or other non-linear metrics, might improve the prediction of spontaneous AF termination.
Abásolo D, Winsky-Sommerer R (2013) Investigating the possible usefulness of Lempel-Ziv complexity for the characterisation of electroencephalogram sleep recordings, Final Programme and Book of Abstracts of BioDynamics 2013pp. 24-25
Approximate entropy (ApEn) is a family of statistics introduced as a quantification of regularity in time series without any a priori knowledge about the system generating them. The aim of this preliminary study was to assess whether a time series analysis of arterial oxygen saturation (SaO2) signals from overnight pulse oximetry by means of ApEn could yield essential information on the diagnosis of obstructive sleep apnea (OSA) syndrome. We analyzed SaO2 signals from 187 subjects: 111 with a positive diagnosis of OSA and 76 with a negative diagnosis of OSA. We divided our data in a training set (44 patients with OSA Positive and 30 patients with OSA Negative) and a test set (67 patients with OSA Positive and 46 patients with OSA Negative). The training set was used for algorithm development and optimum threshold selection. Results showed that recurrence of apnea events in patients with OSA determined a significant increase in ApEn values. This method was assessed prospectively using the test dataset, where we obtained 82.09% sensitivity and 86.96% specificity. We conclude that ApEn analysis of SaO2 from pulse oximetric recording could be useful in the study of OSA.
The goal of this study was to analyze the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD) using the auto mutual information (AMI). Applied to time series, AMI provides a measure of future points predictability from past points. Five minutes of recording were acquired with a 148-channel whole-head magnetometer (MAGNES 2500 WH, 4D neuroimaging) in 12 patients with probable AD and 12 elderly control subjects. Artifact-free epochs of 20 seconds (3392 points, sample frequency of 169.6 Hz) were selected for our study. Our results showed that the absolute values of the averaged decline rate of AMI were lower in AD patients than in control subjects for all channels. In addition, there were statistically significant differences (p<0.01, Student's t-test) in most channels. These preliminary results suggest that neuronal dysfunction in AD is associated with differences in the dynamical processes underlying the MEG recording.
Aboy M, Hornero R, Abásolo D, Alvarez D (2006) Interpretation of the Lempel-Ziv complexity measure in the context of biomedical signal analysis., IEEE Trans Biomed Eng53(11)pp. 2282-2288
Lempel-Ziv complexity (LZ) and derived LZ algorithms have been extensively used to solve information theoretic problems such as coding and lossless data compression. In recent years, LZ has been widely used in biomedical applications to estimate the complexity of discrete-time signals. Despite its popularity as a complexity measure for biosignal analysis, the question of LZ interpretability and its relationship to other signal parameters and to other metrics has not been previously addressed. We have carried out an investigation aimed at gaining a better understanding of the LZ complexity itself, especially regarding its interpretability as a biomedical signal analysis technique. Our results indicate that LZ is particularly useful as a scalar metric to estimate the bandwidth of random processes and the harmonic variability in quasi-periodic signals.
Abásolo D, Jelinek HF (2012) Lempel-Ziv complexity dynamics in early detection of cardiac autonomic neuropathy in diabetes, Proceedings of the 9th IASTED International Conference on Biomedical Engineering, BioMed 2012pp. 16-20
Early detection of cardiac autonomic neuropathy (CAN) in patients with diabetes mellitus (DM) is of prime importance, as it will facilitate the prevention of its serious consequences. In the present work, the non-linear dynamics of electrocardiogram (ECG) recordings in 41 Type 2 DM patients with early CAN and 40 controls without clinical signs and symptoms of CAN were analysed with different implementations of Lempel-Ziv (LZ) complexity. LZ complexity is a non-linear analysis method that estimates the complexity of time series of finite length and reflects the arising rate of new patterns along the sequence. Results suggest that ECG traces are less complex in patients with early CAN than in those with no CAN. Differences were statistically significant (p < 0.05, Kruskal-Wallis test) when the LZ complexity was implemented with a three symbol conversion and the mean used to define the thresholds. Furthermore, the discriminative abilities of the different LZ complexity implementations in the context of CAN were evaluated with ROC curves. Accuracies over 65% were obtained when the mean was used to define the thresholds, with a sensitivity of 75.61% with a two symbol conversion. Our results suggest that LZ complexity might be a useful tool for an early detection of CAN from ECG recordings. Nevertheless, further studies are needed to address the possible usefulness of this methodology in the characterisation and early detection of CAN in type 2 DM patients.
Hornero R, Abásolo D, Jimeno N, Espino P (2003) Applying Approximate Entropy and Central Tendency Measure to Analyze Time Series Generated by Schizophrenic Patients, Annual International Conference of the IEEE Engineering in Medicine and Biology - Proceedings3pp. 2447-2450
Hornero R, Aboy M, Abásolo D (2007) Analysis of intracranial pressure during acute intracranial hypertension using Lempel-Ziv complexity: further evidence., Med Biol Eng Comput45(6)pp. 617-620
We analyzed intracranial pressure (ICP) signals during periods of acute intracranial hypertension (ICH) using the Lempel-Ziv (LZ) complexity measure. Our results indicate the LZ complexity of ICP decreases during periods of ICH. The mean LZ complexity before ICH was 0.20+/-0.04, while the mean LZ complexity during ICH was 0.16+/-0.03 (p<0.05). The mean decrease of the LZ complexity values during the ICH episodes was 19.5%. Additionally, we present preliminary evidence suggesting that periods of ICH may be detectable from non-invasive signals coupled with ICP, such as pulse oximetry (SpO2).
The aim of the present study was to analyze resting-state brain activity in patients with Parkinson's disease (PD), a degenerative disorder of the nervous system. Magnetoencephalography (MEG) signals were recorded with a 151-channel whole-head radial gradiometer MEG system in 18 early-stage untreated PD patients and 20 age-matched control subjects. Artifact-free epochs of 4 s (1250 samples) were analyzed with Lempel-Ziv complexity (LZC), applying two- and three-symbol sequence conversion methods. The results showed that MEG signals from PD patients are less complex than control subjects' recordings. We found significant group differences (p-values <0.01) for the 10 major cortical areas analyzed (e.g., bilateral frontal, central, temporal, parietal, and occipital regions). In addition, using receiver-operating characteristic curves with a leave-one-out cross-validation procedure, a classification accuracy of 81.58% was obtained. In order to investigate the best combination of LZC results for classification purposes, a forward stepwise linear discriminant analysis with leave-one out cross-validation was employed. LZC results (three-symbol sequence conversion) from right parietal and temporal brain regions were automatically selected by the model. With this procedure, an accuracy of 84.21% (77.78% sensitivity, 90.0% specificity) was achieved. Our findings demonstrate the usefulness of LZC to detect an abnormal type of dynamics associated with PD.
In this pilot study, a symbolic sequence decomposition method was used in conjunction with Shannon?s entropy to investigate the changes in electroencephalogram signals of 11 patients with Alzheimer?s disease and 11 age-matched control subjects. Results were statistically analysed by student t-test and later classified with receiver operating curves. Statistically significant differences between both groups were found at electrodes Fp1, O2, P3, T4 and T5. Sensitivity (defined as percentages of correctly classified patients) and
specificity (defined as correctly classified controls) were evaluated using the receiver operating curves method. Accuracy of the methods was calculated according to sensitivity and specificity measures of electrodes showing statistically significant differences between the control group and Alzheimer's disease patients and ranged between 72.73-77.27%. These accuracy values were in agreement with previously published entropy studies on this data set. Although combining these methods did not provide any greater accuracy over previous findings, using a symbolic sequence decomposition method enhanced the data processing.
The goal of this study was to analyze the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD) using a regularity measure: approximate entropy (ApEn). This measure was computed for a broad band (0.5-40 Hz) as well as typical frequency bands (delta, theta, alpha, beta and gamma). Five minutes of recording were acquired with a 148-channel whole-head magnetometer in 15 patients with probable AD and 15 elderly control subjects. Our results showed that AD patients' MEGs were more regular than controls' recordings at all frequency bands, with the exception of beta. Additionally, there were statistically significant differences (p ) 0.01, Student's t-test) at the broad and delta bands. Using receiver operating characteristic curves, the highest accuracy (83.33%) was reached at delta band. These results suggest the usefulness of ApEn to gain a better understanding of dynamical processes underlying the MEG recording.
Alzheimer's disease (AD) is the most common neurodegenerative disorder. Although a definite diagnosis is only possible by necropsy, a differential diagnosis with other types of dementia and with major depression should be attempted. The aim of this study was to analyse the electroencephalogram (EEG) background activity of AD patients to test the hypothesis that the regularity of the AD patients' EEG is higher than that of age-matched controls. We recorded the EEG from 19 scalp electrodes in 11 AD patients and 11 age-matched controls. Two different methods were used to estimate the regularity of the EEG background activity: spectral entropy (SpecEn) and sample entropy (SampEn). We did not find significant differences between AD patients and control subjects' EEGs with SpecEn. On the other hand, AD patients had significantly lower SampEn values than control subjects (p < 0.01) at electrodes P3, P4, O1 and O2. Our results show an increase of EEG regularity in AD patients. These findings suggest that nonlinear analysis of the EEG with SampEn could yield essential information and may contribute to increasing the insight into brain dysfunction in AD in ways which are not possible with more classical and conventional statistical methods.
The aim of this study was to analyse the electroencephalogram (EEG) background activity in Alzheimer's disease (AD) with the Detrended Moving Average (DMA) method, a new approach to quantify correlation properties in non-stationary signals with underlying trends. EEGs were recorded from the 19 scalp loci of the international 10-20 system in 11 AD patients and 11 age-matched controls. Our results showed two scaling regions in all subjects' channels, with a clear bend when their corresponding slopes (alpha(1) and alpha(2)) were distinctly different. With the exception of electrode T4, the alpha(1) values were lower in control subjects than in AD patients, with significant differences at TS, P3, P4 and O1 (p < 0.01, Student's t-test). On the other hand, alpha(2) values were higher in control subjects than in AD patients, with significant differences only at F4. Furthermore, we evaluated the ability of alpha(2) to discriminate AD patients from control subjects at these electrodes using ROC plots. We obtained a maximum accuracy of 81.82% at O1 with alpha(1) and at F4 with alpha(2). These findings suggest that the scaling behaviour of the EEG is sensitive to AD and that the DMA method could help to increase our insight into brain dysfunction in AD.
Body temperature is a classical diagnostic tool for a number of diseases. However, it is usually employed as a plain binary classification function (febrile or not febrile), and therefore its diagnostic power has not been fully developed. In this paper, we describe how body temperature regularity can be used for diagnosis. Our proposed methodology is based on obtaining accurate long-term temperature recordings at high sampling frequencies and analyzing the temperature signal using a regularity metric (approximate entropy). In this study, we assessed our methodology using temperature registers acquired from patients with multiple organ failure admitted to an intensive care unit. Our results indicate there is a correlation between the patient's condition and the regularity of the body temperature. This finding enabled us to design a classifier for two outcomes (survival or death) and test it on a dataset including 36 subjects. The classifier achieved an accuracy of 72%.
The aim of this study was to assess whether independent component analysis (ICA) could be valuable to remove power line noise, cardiac, and ocular artifacts from magnetoencephalogram (MEG) background activity. The MEGs were recorded from 11 subjects with a 148-channel whole-head magnetometer. We used a statistical criterion to estimate the number of independent components. Then, a robust ICA algorithm decomposed the MEG epochs and several methods were applied to detect those artifacts. The whole process had been previously tested on synthetic data. We found that the line noise components could be easily detected by their frequency spectrum. In addition, the ocular artifacts could be identified by their frequency characteristics and scalp topography. Moreover, the cardiac artifact was better recognized by its skewness value than by its kurtosis one. Finally, the MEG signals were compared before and after artifact rejection to evaluate our method.
The aim of this study was to analyze the electroencephalogram (EEG) background activity in Alzheimer's disease (AD) with two non-linear methods: Approximate Entropy (ApEn) and Auto Mutual Information (AMI). ApEn quantifies the regularity in data, while AMI detects linear and non-linear dependencies in time series. EEGs were recorded from the 19 scalp loci of the international 10-20 system in 11 AD patients and 11 age-matched controls. ApEn was significantly lower in AD patients at electrodes O1, O2, P3 and P4 (p<0.01). The AMI of the AD patients decreased significantly more slowly with time delays than the AMI of normal controls at electrodes T5, T6, O1, O2, P3 and P4 (p<0.01). Furthermore, we observed a strong correlation between the results obtained with both non-linear methods, suggesting that the AMI rate of decrease can be used to estimate the regularity in time series. The decreased irregularity found in AD patients suggests that EEG analysis with ApEn and AMI could help to increase our insight into brain dysfunction in AD.
Abásolo D, Gómez N, Santamarta D, García M, Hornero R (2012) Characterisation of Intracranial Pressure in Infusion Studies with Sample Entropy, Proceedings of the 5th International Conference on Medical Signal & Information Processing (MEDSIP 2012)
Neuronal outputs are complex signals of dynamically integrated excitatory and inhibitory components. Decreased synaptic inhibition in a neuronal network increases excitability and multiple spiking in neurons. Synchronized multiple spiking among a neuronal population further generates rhythmic field potentials and this epileptiform activity can propagate in the brain and cause seizures. Pharmacological interventions that reduce rhythmicity of epileptiform activity may have antiepileptic potentials. We evaluated the Lempel-Ziv (LZ) complexity for identifying rhythmicity in population spikes recorded in granule cells of the murine dentate gyrus in vitro. Blocking synaptic inhibition by the GABAA receptor antagonist, bicuculline, caused epileptiform population spikes, and we found that the LZ complexity of the signal was significantly reduced. Moreover, the GABAB receptor agonist, baclofen, reduced the amplitude of the epileptiform population spike and we found that it increased LZ complexity. The results show that LZ complexity is sensitive to pharmacological interventions that apparently alter rhythmicity of neuronal outputs by desynchronizing neuronal population firing. This novel approach in neuronal signal processing may be used to identify new antiepileptic targets. © Springer International Publishing Switzerland 2014.
The advent of new signal processing methods, such as non-linear analysis techniques, represents a new perspective which adds further value to brain signals' analysis. Particularly, Lempel-Ziv's Complexity (LZC) has proven to be useful in exploring the complexity of the brain electromagnetic activity. However, an important problem is the lack of knowledge about the physiological determinants of these measures. Although a correlation between complexity and connectivity has been proposed, this hypothesis was never tested in vivo. Thus, the correlation between the microstructure of the anatomic connectivity and the functional complexity of the brain needs to be inspected. In this study we analyzed the correlation between LZC and fractional anisotropy (FA), a scalar quantity derived from diffusion tensors that is particularly useful as an estimate of the functional integrity of myelinated axonal fibers, in a group of sixteen healthy adults (all female, mean age 65.56±6.06 years, intervals 58-82). Our results showed a positive correlation between FA and LZC scores in regions including clusters in the splenium of the corpus callosum, cingulum, parahipocampal regions and the sagittal stratum. This study supports the notion of a positive correlation between the functional complexity of the brain and the microstructure of its anatomical connectivity. Our investigation proved that a combination of neuroanatomical and neurophysiological techniques may shed some light on the underlying physiological determinants of brain's oscillations.
The Surrey Virtual Rehabilitation System (SVRS) was developed for clinical use in physiotherapy for children with cerebral palsy (CP). The overall aim of this work was to evaluate the practicality of the SVRS for three simple lower extremity exercises. Two children with CP, two clinical engineering trainees and a physiotherapist performed the tests whilst parents and a paediatric physiotherapist observed the exercise sessions. The feedback collected during an open-discussion and the descriptive analysis of responses to 15 closed-ended questions suggest that the participants were satisfied with the practicality of the SVRS. Outcome measures derived from data collected during the sessions indicate that the SVRS may provide clinically relevant feedback on the performance of patients for themselves and their treating clinicians. In conclusion, the SVRS appears practical for rehabilitation purposes and is worthy of further evaluation and development. Copyright © 2014 Inderscience Enterprises Ltd.
The aim of this study was to explore the ability of several spectral entropies to discriminate between spontaneous magnetoencephalographic (MEG) oscillations from 20 Alzheimer's disease (AD) patients and 21 controls. Hence, the relative spectral power (RSP) in classical frequency bands was calculated from the averaged power spectral density. Given the fact that the RSP can be viewed as a probability distribution function, the Shannon spectral entropy, Tsallis spectral entropy, generalized escort-Tsallis spectral entropy and Rényi spectral entropy were calculated from the RSP. Significant differences for each parameter were assessed with Mann-Whitney U test, whereas classification performance was studied using binary logistic regression. Results revealed an increase in the RSP of control subjects at beta and gamma bands, while AD patients showed an increase in the RSP values at delta and theta bands. Entropies obtained statistically significant lower values for AD patients than for controls. This issue suggests a significant decrease in irregularity of AD patients' MEG activity.
Alzheimer's disease (AD) is one of the most prominent neurodegenerative disorders. The aim of this research work is to study the magnetoencephalogram (MEG) background activity in AD patients using two scaling analysis methods: detrended fluctuation analysis (DFA) and backward detrended moving average (BDMA). Both measures have been designed to quantify correlations in noisy and non-stationary signals. Five minutes of recording were acquired with a 148-channel whole-head magnetometer in 15 patients with probable AD and 15 control subjects. Both DFA and BDMA exhibited two scaling regions with different slopes. Significant differences between both groups were found in the second region of DFA and in the first region of BDMA (p < 0.01, Student's t-test). Using receiver operating characteristic curves, accuracies of 83.33% with DFA and of 80% with BDMA were reached. Our findings show the usefulness of these scaling analysis methods to increase our insight into AD.
Alzheimer's disease (AD) is the main cause of dementia in western countries. Although a definite diagnosis of this illness is only possible by necropsy, the analysis of nonlinear dynamics in electroencephalogram (EEG) signals could help physicians in this difficult task In this study we have applied Approximate Entropy (ApEn) to analyze the EEG background activity of patients with a clinical diagnosis of Alzheimer's disease and control subjects. ApEn is a newly introduced statistic that can be used to quantify the complexity (or irregularity) of a time series. We have divided the EEG data into frames to calculate their ApEn. Our results show that the degree of complexity of EEGs from control subjects is higher. Applying the ANOVA test, we have verified that there was a significant difference (p < 0.05) between the EEGs of these groups.
The aim of this study was to analyse the electroencephalogram (EEG) background activity of Alzheimer's disease (AD) patients using multiscale entropy (MSE). MSE is a recently developed method that quantifies the regularity of a signal on different time scales. These time scales are inspected by means of several coarse-grained sequences formed from the analysed signals. We recorded the EEGs from 19 scalp electrodes in 11 AD patients and 11 age-matched controls and estimated the MSE profile for each epoch of the EEG recordings. The shape of the MSE profiles reveals the EEG complexity, and it suggests that the EEG contains information in deeper scales than the smallest one. Moreover, the results showed that the EEG background activity is less complex in AD patients than control subjects. We found significant differences between both subject groups at electrodes F3, F7, Fp1, Fp2, T5, T6, P3, P4, O1 and O2 (p-value < 0.01, Student's t-test). These findings indicate that the EEG complexity analysis performed on deeper time scales by MSE may be a useful tool in order to increase our knowledge of AD.
Epileptic seizures are generated by an abnormal synchronization of neurons unforeseeable for the patients. In this study we analyzed invasive electroencephalogram (EEG) recordings in patients suffering from medically intractable focal epilepsy with two non-linear methods, Approximate Entropy (ApEn) and Lempel-Ziv (LZ) complexity. ApEn and LZ complexity quantify the regularity and complexity of a time series, respectively, and are well suited to the analysis of non-stationary biomedical signals of short length. Our results show an increase in ApEn and LZ complexity values during seizures at the focal electrodes. These changes could also be seen at some extra focal electrodes. After the seizure ends, the values of both non-linear metrics return to values lower than those before the seizure. Moreover, we quantified the changes in LZ complexity, showing the complexity increase during the seizure and its notable decrease after its end. Our results suggest that these techniques are useful to detect changes due to epileptic seizures in the EEG.
Analysis of electrophysiological signals recorded from the brain with Lempel-Ziv (LZ) complexity, a measure based on coarse-graining of the signal, can provide valuable insights into understanding brain activity. LZ complexity of local field potential signals recorded from the neocortex of 11 adult male Wistar-Kyoto rats in different vigilance states - waking, non-rapid-eye movement (NREM) and REM sleep - was estimated with different coarse-graining techniques (median, LZCm, and k-means, LZCkm). Furthermore, surrogate data were used to test the hypothesis that LZ complexity results reveal effects accounted for by temporal structure of the signal, rather than merely its frequency content. LZ complexity values were significantly lower in NREM sleep as compared to waking and REM sleep, for both real and surrogate signals. LZCkm and LZCm values were similar, although in NREM sleep the values deviated in some epochs, where signals also differed significantly in terms of temporal structure and spectral content. Thus, the interpretation of LZ complexity results should take into account the specific algorithm used to coarse-grain the signal. Moreover, the occurrence of high amplitude slow waves during NREM sleep determines LZ complexity to a large extent, but characteristics such as the temporal sequence of slow waves or cross-frequency interactions might also play a role. © Springer International Publishing Switzerland 2014.
In this proof-of-principle study we analyzed intracranial electroencephalogram recordings in patients with intractable focal epilepsy. We contrast two implementations of Independent Component Analysis (ICA) - Ensemble (or spatial) ICA (E-ICA) and Space-Time ICA (ST-ICA) in separating out the ictal components underlying the measurements. In each case we assess the outputs of the ICA algorithms by means of a non-linear method known as the Lempel-Ziv (LZ) complexity. LZ complexity quantifies the complexity of a time series and is well suited to the analysis of non-stationary biomedical signals of short length. Our results show that for small numbers of intracranial recordings, standard E-ICA results in marginal improvements in the separation as measured by the LZ complexity changes. ST-ICA using just 2 recording channels both near and far from the epileptic focus result in more distinct ictal components--although at this stage there is a subjective element to the separation process for ST-ICA. Our results are promising showing that it is possible to extract meaningful information from just 2 recording electrodes through ST-ICA, even if they are not directly over the seizure focus. This work is being further expanded for seizure onset analysis.
Hornero R, Alvarez D, Abásolo D, Gómez C, del Campo F, Zamarrón C (2005) Approximate entropy from overnight pulse oximetry for the obstructive sleep apnea syndrome, 2005 27th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vols 1-7pp. 6157-6160 IEEE
The aim of this preliminary study was to asses whether a time series analysis from overnight pulse oximetry by means of Approximate Entropy (ApEn) could yield essential information on the diagnosis of obstructive sleep apnea (OSA) syndrome. We analyzed the oxygen saturation (SaO(2)) signals of 74 patients (44 with a positive diagnosis of OSA and 30 with a negative diagnosis of OSA) by means of ApEn, which quantified the regularity (or complexity) of time series. Results showed that recurrence of apnea events in patients with OSA determined a significant increase in ApEn values with a mean standard deviation (SD) of 1.07 +/- 0.30. A mean +/- SD ApEn value of 0.47 +/- 0.25 was estimated in patients without OSA. We obtained an area under the ROC curve of 0.94. The optimum threshold was selected at 0.81, where we achieved a 79,5% sensitivity and 90% specificity. Further analyses are necessary with new and larger data set to test the potential value of our methodology prospectively.
Currently accepted input parameter limitations in entropy-based, non-linear signal processing methods, for example, sample entropy (SampEn), may limit the information gathered from tested biological signals. The ability of quadratic sample entropy (QSE) to identify changes in electroencephalogram (EEG) signals of 11 patients with a diagnosis of Alzheimer's disease (AD) and 11 age-matched, healthy controls is investigated. QSE measures signal regularity, where reduced QSE values indicate greater regularity. The presented method allows a greater range of QSE input parameters to produce reliable results than SampEn. QSE was lower in AD patients compared with controls with significant differences (p < 0.01) for different parameter combinations at electrodes P3, P4, O1 and O2. Subject- and epoch-based classifications were tested with leave-one-out linear discriminant analysis. The maximum diagnostic accuracy and area under the receiver operating characteristic curve were 77.27 and more than 80%, respectively, at many parameter and electrode combinations. Furthermore, QSE results across all r values were consistent, suggesting QSE is robust for a wider range of input parameters than SampEn. The best results were obtained with input parameters outside the acceptable range for SampEn, and can identify EEG changes between AD patients and controls. However, caution should be applied because of the small sample size.
Simons S, Abásolo D, Hughes M (2013) Distance Measure Based on Lempel-Ziv Complexity for Analysis of Electroencephalograms in Alzheimer's Disease, Proceedings of the 7th United Kingdom & Republic of Ireland Postgraduate Conference in Biomedical Engineering and Medical Physics (PGBiomed 2013) University of Surrey
de la Rosa R, Alonso A, de la Rosa S, Abásolo D (2008) Myo-Pong: a neuromuscular game for the UVa-Neuromuscular Training System Platform, 2008 VIRTUAL REHABILITATIONpp. 61-61 IEEE
This paper presents the game Myo-Pong, a simple graphical table-tennis game included in the portable UVa Neuromuscular Training System (UVa-NTS). Myo-Pong demonstrates the capabilities of the UVa-NTS as a myoelectric real-time system for training and for playing by means of myoelectric control.
Non-linear analysis of the electroencephalogram (EEG) background activity can help to obtain a better understanding of abnormal dynamics in the brain. The aim of this study was to analyze the regularity of the EEG time series of Alzheimer's disease (AD) Patients to test the hypothesis that the irregularity of the AD patients' EEG is lower than that of age-matched controls. We recorded the EEG from 19 scalp electrodes in 11 AD patients and 11 age-matched controls and estimated the Approximate Entropy (ApEn). ApEn is a non-linear method that can be used to quantify the irregularity of a time series. Larger values correspond to more irregularity. We evaluated different values for input parameters m and r to estimate ApEn and concluded that m=1 and r=0.25 times the SD of the time series were the optimum choices. With these parameters, ApEn was significantly lower in the AD patients at the P3, P41 O1 and O2 (p < 0.01) electrodes. The decreased irregularity found in the EEG of AD patients in the parietal and occipital regions leads us to think that regularity analysis of the EEG with ApEn could be a useful tool to increase our insight into brain dysfunction in Alzheimer's disease.
Alzheimer's disease (AD) is the most frequent form of dementia in western countries. An early detection would be beneficial, but currently diagnostic accuracy is relatively poor. In this study, differences in information content between cortical areas in 12 AD patients and 11 control subjects were assessed with Kullback-Leibler (KL) entropy. KL entropy measures the degree of similarity between two probability distributions. EEGs were recorded from 19 scalp electrodes and KL entropy values of the EEGs in both groups were estimated for the local, distant and interhemispheric electrodes. KL entropy values were lower in AD patients than in age-matched control subjects, with significant effects for diagnosis and brain region (p < 0.05, two-way ANOVA). No significant interaction for diagnosis X region was found (p = 0.7671). Additionally a one-way ANOVA showed that KL entropy values were significantly lower in AD patients (p < 0.05) for the distant electrodes on the right hemisphere. These results suggest that KL entropy highlights information content changes in the EEG due to AD. However, further studies are needed to address the possible usefulness of KL entropy in the characterisation and early detection of AD.
The aim of the present study was to analyse the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD) using the Lempel-Ziv (LZ) complexity. This non-linear method measures the complexity of finite sequences and is related to the number of distinct substrings and the rate of their occurrence along the sequence. The MEGs were recorded with a 148-channel whole-head magnetometer (MAGNES 2500 WH, 4D Neuroimaging) in 21 patients with AD and in 21 age-matched control subjects. Artefact-free epochs were selected for complexity analysis. Results showed that MEG signals from AD patients had lower complexity than control subjects' MEGs and the differences were statistically significant (p<0.01). In order to reduce the dimension of the LZ complexity results, a principal components analysis (PCA) was applied, and only the first principal component was retained. The first component score from PCA was graphically analysed using a box plot and a receiver-operating characteristic (ROC) curve. A specificity of 85.71%, a sensitivity of 80.95% and an area under the ROC curve of 0.9002 were obtained. These preliminary results suggest that cognitive dysfunction in AD is associated with a decreased LZ complexity in the MEG signals.
Sample entropy (SampEn) is a nonlinear regularity index that requires the a priori selection of three parameters: the length of the sequences to be compared, m, the patterns similarity tolerance, r, and the number of samples under analysis, N. Appropriate values for m, r and N have been recommended and widely used in the literature for the application of SampEn to some physiological time series, such as heart rate, hormonal data, etc. However, no guidelines exist for the selection of that values in other cases. Therefore, an optimal parameters study should be required for the application of SampEn to not previously analyzed biomedical signals. In the present work, a thorough analysis on the optimal values for m, r and N is presented within the context of atrial fibrillation (AF) organization estimation, computed from surface electrocardiogram recordings. Recently, the evaluation of AF organization through SampEn, has revealed clinically useful information that could be used for a better treatment of this arrhythmia. The present study analyzed optimal SampEn parameter values within two different scenarios of AF organization estimation, such as the prediction of paroxysmal AF termination and the electrical cardioversion outcome in persistent AF. As a result, interesting recommendations about the selection of m, r and N, together with the relationship between N and the sampling rate (f(s)) were obtained. More precisely, (i) the proportion between N and f(s) should be higher than 1s and f(s) >= 256 Hz, (ii) overlapping between adjacent N-length windows does not improve AF organization estimation with respect to the analysis of non-overlapping windows, and (iii) values of m and r maximizing successful classification for the analyzed AF databases should be considered within a range wider than the proposed in the literature for heart rate analysis, i.e. m = 1 and m = 2 and r between 0.1 and 0.25 times the standard deviation of the data. (C) 2010 Elsevier Ireland Ltd. All rights reserved.
García M, Poza J, Abásolo D, Santamarta D, Hornero R (2013) Analysis of intracranial pressure signals recorded during infusion studies using the spectral entropy., Conf Proc IEEE Eng Med Biol Soc2013pp. 2543-2546
Hydrocephalus includes a range of disorders characterized by clinical symptoms, abnormal brain imaging and altered cerebrospinal fluid (CSF) dynamics. Infusion tests can be used to study CSF circulation in patients with hydrocephalus. In them, intracranial pressure (ICP) is deliberately raised and CSF circulation disorders evaluated through measurements of the resulting ICP. In this study, we analyzed 77 ICP signals recorded during infusion tests using the spectral entropy (SE). Each signal was divided into four artifact-free epochs. The mean SE, , and the standard deviation of SE, SD[SE], were calculated for each epoch. Statistically significant differences were found between phases of the infusion test using and SD[SE] (p<1.7 · 10(-3), Bonferroni-corrected Wilcoxon tests). Furthermore, we found significantly lower and SD[SE] values in the plateau phase than in the basal phase. These findings suggest that the increase in ICP during infusion studies is associated with a significant decrease in irregularity and variability of the spectral content of ICP signals, measured in terms of SE. We conclude that the spectral analysis of ICP signals could be useful for understanding CSF dynamics in hydrocephalus.
Abásolo D, Simons S (2012) Fuzzy Entropy analysis of the electroencephalogram background activity in Alzheimer's disease, BIOMAG 2012 Book of Abstracts - 18th International Conference on Biomagnetism
Hornero R, Aboy M, Abásolo D, McNames J, Wakeland W, Goldstein B (2006) Complex analysis of intracranial hypertension using approximate entropy., Crit Care Med34(1)pp. 87-95
OBJECTIVE: To determine whether decomplexification of intracranial pressure dynamics occurs during periods of severe intracranial hypertension (intracranial pressure >25 mm Hg for >5 mins in the absence of external noxious stimuli) in pediatric patients with intracranial hypertension. DESIGN: Retrospective analysis of clinical case series over a 30-month period from April 2000 through January 2003. SETTING: Multidisciplinary 16-bed pediatric intensive care unit. PATIENTS: Eleven episodes of intracranial hypertension from seven patients requiring ventriculostomy catheter for intracranial pressure monitoring and/or cerebral spinal fluid drainage. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: We measured changes in the intracranial pressure complexity, estimated by the approximate entropy (ApEn), as patients progressed from a state of normal intracranial pressure (<25 mm Hg) to intracranial hypertension. We found the ApEn mean to be lower during the intracranial hypertension period than during the stable and recovering periods in all the 11 episodes (0.5158 +/- 0.0089, 0.3887 +/- 0.077, and 0.5096 +/- 0.0158, respectively, p < .01). Both the mean reduction in ApEn from the state of normal intracranial pressure (stable region) to intracranial hypertension (-0.1271) and the increase in ApEn from the ICH region to the recovering region (0.1209) were determined to be statistically significant (p < .01). CONCLUSIONS: Our results indicate that decreased complexity of intracranial pressure coincides with periods of intracranial hypertension in brain injury. This suggests that the complex regulatory mechanisms that govern intracranial pressure may be disrupted during acute periods of intracranial hypertension. This phenomenon of decomplexification of physiologic dynamics may have important clinical implications for intracranial pressure management.
Sample Entropy (SampEn) is a nonlinear regularity index that requires the a priori selection of three parameters: the length of the sequences to be compared, m, the patterns similarity tolerance, r, and the number of samples under analysis, N. Appropriate values for m, r and N have been recommended in some cases, such as heart rate, hormonal data, etc., but no guidelines exist for the selection of that values. Hence, an optimal parameters study should be required for the application of SampEn to not previously analyzed biomedical signals. In this work, a thorough analysis on the optimal SampEn parameter values within two different scenarios of AF organization estimation, such as the prediction of paroxysmal AF termination and the electrical cardioversion outcome in persistent AF, is presented. Results indicated that, (i) the proportion between N and the sampling rate (fs) should be higher than one second and fs e 256 Hz, (ii) overlapping between adjacent N-length windows does not improve organization estimation and (iii) values of m and r maximizing classification should be considered within a range wider than the proposed in the literature for heart rate analysis.
A Real-time Treadmill Speed Control Algorithm (RTSCA) has been developed for gait
rehabilitation of children with cerebral palsy (CP). The objective of the work described in this
paper was to investigate the feasibility of the RTSCA prior to use by children with CP.
Thirteen healthy subjects aged between 19 and 25 were recruited to walk on the treadmill using
conventional speed buttons without the virtual reality (VR) environment, and the RTSCA with
and without VR. The participants were asked to undertake three treadmill tests and to complete
a questionnaire to provide feedback on the control of the treadmill. The descriptive results
show that for 10 participants changing walking speed from stationary when using the RTSCA
was similar or more comfortable to using conventional treadmill speed control buttons. For
those who found it less comfortable the core issue was insufficient time to practise with the
system. All the participants were satisfied with the safety and the performance of the RTSCA
when incorporated into the VR scenario. A Wilcoxon test was conducted to examine whether
there was a significant difference between walking speeds on the treadmill when using the
conventional speed buttons and the RTSCA. The results showed that participants walked at
significantly higher speeds when using the RTSCA. This may suggest that they walked more
naturally or confidently on the treadmill when using the RTSCA as compared to the use of
conventional treadmill speed control buttons.
Nocturnal oximetry is an attractive option for the diagnosis of obstructive sleep apnoea (OSA) syndrome because of its simplicity and low cost compared to polysomnography (PSG). The present study assesses nonlinear analysis of blood oxygen saturation (SaO(2)) from nocturnal oximetry as a diagnostic test to discriminate between OSA positive and OSA negative patients. A sample of 187 referred outpatients, clinically suspected of having OSA, was studied using nocturnal oximetry performed simultaneously with complete PSG. A positive OSA diagnosis was found for 111 cases, while the remaining 76 cases were classified as OSA negative. The following oximetric indices were obtained: cumulative time spent below a saturation of 90% (CT90), oxygen desaturation indices of 4% (ODI4), 3% (ODI3) and 2% (ODI2) and the delta index (Delta index). SaO(2) records were subsequently processed applying two nonlinear methods: central tendency measure (CTM) and Lempel-Ziv (LZ) complexity. Significant differences (p < 0.01) were found between OSA positive and OSA negative patients. Using CTM we obtained a sensitivity of 90.1% and a specificity of 82.9%, while with LZ the sensitivity was 86.5% and the specificity was 77.6%. CTM and LZ accuracies were higher than those provided by ODI4, ODI3, ODI2 and CT90. The results suggest that nonlinear analysis of SaO(2) signals from nocturnal oximetry could yield useful information in OSA diagnosis.
Background: The complexity of the intracranial pressure (ICP) signal decreases with intracranial hypertension in children with acute brain injury as well as during infusion studies in adults with hydrocephalus. In this study we have analysed the pressure signal obtained in the lumbar subarachnoid space during infusion testing. The pulse amplitude rises when the ICP is increased by additional external volume. Our objective was to determine the relative influence of the pressure range and pulse amplitude on the loss of complexity observed during infusion-related intracranial hypertension. Materials and Methods: The Lempel-Ziv (LZ) complexity of the cerebrospinal fluid pressure (CSFP) signal was analysed in 52 infusion studies performed in patients with normal pressure hydrocephalus (median age 71 years, IQR: 60-78). Four sequences during the baseline, infusion, steady plateau and recovery periods of each infusion study were selected. The mean values of the CSFP (mCSFP), pulse amplitude and LZ complexity in every sequence were measured. Correlations between LZ complexity and CSFP parameters were explored. Results: Significant inverse correlations were found among LZ complexity, pulse amplitude and mCSFP during all periods of infusion testing, except at baseline. Partial correlation analysis controlling the effect of mCSFP emphasised the relationship between pulse amplitude and LZ complexity. When pulse amplitude is held constant the partial correlation between LZ complexity and mCSFP is not significant. Conclusions: The pulse amplitude of the CSFP signal seems to be a major determinant of the waveform complexity. © 2012 Springer-Verlag/Wien.
Simons S, Abásolo D (2014) Can Distance Measures Based on Lempel-Ziv Complexity Help in the Detection of Alzheimer's Disease from Electroencephalograms?, IFMBE Proceedings41pp. 698-701
This pilot study applied three distance-based bivariate Lempel-Ziv complexity (LZC) measures to investigate the changes in electroencephalogram (EEG) signals between 11 patients with Alzheimer's disease (AD) and 11 age matched controls. These methods measure richness of complexity between pairs of signals. Complexity of control subjects' EEGs was richer, i.e. signals were made from a greater number and greater range of subsequences, than those from AD patients in almost all cases in two non-normalized distance- based methods. Only some pairs including electrode T4 (2.1% of the total) occasionally showed the reverse result. Statistically significant differences were found with these two methods in 21 and 18 of 120 tested electrode pairs, respectively (Student's t test, p<0.01). Receiver operating curves were used to calculate the sensitivity (number of correctly classified AD patients) and specificity (number of correctly classified controls). Accuracy is the combined correct classification of controls and AD patients. The maximum sensitivity found was 100%, specificity 90.9% and accuracy 86.4% at various electrode pairs with both non-normalized methods. The normalized method showed many electrode pairs with increased richness of complexity for AD patients than controls (67.5% of the total). It was found that this was due to the normalization procedure modifying the distribution of the original complexities from the electrode pairs. These findings suggest non-normalized distance-based bivariate LZC measures can be reliably applied to complex physiological signals such as human EEGs to further understand the effect of AD on the complexity of brain signals of patients. However, care must be taken when normalization procedures are applied. © Springer International Publishing Switzerland 2014.
Sample Entropy (SampEn) is a nonlinear regularity index that requires the a priori selection of three parameters: the length of the sequences to be compared, m, the patterns similarity tolerance, r, and the number of samples under analysis, N. Appropriate values for m, r and N have been recommended in some cases, such as heart rate, hormonal data, etc., but no guidelines exist for the selection of that values. Hence, an optimal parameters study should be required for the application of SampEn to not previously analyzed biomedical signals. In this work, a thorough analysis on the optimal SampEn parameter values within two different scenarios of AF organization estimation, such as the prediction of paroxysmal AF termination and the electrical cardioversion outcome in persistent AF, is presented. Results indicated that, (i) the proportion between N and the sampling rate ((s)) should be higher than one second and (s) e 256 Hz, (ii) overlapping between adjacent N-length windows does not improve organization estimation and (iii) values of m and r maximizing classification should be considered within a range wider than the proposed in the literature for heart rate analysis, i. e. m = 1 and m = 2 and r between 0.1 and 0.25 times the standard deviation of the data.
OBJECTIVE: The aim of this study was to analyse the regularity of the EEG background activity of Alzheimer's disease (AD) patients to test the hypothesis that the irregularity of the AD patients' EEG is lower than that of age-matched controls. METHODS: We recorded the EEG from 19 scalp electrodes in 10 AD patients and 8 age-matched controls and estimated the Approximate Entropy (ApEn). ApEn is a non-linear statistic that can be used to quantify the irregularity of a time series. Larger values correspond to more complexity or irregularity. A spectral analysis was also performed. RESULTS: ApEn was significantly lower in the AD patients at electrodes P3 and P4 (P < 0.01), indicating a decrease of irregularity. We obtained 70% sensitivity and 100% specificity at P3, and 80% sensitivity and 75% specificity at P4. Results seemed to be complementary to spectral analysis. CONCLUSIONS: The decreased irregularity found in the EEG of AD patients in the parietal region leads us to think that EEG analysis with ApEn could be a useful tool to increase our insight into brain dysfunction in AD. However, caution should be applied due to the small sample size. SIGNIFICANCE: This article represents a first step in demonstrating the feasibility of ApEn for recognition of EEG changes in AD.
PURPOSE: Nonlinear dynamics has enhanced the diagnostic abilities of some physiological signals. Recent studies have shown that the complexity of the intracranial pressure waveform decreases during periods of intracranial hypertension in paediatric patients with acute brain injury. We wanted to assess changes in the complexity of the cerebrospinal fluid (CSF) pressure signal over the large range covered during the study of CSF circulation with infusion studies. METHODS: We performed 37 infusion studies in patients with hydrocephalus of various types and origin (median age 71 years; interquartile range 60-77 years). After 5 min of baseline measurement, infusion was started at a rate of 1.5 ml/min until a plateau was reached. Once the infusion finished, CSF pressure was recorded until it returned to baseline. We analysed CSF pressure signals using the Lempel-Ziv (LZ) complexity measure. To characterise more accurately the behaviour of LZ complexity, the study was segmented into four periods: basal, early infusion, plateau and recovery. RESULTS: The LZ complexity of the CSF pressure decreased in the plateau of the infusion study compared to the basal complexity (p=0.0018). This indicates loss of complexity of the CSF pulse waveform with intracranial hypertension. We also noted that the level of complexity begins to increase when the infusion is interrupted and CSF pressure drops towards the initial values. CONCLUSIONS: The LZ complexity decreases when CSF pressure reaches the range of intracranial hypertension during infusion studies. This finding provides further evidence of a phenomenon of decomplexification in the pulsatile component of the pressure signal during intracranial hypertension.
Hornero R, Abásolo DE, Espino P (2003) Use of wavelet entropy to compare the EEG background activity of epileptic patients and control subjects, Proceedings - 7th International Symposium on Signal Processing and Its Applications, ISSPA 20032pp. 5-8
In this study, we processed the EEG background activities of five epileptic patients by means of the wavelet packet transform (WPT), and compared them with the EEGs of five control subjects. The WPT provides a fast way of choosing the best out of a large number of time-frequency bases according to features of the signal to be analysed. From the wavelet packet coefficients, we calculated an information cost function (ICF), which evaluates the information associated with the energy distribution of the signal: the wavelet entropy. Results suggested that the entropy or degree of disorder in the energy distribution of control subjects' EEG was higher. Applying the ANOVA test, we verified that there was a significant difference (p-ANOV A<0.01 ) between EEG frames of both groups. © 2003 IEEE.
The goal of this study was to analyze the magnetoencephalogram (MEG) background activity in patients with Alzheimer's disease (AD) using a nonlinear forecasting measure. It is a nonparametric method to quantify the predictability of time series. Five minutes of recording were acquired with a 148-channel whole-head magnetometer in 15 patients with probable AD and 15 elderly control subjects. Stationary epochs of 5 seconds (848 points, sample frequency of 169.55 Hz) were selected. Our results showed that AD patients' MEGs were more predictable than controls' recordings. Additionally, an accuracy of 76.7% (80.0% sensitivity; 73.3% specificity) was reached using a receiver operating characteristic curve. These preliminary results suggest the usefulness of nonlinear forecasting to gain a better understanding of dynamical processes underlying the MEG recording.
The analysis of P-wave variability from the electrocardiogram (ECG) has been suggested as an early predictor of the onset of paroxysmal atrial fibrillation (PAF). Hence, a preventive treatment could be used to avoid the loss of normal sinus rhythm, thus minimising health risks and improving the patient's quality of life. In these previous studies the variability of different temporal and morphological P-wave features has been only analysed in a linear fashion. However, the electrophysiological alteration occurring in the atria before the onset of PAF has to be considered as an inherently complex, chaotic and non-stationary process. This work analyses the presence of non-linear dynamics in the P-wave progression before the onset of PAF through the application of the central tendency measure (CTM), which is a non-linear metric summarising the degree of variability in a time series. Two hour-length ECG intervals just before the arrhythmia onset belonging to 46 different PAF patients were analysed. In agreement with the invasively observed inhomogeneous atrial conduction preceding the onset of PAF, CTM for all the considered P-wave features showed higher variability when the arrhythmia was closer to its onset. A diagnostic accuracy around 80% to discern between ECG segments far from PAF and close to PAF was obtained with the CTM of the metrics considered. This result was similar to previous P-wave variability methods based on linear approaches. However, the combination of linear and non-linear methods with a decision tree improved considerably their discriminant ability up to 90%, thus suggesting that both dynamics could coexist at the same time in the fragmented depolarisation of the atria preceding the arrhythmia.
The Cluster-Span Threshold (CST) is a recently introduced unbiased threshold for functional connectivity networks. This binarisation technique offers a natural trade-off of sparsity and density of information by balancing the ratio of closed to open triples in the network topology. Here we present findings comparing it with the Union of Shortest Paths (USP), another recently proposed objective method. We analyse standard network metrics of binarised networks for sensitivity to clinical Alzheimer?s disease in the Beta band of Electroencephalogram activity. We find that the CST outperforms the USP, as well as subjective thresholds based on fixing the network density, as a sensitive threshold for distinguishing differences in the functional connectivity between Alzheimer?s disease patients and control. This study provides the first evidence of the usefulness of the CST for clinical research purposes.
Functional connectivity has proven useful to characterise electroencephalogram (EEG) activity in Alzheimer?s disease (AD). However, most current functional connectivity analyses have been static, disregarding any potential variability of the connectivity with time. In this pilot study, we compute short-time resting state EEG functional connectivity based on the imaginary part of coherency for 12 AD patients and 11 controls. We derive binary unweighted graphs using the cluster-span threshold, an objective binary threshold. For each short-time binary graph, we calculate its local clustering coefficient (Cloc), degree (K), and efficiency (E). The distribution of these graph metrics for each participant is then characterised with four statistical moments: mean, variance, skewness, and kurtosis. The results show significant differences between groups in the mean of K and E, and the kurtosis of Cloc and K. Although not significant when considered alone, the skewness of Cloc is the most frequently selected feature for the discrimination of subject groups. These results suggest that the variability of EEG functional connectivity may convey useful information about AD.
Alzheimer?s disease (AD) is a degenerative brain disorder leading to memory loss and changes in other cognitive abilities. The complexity of electroencephalogram (EEG) signals may help to characterise AD. To this end, we propose an extension of multiscale entropy based on variance (MSEÃ2) to multichannel signals, termed multivariate MSEÃ2 (mvMSEÃ2), to take into account both the spatial and time domains of time series. Then, we investigate the mvMSEÃ2 of EEGs at different frequency bands, including the broadband signals filtered between 1 and 40Hz, ¸, ±, and ² bands, and compare it with the previously proposed multiscale entropy based on mean (MSEµ), multivariate MSEµ (mvMSEµ), and MSEÃ2, to distinguish different kinds of dynamical properties of spread and mean in the signals. Results from 11 AD patients and 11 age-matched controls suggest that the presence of broadband activity of EEGs is required for a proper evaluation of complexity. MSEÃ2 and mvMSEÃ2 results, showing a loss of complexity in AD signals, led to smaller p-values in comparison with MSEµ and mvMSEµ ones, suggesting that the variance-based MSE and mvMSE can characterise changes in EEGs as a result of AD in a more detailed way. The p-values for the slope values of mvMSE curves were smaller than for MSE at large scale factors, also showing the possible usefulness of multivariate techniques.
We propose a novel complexity measure to overcome the deficiencies of the widespread and powerful multiscale entropy (MSE), including, MSE values may be undefined for short signals, and MSE is slow for real-time applications. Methods:
We introduce multiscale dispersion entropy (DisEn - MDE) as a very fast and powerful method to quantify the complexity of signals. MDE is based on our recently developed DisEn, which has a computation cost of O(N), compared with O(N2) for sample entropy used in MSE. We also propose the refined composite MDE (RCMDE) to improve the stability of MDE. Results:
We evaluate MDE, RCMDE, and refined composite MSE (RCMSE) on synthetic signals and three biomedical datasets. The MDE, RCMDE, and RCMSE methods show similar results, although the MDE and RCMDE are faster, lead to more stable results, and discriminate different types of physiological signals better than MSE and RCMSE. Conclusion:
For noisy short and long time series, MDE and RCMDE are noticeably more stable than MSE and RCMSE, respectively. For short signals, MDE and RCMDE, unlike MSE and RCMSE, do not lead to undefined values. The proposed MDE and RCMDE are significantly faster than MSE and RCMSE, especially for long signals, and lead to larger differences between physiological conditions known to alter the complexity of the physiological recordings.
The analysis of electroencephalograms (EEGs) of patients with Alzheimer?s disease (AD) could contribute to the diagnosis of this dementia. In this study, a new non-linear signal processing metric, distance-based Lempel-Ziv complexity (dLZC), is introduced to characterise changes between pairs of electrodes in EEGs in AD. When complexity in each signal arises from different sub-sequences dLZC would be greater than when similar sub-sequences are present in each signal. EEGs from 11 AD patients and 11 age-matched control subjects were analysed. The dLZC values for AD patients were lower than for control subjects for most electrode pairs, with statistically significant differences (p<0.01, Student?s t-test) in 17 electrode pairs in the distant left, local posterior left, and interhemispheric regions. Maximum diagnostic accuracies with leave-one-out cross-validation were 77.27% for subject-based classification and 78.25% for epoch-based classification. These findings suggest not only that EEGs from AD patients are less complex than those from controls, but also that the richness of the information contained in pairs of EEGs from patients is also lower than in age-matched controls. The analysis of EEGs in AD with dLZC may increase the insight into brain dysfunction providing complementary information to that obtained with other complexity and synchrony methods.
This study focuses on the resting state network analysis of the brain, as well as how these networks change both in topology and location throughout life. The magnetoencephalogram (MEG) background activity from 220 healthy volunteers (age 7-84 years), was analysed combining complex network analysis principles of graph theory with both linear and non-linear methods to evaluate the changes in the brain. Granger Causality (GC) (linear method) and Phase Slope Index (PSI) (non-linear method) were used to observe the connectivity in the brain during rest, and as a function of age by analysing the degree, clustering coefficient, efficiency, betweenness, modularity and maximised modularity of the observed complex brain networks. Our results showed that GC showed little linear causal activity in the brain at rest, with small world topology, while PSI showed little information flow in the brain, with random network topology. However, both analyses produced complementary results pertaining to the resting state of the brain.
The characterisation of healthy ageing of the brain could help create a fingerprint of normal ageing that might assist in the early diagnosis of neurodegenerative conditions. This study examined changes in resting state MEG permutation entropy due to age and gender in a sample of 220 healthy participants (98 males and 122 females, ages ranging between 7 and 84). Entropy was quantified using normalised permutation entropy and modified permutation entropy, with an embedding dimension of 5 and a lag of 1 as the input parameters for both algorithms. Effects of age were observed over the 5 regions of the brain i.e. anterior, central, posterior, and left and right lateral, with the anterior and central regions containing the highest permutation entropy. Statistically significant differences due to age were observed in the different brain regions for both genders, with the evolutions described using the fitting of polynomial regressions. Nevertheless, no significant differences between the genders were observed across all ages. These results suggest that the evolution of entropy in the background brain activity, quantified with permutation entropy algorithms, might be considered an alternative illustration of a ?nominal? physiological rhythm.
Hydrocephalus is a condition characterised by enlarged cerebral ventricles, which in turn affects intracranial pressure (ICP); however, the mechanisms regulating ICP are not fully understood. A non-linear signal processing approach was applied to ICP signals measured during infusion studies from patients with two forms of hydrocephalus, in a bid to compare the differences. This is the first study of its kind. The two forms of hydrocephalus were idiopathic normal pressure hydrocephalus (iNPH) and secondary normal pressure hydrocephalus (SH). Following infusion tests, the Lempel-Ziv (LZ) complexity was calculated from the iNPH and SH ICP signals. The LZ complexity values were averaged for the baseline, infusion, plateau and recovery stages of the tests. It was found that as the ICP increased from basal levels, the LZ complexities decreased, reaching their lowest during the plateau stage. However, the complexities computed from the SH ICP signals decreased to a lesser extent when compared to the iNPH ICP signals. Furthermore, statistically significant differences were found between the plateau and recovery stage complexities when comparing the iNPH and SH results (p << 0.05). This study suggests that advanced signal processing of ICP signals with LZ complexity can help characterise different types of hydrocephalus in more detail.
Alzheimer?s disease (AD) is neurodegenerative,
caused by the progressive death of brain cells over time. One
non-invasive approach to investigate AD is to use electroencephalogram
(EEG) signals. The data are usually non-stationary
with a strong background activity and noise which makes the
analysis difficult leading to low performance in many real
world applications including the detection of AD. In this study,
we present a method based on local texture changes of EEG
signals to differentiate AD patients from the healthy ones, using
one-dimensional local binary patterns (1D-LBPs) coupled with
support vector machines (SVM). Our proposed method maps
the EEG data into a less detailed representation which is less
sensitive to noise. A 10 fold cross validation performed at both
the epoch and subject level show the discriminancy power of
1D-LBP feature vectors with application to AD data.
Alzheimer?s disease (AD) is a neurodegenerative
disease caused by the progressive death of brain cells over
time. It represents the most frequent cause of dementia in the
western world, and affects an individual?s cognitive ability and
psychological capacity. While clinical diagnoses of AD are made
primarily on the basis of clinical evaluation and mental health
tests, diagnostic certainty is only possible through necropsy. One
non-invasive approach to investigating AD is to use electroencephalograms
(EEGs), which reflect brain electrical activity and
so can be used to detect electrical abnormalities in brain signals
with non-invasive cranial surface electrodes. Generally EEGs
in AD patients show a shift to lower frequencies in spectral
analysis and display less complexity and contain more regular
patterns compared to those of control subjects. Here we present
a method for differentiating AD patients from healthy ones based
on their EEG signals using Benford?s law and support vector
machines (SVMs) with a radial basis function (RBF) kernel.
EEG signals from eleven AD and eleven age-matched controls
were divided into artefact-free 5-sec epochs and used to train an
SVM. 10 fold cross validation was performed at both the epochand
subject-level to evaluate the importance of each electrode
in discriminating between AD and healthy subjects. Substantive
variability was seen across the different electrodes, with electrodes
O1, O2 and C4 particularly being important. Performance across
the electrodes was reduced when subject-level cross validation
was performed, but relative performance across the electrodes
was consistent with that found using epoch-level cross validation.
Research into binary network analysis of brain function faces a methodological challenge
in selecting an appropriate threshold to binarise edge weights. For EEG phase-based
functional connectivity, we test the hypothesis that such binarisation should take into
account the complex hierarchical structure found in functional connectivity. We explore
the density range suitable for such structure and provide a comparison of
state-of-the-art binarisation techniques, the recently proposed Cluster-Span Threshold
(CST), minimum spanning trees, efficiency-cost optimisation and union of shortest path
graphs, with arbitrary proportional thresholds and weighted networks. We test these
techniques on weighted complex hierarchy models by contrasting model realisations
with small parametric differences. We also test the robustness of these techniques to
random and targeted topological attacks.We find that the CST performs consistenty
well in state-of-the-art modelling of EEG network topology, robustness to topological
network attacks, and in three real datasets, agreeing with our hypothesis of hierarchical
complexity. This provides interesting new evidence into the relevance of considering a
large number of edges in EEG functional connectivity research to provide informational
density in the topology.
Specific patterns of brain activity during sleep and waking are recorded in the electroencephalogram (EEG). Time-frequency analyses methods have been widely used to analyse the EEG and identified characteristic oscillations for each vigilance state (VS) i.e., wakefulness, rapid-eye movement (REM) and non-rapid-eye movement (NREM) sleep. However, other aspects such as change of patterns associated with brain dynamics may not be captured unless a non-linear-based analysis method is used. In this pilot study, Permutation Lempel-Ziv Complexity (PLZC), a novel symbolic dynamics analysis method, was used to characterise the changes in the EEG in sleep and wakefulness during baseline and recovery from sleep deprivation. The results obtained with PLZC were contrasted with a related non-linear method, Lempel-Ziv complexity (LZC). Both measure the emergence of new patterns. However, LZC is dependent on the absolute amplitude of the EEG, while PLZC is only dependent on the relative amplitude due to symbolisation procedure and thus, more resistant to noise. We showed that PLZC discriminates activated brain states associated with wakefulness and REM sleep, which both displayed higher complexity, compared to NREM sleep. Additionally, significantly lower PLZC values were measured in NREM sleep during the recovery period following SD compared to baseline, suggesting a reduced emergence of new activity patterns in the EEG. These findings were validated using PLZC on surrogate data. By contrast, LZC was merely reflecting changes in the spectral composition of the EEG. Overall, this study implies that PLZC is a robust non-linear complexity measure, which is not dependent on amplitude variations in the signal, and which may be useful to further assess EEG alterations induced by environmental or pharmacological manipulations.
Alzheimer?s disease (AD) is the most prevalent form of dementia in the world, which is
characterised by the loss of neurones and the build-up of plaques in the brain, causing progressive
symptoms of memory loss and confusion. Although definite diagnosis is only possible by necropsy,
differential diagnosis with other types of dementia is still needed. An electroencephalogram (EEG)
is a cheap, portable, non-invasive method to record brain signals. Previous studies with non-linear
signal processing methods have shown changes in the EEG due to AD, which is characterised
reduced complexity and increased regularity. EEGs from 11 AD patients and 11 age-matched control
subjects were analysed with Fuzzy Entropy (FuzzyEn), a non-linear method that was introduced as an
improvement over the frequently used Approximate Entropy (ApEn) and Sample Entropy (SampEn)
algorithms. AD patients had significantly lower FuzzyEn values than control subjects (p < 0.01) at
electrodes T6, P3, P4, O1, and O2. Furthermore, when diagnostic accuracy was calculated using
Receiver Operating Characteristic (ROC) curves, FuzzyEn outperformed both ApEn and SampEn,
reaching a maximum accuracy of 86.36%. These results suggest that FuzzyEn could increase the
insight into brain dysfunction in AD, providing potentially useful diagnostic information. However,
results depend heavily on the input parameters that are used to compute FuzzyEn.
Legislation in Europe requires children younger than 12 or shorter than 150cm to use a child restraint system (CRS) when travelling in motor vehicles. These have decreased the number of children fatalities in car accidents, yet certain body regions such as the abdomen are still prone to severe trauma. One of the mechanisms by which injuries to this body region occur is called "submarining". It involves the seat-belt slipping over the pelvis onto the abdomen during rapid deceleration, therefore afflicting the abdomen with injurious loads.
The motion characterising submarining with CRS is as of yet poorly defined, and although CRS are assessed for protection level, there is currently no established identification criteria for submarining.
As part of the Enabling Protection for Older Children project (EPOCh), standard frontal impact sled tests of 10 CRS (6 high back booster seats; 4 booster cushions) with 10 year old anthropomorphic testing devices (ATDs) were analysed qualitatively, transversally and longitudinally for submarining detection. The methods used included video analysis, descriptive and inferential statistics, principal component analysis, time series analysis, as well as multiple linear regression and logistic regression, applied on both ATD trajectories and ATD instrumentation recordings (dynamic data).
From the videos, trajectories and dynamic data, submarining motion is shown to embody an exaggerated slouching movement regardless of the CRS type. The observations and quantitative results confirm that exaggerated forward knee displacement is characteristic of submarining, as well as a very strong criterion for the latter's detection. No other individual trajectory or dynamic variable distinguish clear submarining behaviour, however the regression analyses on dynamic variables establish the association of the pelvis, lumbar and chest as representative of the knee displacement, and demonstrate the combination's capacity to distinguish submarining cases.
These findings establish the complexity of the movement involved in submarining and the potential of using current ATD instrumentation for its assessment with CRS. This opens a path for an integral approach to ATD movement in CRS appraisal and suggests considering pelvis, lumbar and chest motion control for submarining prevention.
Sleep is an essential physiological phenomenon which is regulated by fundamental sleep-wake cycles. Sleep is formed of non-rapid-eye-movement (NREM) and REM stages where NREM and REM sleep stages alternate with wakefulness in a whole night sleep (i.e., typically consists of several sleep-wake cycles). During sleep, brain is active and the activity also alters with changing vigilance states (VS). Furthermore, physiological or external changes in the brain structure might influence brain activity during sleep. Effects of ageing, sex differences, and pharmacological manipulations have been widely investigated in sleep research using Fourier Transform. However, the use of non-linear analysis techniques might be more suitable in analysing non-linear and non-stationary signals (e.g., electroencephalogram (EEG)). Therefore, non-linear analysis has been used within this PhD with the hypothesis that these methods might reveal hidden characteristics in the changing brain signals that are difficult to detect with traditional EEG power spectral density analysis. The use of non-linear analysis techniques (e.g., symbolic dynamic analysis (SDA)) will allow to further dissect the physiological significance of activity-dependent changes of neuronal networks across sleep-wake cycles, as well as the significance of brain activity patterns during waking, sleep, sleep deprivation (SD), or induced by sleep-promoting drugs and pharmacological treatments.
In this PhD, rodent and human sleep EEG recordings were analysed using SDA methods: Lempel-Ziv complexity (LZC), Permutation Entropy (PE) and Permutation Lempel-Ziv complexity (PLZC). All the methods were able characterise different VS with wakefulness and REM sleep resulting in higher measures of complexity compared to NREM sleep suggesting an active state of the brain in these VS. This was measured in all datasets assisting the hypothesis on the usefulness of these techniques in sleep research by providing the minimum requirement for sleep analysis. In addition to this, SD significantly reduced complexity in the following sleep period supporting the compensation process for the lost sleep by the increased in slow wave activity which was reflected as reduced complexity in this study. Furthermore, a low dose tiagabine administration?s sleep compensation promoting effect was found in mice.
Moreover, ageing was identified as a main effect on changes in brain activity. These changes were more pronounced in the old age where complexity was significantly lower compared to young age. On the one hand, this was found with all three methods and contributing to the hypothesis that these techniques reveal structural dynamic changes due to physiological alterations. On the other hand, no significant differences in complexity across genders were found suggesting the underlying mechanisms to maintain sleep-wake cycles are similar for men and women. This finding with further investigation might corroborate to question the need to use both genders in drug trials. Furthermore, significant changes in brain activity were found at different times of the sleep period highlighting the changes occurring within VS as sleep progresses. This also has an impact on the way sleep stages are scored and investigated which are influenced by different brain activity levels within each VS throughout the entire sleep.
All in all, this study achieved to support its hypothesis of determining the changes in brain activity as a complexity measure by characterising sleep under physiological and pharmacologically induced EEG datasets in mice and in humans. The study was a novel application to analyse sleep in these conditions. However, with further analysis performed on larger datasets, its findings together with surrogate data analysis proved SDA techniques? robust usability which can complement the gold standard FT analysis in sleep research.
Maturation and ageing, which can be characterised by the dynamic changes in brain morphology, can have an impact on the physiology of the brain. As such, it is possible that these changes can have an impact on the magnetic activity of the brain recorded using magnetoencephalography. In this study changes in the resting state brain (magnetic) activity due to healthy ageing were investigated by estimating the complexity of magnetoencephalogram (MEG) signals. The main aim of this study was to identify if the complexity of background MEG signals changed significantly across the human lifespan for both males and females. A sample of 177 healthy participants (79 males and 98 females aged between 21 and 80 and grouped into 3 categories i.e., early-, mid- and late-adulthood) was used in this investigation. This investigation also extended to evaluating if complexity values remained relatively stable during the 5 min recording. Complexity was estimated using permutation Lempel-Ziv complexity, a recently introduced complexity metric, with a motif length of 5 and a lag of 1. Effects of age and gender were investigated in the MEG channels over 5 brain regions, i.e., anterior, central, left lateral, posterior, and, right lateral, with highest complexity values observed in the signals recorded by the channels over the anterior and central regions of the brain. Results showed that while changes due to age had a significant effect on the complexity of the MEG signals recorded over 5 brain regions, gender did not have a significant effect on complexity values in all age groups investigated. Moreover, although some changes in complexity were observed between the different minutes of recording, due to the small magnitude of the changes it was concluded that practical significance might outweigh statistical significance in this instance. The results from this study can contribute to form a fingerprint of the characteristics of healthy ageing in MEGs that could be useful when investigating changes to the resting state activity due to pathology.
Symbolic dynamic analysis (SDA) methods have been applied to biomedical signals and have been proven efficient in characterising differences in the electroencephalogram (EEG) in various conditions (e.g., epilepsy, Alzheimer?s and Parkinson?s diseases). In this study, we investigated the use of SDA on EEGs recorded during sleep. Lempel-Ziv Complexity (LZC), Permutation Entropy (PE), Permutation Lempel-Ziv Complexity (PLZC), as well as power spectral analysis based on the fast Fourier transform (FFT), were applied to 8-h sleep EEG recordings in healthy men (n=31) and women (n=29), aged 20-74 years. The results of the SDA methods and FFT analysis were compared and the effects of age and sex were investigated. Surrogate data were used to determine whether the findings with SDA methods truly reflected changes in non-linear dynamics of the EEG and not merely changes in the power spectrum. The surrogate data analysis showed that LZC merely reflected spectral changes in EEG activity, whereas PE and PLZC reflected genuine changes in the non-linear dynamics of the EEG. All three SDA techniques distinguished the vigilance states (i.e. wakefulness, REM sleep, NREM sleep and its sub stages: stage 1, stage 2 and slow wave sleep). Complexity of the sleep EEG increased with ageing. Sex on the other hand did not affect the complexity values assessed with any of these three SDA methods, even though FFT detected sex differences. This study shows that SDA provides additional insights into the dynamics of sleep EEG and how it is affected by ageing.
Symbolic dynamic analysis (SDA) methods have been applied to biomedical signals and have been proven efficient in characterising differences in the electroencephalogram (EEG) in various conditions (e.g., epilepsy, Alzheimer?s, and Parkinson?s diseases). In this study, we investigated the use of SDA on EEGs recorded during sleep. Lempel-Ziv complexity (LZC), permutation entropy (PE), and permutation Lempel-Ziv complexity (PLZC), as well as power spectral analysis based on the fast Fourier transform (FFT), were applied to 8-h sleep EEG recordings in healthy men (n=31) and women (n=29), aged 20-74 years. The results of the SDA methods and FFT analysis were compared and the effects of age and sex were investigated. Surrogate data were used to determine whether the findings with SDA methods truly reflected changes in nonlinear dynamics of the EEG and not merely changes in the power spectrum. The surrogate data analysis showed that LZC merely reflected spectral changes in EEG activity, whereas PE and PLZC reflected genuine changes in the nonlinear dynamics of the EEG. All three SDA techniques distinguished the vigilance states (i.e., wakefulness, REM sleep, NREM sleep, and its sub-stages: stage 1, stage 2, and slow wave sleep). Complexity of the sleep EEG increased with ageing. Sex on the other hand did not affect the complexity values assessed with any of these three SDA methods, even though FFT detected sex differences. This study shows that SDA provides additional insights into the dynamics of sleep EEG and how it is affected by ageing.
Cuesta-Frau David, Novák Daniel, Burda Vacláv, Abasolo Daniel, Adjei Tricia, Varela Manuel, Vargas Borja, Mraz Milos, Kavalkova Petra, Benes Marek, Haluzik Martin (2019) Influence of Duodenal?Jejunal Implantation on Glucose Dynamics: A Pilot Study Using Different Nonlinear Methods,Complexity20196070518pp. 1-10
Hindawi Publishing Corporation
Diabetes is a disease of great and rising prevalence, with the obesity epidemic being a significant contributing risk factor. Duodenal?jejunal bypass liner (DJBL) is a reversible implant that mimics the effects of more aggressive surgical procedures, such as gastric bypass, to induce weight loss. We hypothesized that DJBL also influences the glucose dynamics in type II diabetes, based on the induced changes already demonstrated in other physiological characteristics and parameters. In order to assess the validity of this assumption, we conducted a quantitative analysis based on several nonlinear algorithms (Lempel?Ziv Complexity, Sample Entropy, Permutation Entropy, and modified Permutation Entropy), well suited to the characterization of biomedical time series. We applied them to glucose records drawn from two extreme cases available of DJBL implantation: before and after 10 months. The results confirmed the hypothesis and an accuracy of 86.4% was achieved with modified Permutation Entropy. Other metrics also yielded significant classification accuracy results, all above 70%, provided a suitable parameter configuration was chosen. With the Leave?One?Out method, the results were very similar, between 72% and 82% classification accuracy. There was also a decrease in entropy of glycaemia records during the time interval studied. These findings provide a solid foundation to assess how glucose metabolism may be influenced by DJBL implantation and opens a new line of research in this field.
The analysis of resting-state brain activity recording in magnetoencephalograms (MEGs) with new algorithms of symbolic dynamics analysis could help obtain a deeper insight into the functioning of the brain and identify potential differences between males and females. Permutation Lempel-Ziv complexity (PLZC), a recently introduced non-linear signal processing algorithm based on symbolic dynamics, was used to evaluate the complexity of MEG signals in source space. PLZC was estimated in a broad band of frequencies (2?45 Hz), as well as in narrow bands (i.e., theta (4?8 Hz), alpha (8?12 Hz), low beta (12?20 Hz), high beta (20?30 Hz), and gamma (30?45 Hz)) in a sample of 98 healthy elderly subjects (49 males, 49 female) aged 65?80 (average age of 72.71 ± 4.22 for males and 72.67 ± 4.21 for females). PLZC was significantly higher for females than males in the high beta band at posterior brain regions including the precuneus, and the parietal and occipital cortices. Further statistical analyses showed that higher complexity values over highly overlapping regions than the ones mentioned above were associated with larger hippocampal volumes only in females. These results suggest that sex differences in healthy aging can be identified from the analysis of magnetoencephalograms with novel signal processing methods.
The brain is a very delicate, but sophisticated and complex organ of the body. During life, the brain grows and develops, matures, and then ages like all other organs of the body. This maturation and ageing comes with various physiological and anatomical changes which have an impact on the background activity of the brain. Brain activity can be recorded using various techniques including magnetoencephalography and magnetic resonance imaging, and these recordings combined with signal processing can be useful to characterise the changes in the brain that can be a result of the ageing process. By analysing various aspects of the brain, such as functional and effective connectivity, entropy and complexity, the state of the resting brain at different ages can be understood with greater detail. Furthermore, these analyses can be used to investigate the resting state brain networks that are present in the brain, as well as their topology. This information can be combined with network analysis techniques such as graph theory and used to understand the manner in which the brain both matures and ages.
This research made use of two magnetoencephalogram (MEG) databases recorded from both males and females. The first, containing resting state MEGs (rMEGs) from 220 healthy volunteers (aged 7-84), was used as the main database in this thesis to investigate the effects of age on rMEG signals throughout life. The aim of this research was to make use of rMEG signals and signal processing techniques to determine the effects of healthy ageing on the brain throughout life. It was hypothesised that the effects of age are identifiable using advanced signal processing techniques. Thus, the effects of age on linear interactions, causality, synchronisation, information flow, entropy and complexity, were investigated in the 148 MEG channels lying over 5 brain regions (anterior, central, left lateral, posterior, and, right lateral) using multiple linear and non-linear analysis techniques (namely: Pearson?s correlation, coherence, Granger causality, phase slope index, rho index, transfer entropy, synchronisation likelihood, Lempel-Ziv complexity, permutation Lempel-Ziv complexity, permutation entropy, and, modified permutation entropy). Additionally, graph theory principles were used to evaluate different network components such as integration (global efficiency), segregation (clustering coefficient and modularity), centrality (betweenness), and resilience (strength and assortativity) so as to obtain an understanding of the construct of the resting brain network. Moreover, complex network analysis was also used to determine the overall network topology of the brain network and how this changed at different stages in life. Gender effects were also studied so as to identify if there were any significant differences between males and females at different stages of life. Results from these analyses showed that the healthy resting brain has low effective and functional connectivity, relatively low complexity and entropy, as well as no significant detectable direction of information flow. Therefore this showed that there is very little synchronous or simultaneously occurring information in the rMEG time series. Thus, during rest, the brain resembles a system in limbo/phase transition, with low effective and functional connectivity, relatively low complexity and entropy, and no significant detectable direction of information flow.
The second database used in this research project was obtained during a collaboration visit to the Cognitive and Computational Neuroscience laboratory at the Centre for Biomedical Technology- Universidad Politécnica de Madrid (CTB-UPM). This database was made up of rMEGs recorded from 199 healthy volunteers (aged 60-80), and the focus of this additional set of analyses was to identify differences between the rMEG signals recorded from healthy individuals, those with subjective cognitive decline as well as those with
Biomedical signals are measurable time series that describe a physiological state of a biological system. Entropy algorithms have been previously used to quantify the complexity of biomedical signals, but there is a need to understand the relationship of entropy to signal processing concepts. In this study, ten synthetic signals that represent widely encountered signal structures in the field of signal processing were created to interpret permutation, modified permutation, sample, quadratic sample and fuzzy entropies. Subsequently, the entropy algorithms were applied to two different databases containing electroencephalogram (EEG) signals from epilepsy studies. Transitions from randomness to periodicity were successfully detected in the synthetic signals, while significant differences in EEG signals were observed based on different regions and states of the brain. In addition, using results from one entropy algorithm as features and the k-nearest neighbours algorithm, maximum classification accuracies in the first EEG database ranged from 63% to 73.5%, while these values increased by approximately 20% when using two different entropies as features. For the second database, maximum classification accuracy reached 62.5% using one entropy algorithm, while using two algorithms as features further increased that by 10%. Embedding entropies (sample, quadratic sample and fuzzy entropies) are found to outperform the rest of the algorithms in terms of sensitivity and show greater potential by considering the fine-tuning possibilities they offer. On the other hand, permutation and modified permutation entropies are more consistent across different input parameter values and considerably faster to calculate.
To analyse magnetoencephalogram (MEG) signals with Lempel-Ziv Complexity (LZC) to identify the regions of the brain showing changes related to cognitive decline and AlzheimeUs Disease (AD).
LZC was used to study MEG signals in the source space from 99 participants (36 male, 63 female, average age: 71.82 ± 4.06) in three groups (33 subjects per group): healthy (control) older adults, older adults with subjective cognitive decline (SCD), and adults with mild cognitive impairment (MCI). Analyses were performed in broadband (2?45 Hz) and in classic narrow bands (theta (4?8 Hz), alpha (8?12 Hz), low beta (12?20 Hz), high beta (20?30 Hz), and, gamma (30?45 Hz)).
LZC was significantly lower in subjects with MCI than in those with SCD. Moreover, subjects with MCI had significantly lower MEG complexity than controls and SCD subjects in the beta frequency band. Lower complexity was correlated with smaller hippocampal volumes.
Brain complexity ? measured with LZC ? decreases in MCI patients when compared to SCD and healthy controls. This decrease is associated with a decrease in hippocampal volume, a key feature in AD progression.
This is the first study to date characterising the changes of brain activity complexity showing the specific spatial pattern of the alterations as well as the morphological correlations throughout preclinical stages of AD.