Interpretation of the Lempel-Ziv complexity measure in the context of biomedical signal analysis

The effect of age on balancing behavior: complexity analysis of mediolateral force trajectories

OBJECTIVE

We quantified, via complexity analysis, the postural stability of healthy people from a wide age range.

APPROACH

Thirty-five healthy people aged 18-72 performed three tasks while balancing on one foot on a force plate: standard balancing task, mental task (balancing while answering basic arithmetic questions), and knot-tying task (balancing while tying two knots in a piece of ribbon). Mediolateral force trajectories were analyzed to determine control strategy via Hurst exponents, Lyapunov exponents, Kolmogorov complexity, root mean square, and phase-space plots.

MAIN RESULTS

We found increased pattern repetition in balancing with increased age, as evidenced by the emergence of a double attractor pattern in phase-space plots and the increase of Hurst exponents with age from approximately 0.3 to 0.8.

SIGNIFICANCE

As people age, they tend to develop strong feed-forward control strategies for balancing, and lose the complexity of micro movements intrinsic to young age. There is an open-loop control strategy for balancing that emerges in older adulthood, and there are attractors inherent to balancing which begin to develop in middle age.

Comparing EEG Nonlinearity in Deficit and Nondeficit Schizophrenia Patients: Preliminary Data

Electroencephalogram (EEG) contains valuable information obtained noninvasively that can be used for assessment of brain's processing capacity of patients with psychiatric disorders. The purpose of the present work was to evaluate possible differences in EEG complexity between deficit (DS) and nondeficit (NDS) subtypes of schizophrenia as a reflection of the cognitive processing capacities in these groups. A particular nonlinear metric known as Lempel-Ziv complexity (LZC) was used as a computational tool in order to determine the randomness in EEG alpha band time series from 3 groups (deficit schizophrenia [n = 9], nondeficit schizophrenia [n = 10], and healthy controls [n = 10]) according to time series randomness. There was a significant difference in frontal EEG complexity between the DS and NDS subgroups ( p = .013), with DS group showing less complexity. A significant positive correlation was found between LZC values and Positive and Negative Syndrome Scale (PANSS) general psychopathology scores (ie, larger frontal EEG complexity correlated with more severe psychopathology), explained partially by the emotional component subscore of the PANSS. These findings suggest that cognitive processing occurring in the frontal networks in DS is less complex compared to NDS patients as reflected by EEG complexity measures. The data also suggest that there may be a relationship between the degree of emotionality and the complexity of the frontal EEG signal.

Altered cerebral blood flow velocity features in fibromyalgia patients in resting-state conditions

The aim of this study is to characterize in resting-state conditions the cerebral blood flow velocity (CBFV) signals of fibromyalgia patients. The anterior and middle cerebral arteries of both hemispheres from 15 women with fibromyalgia and 15 healthy women were monitored using Transcranial Doppler (TCD) during a 5-minute eyes-closed resting period. Several signal processing methods based on time, information theory, frequency and time-frequency analyses were used in order to extract different features to characterize the CBFV signals in the different vessels. Main results indicated that, in comparison with control subjects, fibromyalgia patients showed a higher complexity of the envelope CBFV and a different distribution of the power spectral density. In addition, it has been observed that complexity and spectral features show correlations with clinical pain parameters and emotional factors. The characterization features were used in a lineal model to discriminate between fibromyalgia patients and healthy controls, providing a high accuracy. These findings indicate that CBFV signals, specifically their complexity and spectral characteristics, contain information that may be relevant for the assessment of fibromyalgia patients in resting-state conditions.

A comparison between swallowing sounds and vibrations in patients with dysphagia

The cervical auscultation refers to the observation and analysis of sounds or vibrations captured during swallowing using either a stethoscope or acoustic/vibratory detectors. Microphones and accelerometers have recently become two common sensors used in modern cervical auscultation methods. There are open questions about whether swallowing signals recorded by these two sensors provide unique or complementary information about swallowing function; or whether they present interchangeable information. This study aims to compare of swallowing signals recorded by a microphone and a tri-axial accelerometer from 72 patients (mean age 63.94 ± 12.58 years, 42 male, 30 female), who had videofluoroscopic examination. The participants swallowed one or more boluses of thickened liquids of different consistencies, including thin liquids, nectar-thick liquids, and pudding. A comfortable self-selected volume from a cup or a controlled volume by the examiner from a 5 ml spoon was given to the participants. A broad feature set was extracted in time, information-theoretic, and frequency domains from each of 881 swallows presented in this study. The swallowing sounds exhibited significantly higher frequency content and kurtosis values than the swallowing vibrations. In addition, the Lempel-Ziv complexity was lower for swallowing sounds than those for swallowing vibrations. To conclude, information provided by microphones and accelerometers about swallowing function are unique and these two transducers are not interchangeable. Consequently, the selection of transducer would be a vital step in future studies.

Seizure prediction in patients with focal hippocampal epilepsy

OBJECTIVE

We evaluated the performance of our previously developed seizure prediction approach on thirty eight seizures from ten patients with focal hippocampal epilepsy.

METHODS

The seizure prediction system was developed based on the extraction of correlation dimension, correlation entropy, noise level, Lempel-Ziv complexity, largest Lyapunov exponent, and nonlinear interdependence from segments of intracranial EEG.

RESULTS

Our results showed an average sensitivity of 86.7% and 92.9%, an average false prediction rate of 0.126 and 0.096/h, and an average minimum prediction time of 14.3 and 33.3min, respectively, using seizure occurrence periods of 30 and 50min and a seizure prediction horizon of 10s. Two-third of the analyzed seizures showed significantly increased complexity in periods prior to the seizures in comparison with baseline. In four patients, strong bidirectional connectivities between epileptic contacts and the surrounding areas were observed. However, in five patients, unidirectional functional connectivities in preictal periods were observed from remote areas to epileptogenic zones.

CONCLUSIONS

Overall, preictal periods in patients with focal hippocampal epilepsy were characterized with patient-specific changes in univariate and bivariate nonlinear measures.

SIGNIFICANCE

The spatio-temporal characterization of preictal periods may help to better understand the mechanism underlying seizure generation in patients with focal hippocampal epilepsy.

Refined multiscale fuzzy entropy based on standard deviation for biomedical signal analysis

Multiscale entropy (MSE) has been a prevalent algorithm to quantify the complexity of biomedical time series. Recent developments in the field have tried to alleviate the problem of undefined MSE values for short signals. Moreover, there has been a recent interest in using other statistical moments than the mean, i.e., variance, in the coarse-graining step of the MSE. Building on these trends, here we introduce the so-called refined composite multiscale fuzzy entropy based on the standard deviation (RCMFE_σ) and mean (RCMFE_μ) to quantify the dynamical properties of spread and mean, respectively, over multiple time scales. We demonstrate the dependency of the RCMFE_σ and RCMFE_μ, in comparison with other multiscale approaches, on several straightforward signal processing concepts using a set of synthetic signals. The results evidenced that the RCMFE_σ and RCMFE_μ values are more stable and reliable than the classical multiscale entropy ones. We also inspect the ability of using the standard deviation as well as the mean in the coarse-graining process using magnetoencephalograms in Alzheimer’s disease and publicly available electroencephalograms recorded from focal and non-focal areas in epilepsy. Our results indicated that when the RCMFE_μ cannot distinguish different types of dynamics of a particular time series at some scale factors, the RCMFE_σ may do so, and vice versa. The results showed that RCMFE_σ-based features lead to higher classification accuracies in comparison with the RCMFE_μ-based ones. We also made freely available all the Matlab codes used in this study at 10.7488/ds/1477.

Dispersion entropy for the analysis of resting-state MEG regularity in Alzheimer's disease

Alzheimer's disease (AD) is a progressive degenerative brain disorder affecting memory, thinking, behaviour and emotion. It is the most common form of dementia and a big social problem in western societies. The analysis of brain activity may help to diagnose this disease. Changes in entropy methods have been reported useful in research studies to characterize AD. We have recently proposed dispersion entropy (DisEn) as a very fast and powerful tool to quantify the irregularity of time series. The aim of this paper is to evaluate the ability of DisEn, in comparison with fuzzy entropy (FuzEn), sample entropy (SampEn), and permutation entropy (PerEn), to discriminate 36 AD patients from 26 elderly control subjects using resting-state magnetoencephalogram (MEG) signals. The results obtained by DisEn, FuzEn, and SampEn, unlike PerEn, show that the AD patients' signals are more regular than controls' time series. The p-values obtained by DisEn, FuzEn, SampEn, and PerEn based methods demonstrate the superiority of DisEn over PerEn, SampEn, and PerEn. Moreover, the computation time for the newly proposed DisEn-based method is noticeably less than for the FuzEn, SampEn, and PerEn based approaches.

Changes in the EEG amplitude as a biomarker for early detection of Alzheimer's disease

The rapid increase in the number of older people with Alzheimer's disease (AD) and other forms of dementia represents one of the major challenges to the health and social care systems. Early detection of AD makes it possible for patients to access appropriate services and to benefit from new treatments and therapies, as and when they become available. The onset of AD starts many years before the clinical symptoms become clear. A biomarker that can measure the brain changes in this period would be useful for early diagnosis of AD. Potentially, the electroencephalogram (EEG) can play a valuable role in early detection of AD. Damage in the brain due to AD leads to changes in the information processing activity of the brain and the EEG which can be quantified as a biomarker. The objective of the study reported in this paper is to develop robust EEG-based biomarkers for detecting AD in its early stages. We present a new approach to quantify the slowing of the EEG, one of the most consistent features at different stages of dementia, based on changes in the EEG amplitudes (ΔEEGA). The new approach has sensitivity and specificity values of 100% and 88.88%, respectively, and outperformed the Lempel-Ziv Complexity (LZC) approach in discriminating between AD and normal subjects.

Comparison of oscillatory activity in subthalamic nucleus in Parkinson's disease and dystonia

OBJECTIVES

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) has been successfully used to treat both Parkinson's disease (PD) and dystonia. Local field potentials (LFPs) recorded from the STN of PD patients demonstrate prominent beta frequency band activity. It is unclear whether such activity occurs in the STN in dystonia, and, if not, whether dystonia has another distinctive neural population activity in the STN.

METHODS

Twelve patients with PD, and eight patients with dystonia underwent DBS electrode implantation targeting the STN. Seven dystonia patients were off medication and one was on aripiprazole and clonazepam. LFPs were recorded from the DBS electrodes in PD in the on/off medication states and in dystonia. Power spectra and temporal dynamics measured by the with Lempel-Ziv complexity of the LFPs were compared among these states.

RESULTS

Normalised power spectra and Lempel-Ziv complexity of subthalamic LFPs differed between dystonia off and PD on/off, and between PD off and on over the low frequency, beta and high gamma bands. Patients with dystonia and off medication had lower beta power but higher low frequency and high gamma power than PD. Spectral power in the low beta frequency (11-20Hz) range was attenuated in medicated PD.

CONCLUSION

The results suggest that dystonia and PD are characterized by different patterns of oscillatory activities even within the same nucleus, and exaggerated beta activity may relate to hypo-dopaminergic status.

Aging and cardiovascular complexity: effect of the length of RR tachograms

As we age, our hearts undergo changes that result in a reduction in complexity of physiological interactions between different control mechanisms. This results in a potential risk of cardiovascular diseases which are the number one cause of death globally. Since cardiac signals are nonstationary and nonlinear in nature, complexity measures are better suited to handle such data. In this study, three complexity measures are used, namely Lempel-Ziv complexity (LZ), Sample Entropy (SampEn) and Effort-To-Compress (ETC). We determined the minimum length of RR tachogram required for characterizing complexity of healthy young and healthy old hearts. All the three measures indicated significantly lower complexity values for older subjects than younger ones. However, the minimum length of heart-beat interval data needed differs for the three measures, with LZ and ETC needing as low as 10 samples, whereas SampEn requires at least 80 samples. Our study indicates that complexity measures such as LZ and ETC are good candidates for the analysis of cardiovascular dynamics since they are able to work with very short RR tachograms.

Anesthesia-related changes in information transfer may be caused by reduction in local information generation

In anesthesia research it is an open question how general anesthetics lead to loss of consciousness (LOC). It has been proposed that LOC may be caused by the disruption of cortical information processing, preventing information integration. Therefore, recent studies investigating information processing under anesthesia focused on changes in information transfer, measured by transfer entropy (TE). However, often this complex technique was not applied rigorously, using time series in symbolic representation, or using TE differences without accounting for neural conduction delays, or without accounting for signal history. Here, we used current best-practice in TE estimation to investigate information transfer under anesthesia: We conducted simultaneous recordings in primary visual cortex (V1) and prefrontal cortex (PFC) of head-fixed ferrets in a dark environment under different levels of anesthesia (awake, 0.5% isoflurane, 1.0 % isoflurane). To elucidate reasons for changes in TE, we further quantified information processing within brain areas by estimating active information storage (AIS) as an estimator of predictable information, and Lempel-Ziv complexity (LZC) as an estimator of signal entropy. Under anesthesia, we found a reduction in information transfer (TE) between PFC and V1 with a stronger reduction for the feedback direction (PFC to V1), validating previous results. Furthermore, entropy (LZC) was reduced and activity became more predictable as indicated by higher values of AIS. We conclude that higher anesthesia concentrations indeed lead to reduced inter-areal information transfer, which may be partly caused by decreases in local entropy and increases in local predictability. In revealing a possible reason for reduced TE that is potentially independent of inter-areal coupling, we demonstrate the value of directly quantifying information processing in addition to focusing on dynamic properties such as coupling strength.

Using Lempel-Ziv Complexity to Assess ECG Signal Quality

The poor quality of wireless electrocardiography (ECG) recordings can lead to misdiagnosis and waste of medical resources. This study presents an interpretation of Lempel–Ziv (LZ) complexity in terms of ECG quality assessment, and verifies its performance on real ECG signals. Firstly, LZ complexities for typical signals, namely high-frequency (HF) noise, low-frequency (LF) noise, power-line (PL) noise, impulse (IM) noise, clean artificial ECG signals, and ECG signals with various types of noise added (ECG plus HF, LF, PL, and IM noise, respectively) were analyzed. Then, the effects of noise, signal length, and signal-to-noise ratio (SNR) on the LZ complexity of ECG signals were analyzed. The simulation results show that LZ complexity for HF noise was obviously different from those for PL and LF noise. The LZ value can be used to determine the presence of HF noise. ECG plus HF noise had the highest LZ values. Other types of noise had low LZ values. Signal lengths of over 40 s had only a small effect on LZ values. The LZ values for ECG plus all types of noise increased monotonically with decreasing SNR except for LF and PL noise. For the test of real ECG signals plus three types of noise, namely muscle artefacts (MAs), baseline wander (BW), and electrode motion (EM) artefacts, LZ complexity varied obviously with increasing MA but not for BW and EM noise. This study demonstrates that LZ complexity is sensitive to noise level (especially for HF noise) and can thus be a valuable reference index for the assessment of ECG signal quality.

Lempel-Ziv and multiscale Lempel-Ziv complexity in depression

There is a high demand for objective indicators in diagnosis of depression as diagnosis of depression is still based on psychiatrist's subjective judgment. A nonlinear method Lempel Ziv Complexity (LZC) has been previously successfully used for detection of neuronal or mental disorders based on electroencephalographic (EEG) signals. However, the method overlooks the high frequency content of EEG signals. Therefore, this study is aimed to find out whether the use of Multiscale Lempel Ziv Complexity (MLZC), considering also high frequencies, could overcome the limitations of LZC and better differentiate depression. In current study the EEG recordings were carried out on the groups of depressive and healthy subjects of 11 volunteers each. The LZC and MLZC were calculated on resting EEG signals in eyes open condition from 30 channels at a length of 2 minutes. The results revealed the incapability of traditional LZC to differentiate depressive subjects from healthy controls in eyes open condition, while MLZC differentiated two groups in numerous channels at different frequencies, giving the highest classification accuracy in channel F3 (86 %) at frequencies 9 and 15.5 Hz. The results indicate that the high frequency information, which is lost in calculation of traditional LZC, has a great value in differentiating between depressive and control groups.

A novel encoding Lempel-Ziv complexity algorithm for quantifying the irregularity of physiological time series

BACKGROUND AND OBJECTIVE

The Lempel-Ziv (LZ) complexity and its variants have been extensively used to analyze the irregularity of physiological time series. To date, these measures cannot explicitly discern between the irregularity and the chaotic characteristics of physiological time series. Our study compared the performance of an encoding LZ (ELZ) complexity algorithm, a novel variant of the LZ complexity algorithm, with those of the classic LZ (CLZ) and multistate LZ (MLZ) complexity algorithms.

METHODS AND RESULTS

Simulation experiments on Gaussian noise, logistic chaotic, and periodic time series showed that only the ELZ algorithm monotonically declined with the reduction in irregularity in time series, whereas the CLZ and MLZ approaches yielded overlapped values for chaotic time series and time series mixed with Gaussian noise, demonstrating the accuracy of the proposed ELZ algorithm in capturing the irregularity, rather than the complexity, of physiological time series. In addition, the effect of sequence length on the ELZ algorithm was more stable compared with those on CLZ and MLZ, especially when the sequence length was longer than 300. A sensitivity analysis for all three LZ algorithms revealed that both the MLZ and the ELZ algorithms could respond to the change in time sequences, whereas the CLZ approach could not. Cardiac interbeat (RR) interval time series from the MIT-BIH database were also evaluated, and the results showed that the ELZ algorithm could accurately measure the inherent irregularity of the RR interval time series, as indicated by lower LZ values yielded from a congestive heart failure group versus those yielded from a normal sinus rhythm group (p < 0.01).

Effect of the average delay and mean connectivity of the Kuramoto model on the complexity of the output electroencephalograms

Cognitive functions result from the interplay of distributed brain areas operating in large-scale networks. These networks can be modelled with a number of parameters that represent their underlying dynamics. One particularly fruitful model to simulate key aspects of the large-scale brain networks is the Kuramoto model, which simulates the phase evolution of several weakly coupled oscillators that represent the mean oscillatory behavior of different cortical regions. Here, we inspected the dependency of two widespread nonlinear complexity markers, Sample Entropy (SampEn) and Lempel-Ziv Complexity (LZC), on EEG activity generated with a Kuramoto phase model where the time delay and connectivity strength among oscillators varied. We also added different levels of noise to the electroencephalogram (EEG) signals. Our results indicated that both complexity metrics reflected the changes in the delays and global synchrony levels, but we found that SampEn was slightly more sensitive to the state transition and its results were less affected by the presence of noise. These results help in the effort to understand the dynamics of EEG recordings and their relationship to large-scale networks.

EEG analysis via multiscale Lempel-Ziv complexity for seizure detection

Robust seizure detection and seizure prediction continues to be a challenge. Lempel-Ziv Complexity (LZC) is one of the features that has shown to be relevant in seizure detection. Recent work has shown that augmenting LZC can be beneficial to emphasize variations in amplitude or frequency when analyzing biomedical signals. In this paper, we present a first look into evaluating the feasibility of using a recently proposed feature stemmed from LZC, namely the Multiscale Lempel-Ziv Complexity (MLZC) for seizure detection. MLZC does not allow the high-frequency signal components to be overwhelmed by the low frequency signal components when calculating complexity values. We compare MLZC and LZC for identifying seizures for three cases and show MLZC can provide a clear separation between non-ictal and ictal periods for all three cases using a single threshold over 7 recordings and 7 seizures per patient, whereas LZC provided such a clear separation for only one of the patients.

Anatomical Directional Dissimilarities in Tri-axial Swallowing Accelerometry Signals

Swallowing accelerometry is a noninvasive approach currently under consideration as an instrumental screening test for swallowing difficulties, with most current studies focusing on the swallowing vibrations in the anterior-posterior (A-P) and superior-inferior (S-I) directions. However, the displacement of the hyolaryngeal structure during the act of swallowing in patients with dysphagia involves declination of the medial-lateral (M-L), which suggests that the swallowing vibrations in the M-L direction have the ability to reveal additional details about the swallowing function. With this motivation, we performed a broad comparison of the swallowing vibrations in all three anatomical directions. Tri-axial swallowing accelerometry signals were concurrently collected from 72 dysphagic patients undergoing videofluoroscopic evaluation of swallowing (mean age: 63.94 ± 12.58 years period). Participants swallowed one or more thickened liquids with different consistencies including thin-thick liquids, nectar-thick liquids, and pudding-thick liquids with either a comfortable self-selected volume from a cup or a controlled volume by the examiner from a 5-ml spoon. Swallows were grouped based on the viscosity of swallows and the participant's stroke history. Then, a comprehensive set of features was extracted in multiple signal domains from 881 swallows. The results highlighted inter-axis dissimilarities among tri-axial swallowing vibrations including the extent of variability in the amplitude of signals, the degree of predictability of signals, and the extent of disordered behavior of signals in time-frequency domain. First, the upward movement of the hyolaryngeal structure, representing the S-I signals, were actually more variable in amplitude and showed less predictable behavior than the sideways and forward movements, representing the A-P and M-L signals, during swallowing. Second, the S-I signals, which represent the upward movement of the hyolaryngeal structure, behaved more disordered in the time-frequency domain than the sideways movement, M-L signals, in all groups of study except for the pudding swallows in the stroke group. Third, considering the viscosity and the participant's pathology, thin liquid swallows in the nonstroke group presented the most directional differences among all groups of study. In summary, despite some directional dissimilarities, M-L axis accelerometry characteristics are similar to those of the two other axes. This indicates that M-L axis characteristics, which cannot be observed in videofluoroscopic images, can be adequately derived from the A-P and S-I axes.

Spontaneous neural activity in the primary visual cortex of retinal degenerated rats

Retinal degeneration (RD) models have been widely used to study retinal degenerative diseases for a long time. The biological and electrophysiological presentations of changes in the retina during degeneration progress have been well investigated; thus, the present study is aimed at investigating the electrophysiological effects of RD in the primary visual cortex. We extracellularly recorded the spontaneous neural activities in the primary visual cortex of RD rats. The firing rate, interspike interval (ISI) and Lempel-Ziv (LZ) complexity of spontaneous neural activities were subsequently analyzed. When compared to the control group, it was found that the neurons in primary visual cortex of the RD model fired more frequently. In addition, there was a decrease in LZ complexity of spontaneous neural firing in the RD model. These results suggest that the progress of RD may not only affect the retina itself but also the primary visual cortex, which may result in an unbalanced inhibition-excitation system as well as the decreased arising rate of new patterns of spontaneous activities.

Cognitive tasks and cerebral blood flow through anterior cerebral arteries: a study via functional transcranial Doppler ultrasound recordings

Background

Functional transcanial Doppler ultrasound (fTCD) is a convenient approach to examine cerebral blood flow velocity (CBFV) in major cerebral arteries.

Methods

In this study, the anterior cerebral artery (ACA) was insonated on both sides, that is, right ACA (R-ACA) and left ACA (L-ACA). The envelope signals (the maximum velocity) and the raw signals were analyzed during cognitive processes, i.e. word-generation tasks, geometric tasks and resting state periods separating each task. Data which were collected from 20 healthy participants were used to investigate the changes and the hemispheric functioning while performing cognitive tasks. Signal characteristics were analyzed in time domain, frequency domain and time-frequency domain.

Results

Significant results have been obtained through the use of both classic/modern methods (i.e. envelope/raw, time and frequency/information-theoretic and time-frequency domains). The frequency features extracted from the raw signals highlighted sex effects on cerebral blood flow which revealed distinct brain response during each process and during resting periods. In the time-frequency analysis, the distribution of wavelet energies on the envelope signals moved around the low frequencies during mental processes and did not experience any lateralization during cognitive tasks.

Conclusions

Even if no lateralization effects were noticed during resting-state, verbal and geometric tasks, understanding CBFV in ACA during cognitive tasks could complement information extracted from cerebral blood flow in middle cerebral arteries during similar cognitive tasks (i.e. sex effects).

A statistical analysis of cervical auscultation signals from adults with unsafe airway protection

Background

Aspiration, where food or liquid is allowed to enter the larynx during a swallow, is recognized as the most clinically salient feature of oropharyngeal dysphagia. This event can lead to short-term harm via airway obstruction or more long-term effects such as pneumonia. In order to non-invasively identify this event using high resolution cervical auscultation there is a need to characterize cervical auscultation signals from subjects with dysphagia who aspirate.

Methods

In this study, we collected swallowing sound and vibration data from 76 adults (50 men, 26 women, mean age 62) who underwent a routine videofluoroscopy swallowing examination. The analysis was limited to swallows of liquid with either thin (<5 cps) or viscous (≈300 cps) consistency and was divided into those with deep laryngeal penetration or aspiration (unsafe airway protection), and those with either shallow or no laryngeal penetration (safe airway protection), using a standardized scale. After calculating a selection of time, frequency, and time-frequency features for each swallow, the safe and unsafe categories were compared using Wilcoxon rank-sum statistical tests.

Results

Our analysis found that few of our chosen features varied in magnitude between safe and unsafe swallows with thin swallows demonstrating no statistical variation. We also supported our past findings with regard to the effects of sex and the presence or absence of stroke on cervical ausculation signals, but noticed certain discrepancies with regards to bolus viscosity.

Conclusions

Overall, our results support the necessity of using multiple statistical features concurrently to identify laryngeal penetration of swallowed boluses in future work with high resolution cervical auscultation.

detectMITE: A novel approach to detect miniature inverted repeat transposable elements in genomes

Miniature inverted repeat transposable elements (MITEs) are prevalent in eukaryotic genomes, including plants and animals. Classified as a type of non-autonomous DNA transposable elements, they play important roles in genome organization and evolution. Comprehensive and accurate genome-wide detection of MITEs in various eukaryotic genomes can improve our understanding of their origins, transposition processes, regulatory mechanisms, and biological relevance with regard to gene structures, expression, and regulation. In this paper, we present a new MATLAB-based program called detectMITE that employs a novel numeric calculation algorithm to replace conventional string matching algorithms in MITE detection, adopts the Lempel-Ziv complexity algorithm to filter out MITE candidates with low complexity, and utilizes the powerful clustering program CD-HIT to cluster similar MITEs into MITE families. Using the rice genome as test data, we found that detectMITE can more accurately, comprehensively, and efficiently detect MITEs on a genome-wide scale than other popular MITE detection tools. Through comparison with the potential MITEs annotated in Repbase, the widely used eukaryotic repeat database, detectMITE has been shown to find known and novel MITEs with a complete structure and full-length copies in the genome. detectMITE is an open source tool (https://sourceforge.net/projects/detectmite).

C. elegans locomotion analysis using algorithmic information theory

This article investigates the use of algorithmic information theory to analyse C. elegans datasets. The ability of complexity measures to detect similarity in animals' behaviours is demonstrated and their strengths are compared to methods such as histograms. Introduced quantities are illustrated on a couple of real two-dimensional C. elegans datasets to investigate the thermotaxis and chemotaxis behaviours.

Lempel-Ziv complexity analysis of one dimensional cellular automata

Lempel-Ziv complexity measure has been used to estimate the entropy density of a string. It is defined as the number of factors in a production factorization of a string. In this contribution, we show that its use can be extended, by using the normalized information distance, to study the spatiotemporal evolution of random initial configurations under cellular automata rules. In particular, the transfer information from time consecutive configurations is studied, as well as the sensitivity to perturbed initial conditions. The behavior of the cellular automata rules can be grouped in different classes, but no single grouping captures the whole nature of the involved rules. The analysis carried out is particularly appropriate for studying the computational processing capabilities of cellular automata rules.

Neural complexity in patients with poststroke depression: A resting EEG study

BACKGROUND

Poststroke depression (PSD) is one of the most common emotional disorders affecting post-stroke patients. However, the neurophysiological mechanism remains elusive. This study was aimed to study the relationship between complexity of neural electrical activity and PSD.

METHODS

Resting state eye-closed electroencephalogram (EEG) signals of 16 electrodes were recorded in 21 ischemic poststroke depression (PSD) patients, 22 ischemic poststroke non-depression (PSND) patients and 15 healthy controls (CONT). Lempel-Ziv Complexity (LZC) was used to evaluate changes in EEG complexity in PSD patients. Statistical analysis was performed to explore difference among different groups and electrodes. Correlation between the severity of depression (HDRS) and EEG complexity was determined with pearson correlation coefficients. Receiver operating characteristic (ROC) and binary logistic regression analysis were conducted to estimate the discriminating ability of LZC for PSD in specificity, sensitivity and accuracy.

RESULTS

PSD patients showed lower neural complexity compared with PSND and CONT subjects in the whole brain regions. There was no significant difference among different brain regions, and no interactions between group and electrodes. None of the LZC significantly correlated with overall depression severity or differentiated symptom severity of 7 items in PSD patients, but in stroke patients, significant correlation was found between HDRS and LZC in the whole brain regions, especially in frontal and temporal. LZC parameters used for PSD recognition possessed more than 85% in specificity, sensitivity and accuracy, suggesting the feasibility of LZC to serve as screening indicators for PSD. Increased slow wave rhythms were found in PSD patients and clearly correlation was confirmed between neuronal complexity and spectral power of the four EEG rhythms.

LIMITATIONS

Lesion location of stroke patients in the study distributed in different brain regions, and most of the PSD patients were mild or moderate in depressive severity.

CONCLUSIONS

Compared with conventional spectral analysis, complexity of neural activity using LZC was more sensitive and stationary in the measurement of abnormal brain activity in PSD patients and may offer a potential approach to facilitate clinical screening of this disease.

Assessing interactions among multiple physiological systems during walking outside a laboratory: An Android based gait monitor

Gait function is traditionally assessed using well-lit, unobstructed walkways with minimal distractions. In patients with subclinical physiological abnormalities, these conditions may not provide enough stress on their ability to adapt to walking. The introduction of challenging walking conditions in gait can induce responses in physiological systems in addition to the locomotor system. There is a need for a device that is capable of monitoring multiple physiological systems in various walking conditions. To address this need, an Android-based gait-monitoring device was developed that enabled the recording of a patient's physiological systems during walking. The gait-monitoring device was tested during self-regulated overground walking sessions of fifteen healthy subjects that included 6 females and 9 males aged 18-35 years. The gait-monitoring device measures the patient's stride interval, acceleration, electrocardiogram, skin conductance and respiratory rate. The data is stored on an Android phone and is analyzed offline through the extraction of features in the time, frequency and time-frequency domains. The analysis of the data depicted multisystem physiological interactions during overground walking in healthy subjects. These interactions included locomotion-electrodermal, locomotion-respiratory and cardiolocomotion couplings. The current results depicting strong interactions between the locomotion system and the other considered systems (i.e., electrodermal, respiratory and cardiovascular systems) warrant further investigation into multisystem interactions during walking, particularly in challenging walking conditions with older adults.

The psychosis-like effects of Δ(9)-tetrahydrocannabinol are associated with increased cortical noise in healthy humans

BACKGROUND

Drugs that induce psychosis may do so by increasing the level of task-irrelevant random neural activity or neural noise. Increased levels of neural noise have been demonstrated in psychotic disorders. We tested the hypothesis that neural noise could also be involved in the psychotomimetic effects of delta-9-tetrahydrocannabinol (Δ(9)-THC), the principal active constituent of cannabis.

METHODS

Neural noise was indexed by measuring the level of randomness in the electroencephalogram during the prestimulus baseline period of an oddball task using Lempel-Ziv complexity, a nonlinear measure of signal randomness. The acute, dose-related effects of Δ(9)-THC on Lempel-Ziv complexity and signal power were studied in humans (n = 24) who completed 3 test days during which they received intravenous Δ(9)-THC (placebo, .015 and .03 mg/kg) in a double-blind, randomized, crossover, and counterbalanced design.

RESULTS

Δ(9)-THC increased neural noise in a dose-related manner. Furthermore, there was a strong positive relationship between neural noise and the psychosis-like positive and disorganization symptoms induced by Δ(9)-THC, which was independent of total signal power. Instead, there was no relationship between noise and negative-like symptoms. In addition, Δ(9)-THC reduced total signal power during both active drug conditions compared with placebo, but no relationship was detected between signal power and psychosis-like symptoms.

CONCLUSIONS

At doses that produced psychosis-like effects, Δ(9)-THC increased neural noise in humans in a dose-dependent manner. Furthermore, increases in neural noise were related with increases in Δ(9)-THC-induced psychosis-like symptoms but not negative-like symptoms. These findings suggest that increases in neural noise may contribute to the psychotomimetic effects of Δ(9)-THC.

A cerebral blood flow evaluation during cognitive tasks following a cervical spinal cord injury: a case study using transcranial Doppler recordings

A spinal cord injury (SCI) is one of the most common neurological disorders. In this paper, we examined the consequences of upper SCI in a male participant on the cerebral blood flow velocity. In particular, transcranial Doppler was used to study these effects through middle cerebral arteries (MCA) during resting-state periods and during cognitive challenges (non-verbal word-generation tasks and geometric-rotation tasks). Signal characteristics were analyzed from raw signals and envelope signals (maximum velocity) in the time domain, the frequency domain and the time-frequency domain. The frequency features highlighted an increase of the peak frequency in L-MCA and R-MCA raw signals, which revealed stronger cerebral blood flow during geometric/verbal processes respectively. This underlined a slight dominance of the right hemisphere during word-generation periods and a slight dominance of the left hemisphere during geometric processes. This finding was confirmed by cross-correlation in the time domain and by the entropy rate in information-theoretic domain. A comparison of our results to other neurological disorders (Alzheimer's disease, Parkinson's disease, autism, epilepsy, traumatic brain injury) showed that the SCI had similar effects such as general decreased cerebral blood flow and similar regular hemispheric dominance in a few cases.

Swallowing accelerometry signal feature variations with sensor displacement

Dual-axis accelerometry has recently shown promise as a non-invasive method for detecting swallowing impairment using signal processing and pattern classification algorithms. However, it is unknown whether variations in sensor placement alter signal characteristics, threatening the accuracy of signal processing classifiers for aspiration detection. To address this question, water swallows were recorded in 14 healthy adults using a dual-axis accelerometer in 13 different positions (baseline, and 2, 4, 6 and 8 mm above, below and to the right of baseline). The baseline position was midline, immediately below the thyroid cartilage during quiet breathing. After segmentation and pre-processing, signal features were extracted in multiple domains (time, frequency, time-frequency). The effect of sensor position on signal feature distributions was examined with non-parametric statistical analysis. The analysis showed that the sensor could be displaced by as much as 4 mm inferior and lateral to the baseline position and by up to 6 mm above the baseline location without significantly altering time-frequency features. In other words, when considering the baseline position as the origin, the admissible region for sensor placement spans 10 mm in the superior-inferior axis and 8 mm in the medial-lateral direction. Results of this study suggest that time-frequency representations of accelerometry signals are most robust to sensor placement variations around the baseline position. The implication of this finding is that a swallowing accelerometry classifier based on time-frequency features can likely tolerate small variations in sensor location without degradation in classification performance.

Characteristics of Dry Chin-Tuck Swallowing Vibrations and Sounds

OBJECTIVE

The effects of the chin-tuck maneuver, a technique commonly employed to compensate for dysphagia, on cervical auscultation are not fully understood. Characterizing a technique that is known to affect swallowing function is an important step on the way to developing a new instrumentation-based swallowing screening tool.

METHODS

In this study, we recorded data from 55 adult participants who each completed five saliva swallows in a chin-tuck position. The resulting data were processed using previously designed filtering and segmentation algorithms. We then calculated nine time-, frequency-, and time-frequency-domain features for each independent signal.

RESULTS

We found that multiple frequency- and time-domain features varied significantly between male and female subjects as well as between swallowing sounds and vibrations. However, our analysis showed that participant age did not play a significant role on the values of the extracted features. Finally, we found that various frequency features corresponding to swallowing vibrations did demonstrate statistically significant variation between the neutral and chin-tuck positions but sounds showed no changes between these two positions.

CONCLUSION

The chin-tuck maneuver affects many facets of swallowing vibrations and sounds and its effects can be monitored via cervical auscultation.

SIGNIFICANCE

These results suggest that a subject's swallowing technique does need to be accounted for when monitoring their performance with cervical auscultation-based instrumentation.

A comparative analysis of swallowing accelerometry and sounds during saliva swallows

BACKGROUND

Accelerometry (the measurement of vibrations) and auscultation (the measurement of sounds) are both non-invasive techniques that have been explored for their potential to detect abnormalities in swallowing. The differences between these techniques and the information they capture about swallowing have not previously been explored in a direct comparison.

METHODS

In this study, we investigated the differences between dual-axis swallowing accelerometry and swallowing sounds by recording data from adult participants and calculating a number of time and frequency domain features. During the experiment, 55 participants (ages 18-65) were asked to complete five saliva swallows in a neutral head position. The resulting data was processed using previously designed techniques including wavelet denoising, spline filtering, and fuzzy means segmentation. The pre-processed signals were then used to calculate 9 time, frequency, and time-frequency domain features for each independent signal. Wilcoxon signed-rank and Wilcoxon rank-sum tests were utilized to compare feature values across transducers and patient demographics, respectively.

RESULTS

In addition to finding a number of features that varied between male and female participants, our statistical analysis determined that the majority of our chosen features were statistically significantly different across the two sensor methods and that the dependence on within-subject factors varied with the transducer type. However, a regression analysis showed that age accounted for an insignificant amount of variation in our signals.

CONCLUSIONS

We conclude that swallowing accelerometry and swallowing sounds provide different information about deglutition despite utilizing similar transduction methods. This contradicts past assumptions in the field and necessitates the development of separate analysis and processing techniques for swallowing sounds and vibrations.

An SVM-based classifier for estimating the state of various rotating components in agro-industrial machinery with a vibration signal acquired from a single point on the machine chassis

The goal of this article is to assess the feasibility of estimating the state of various rotating components in agro-industrial machinery by employing just one vibration signal acquired from a single point on the machine chassis. To do so, a Support Vector Machine (SVM)-based system is employed. Experimental tests evaluated this system by acquiring vibration data from a single point of an agricultural harvester, while varying several of its working conditions. The whole process included two major steps. Initially, the vibration data were preprocessed through twelve feature extraction algorithms, after which the Exhaustive Search method selected the most suitable features. Secondly, the SVM-based system accuracy was evaluated by using Leave-One-Out cross-validation, with the selected features as the input data. The results of this study provide evidence that (i) accurate estimation of the status of various rotating components in agro-industrial machinery is possible by processing the vibration signal acquired from a single point on the machine structure; (ii) the vibration signal can be acquired with a uniaxial accelerometer, the orientation of which does not significantly affect the classification accuracy; and, (iii) when using an SVM classifier, an 85% mean cross-validation accuracy can be reached, which only requires a maximum of seven features as its input, and no significant improvements are noted between the use of either nonlinear or linear kernels.

Nonlinear analysis of saccade speed fluctuations during combined action and perception tasks

BACKGROUND

Saccades are rapid eye movements used to gather information about a scene which requires both action and perception. These are usually studied separately, so that how perception influences action is not well understood. In a dual task, where the subject looks at a target and reports a decision, subtle changes in the saccades might be caused by action-perception interactions. Studying saccades might provide insight into how brain pathways for action and for perception interact.

NEW METHOD

We applied two complementary methods, multifractal detrended fluctuation analysis and Lempel-Ziv complexity index to eye peak speed recorded in two experiments, a pure action task and a combined action-perception task.

RESULTS

Multifractality strength is significantly different in the two experiments, showing smaller values for dual decision task saccades compared to simple-task saccades. The normalized Lempel-Ziv complexity index behaves similarly i.e. is significantly smaller in the decision saccade task than in the simple task.

COMPARISON WITH EXISTING METHODS

Compared to the usual statistical and linear approaches, these analyses emphasize the character of the dynamics involved in the fluctuations and offer a sensitive tool for quantitative evaluation of the multifractal features and of the complexity measure in the saccades peak speeds when different brain circuits are involved.

CONCLUSION

Our results prove that the peak speed fluctuations have multifractal characteristics with lower magnitude for the multifractality strength and for the complexity index when two neural pathways are simultaneously activated, demonstrating the nonlinear interaction in the brain pathways for action and perception.

Multiscale Lempel-Ziv complexity for EEG measures

OBJECTIVE

To demonstrate that the classical calculation of Lempel-Ziv complexity (LZC) has an important limitation when applied to EEGs with rapid rhythms, and to propose a multiscale approach that overcomes this limitation.

METHODS

We have evaluated, both with simulated and real EEGs, whether LZC calculation neglects functional characteristics of rapid EEG rhythms. In addition, we have proposed a procedure to obtain multiple binarization sequences that yield a spectrum of LZC, and we have explored whether complexity would be better captured using this computation.

RESULTS

In our simulated signals, classical LZC did not capture modulations of a rapid component when a slower component of more amplitude was included in the signal. In real EEGs from healthy participants with eyes closed and eyes open, classical LZC calculation failed to show any difference between these two conditions. However, a multiscale LZC showed that complexity was lower for eyes closed than for eyes open conditions.

CONCLUSIONS

As hypothesized, our new approximation captures the complexity of series with fast components masked by slower rhythms.

SIGNIFICANCE

The method we introduce significantly improves LZC calculation, and it allows a better characterization of complexity of EEG signals.

60 Hz electric field changes the membrane potential during burst phase in pancreatic β-cells: in silico analysis

The production, distribution and use of electricity can generate low frequency electric and magnetic fields (50-60 Hz). Considering that some studies showed adverse effects on pancreatic β-cells exposed to these fields; the present study aimed to analyze the effects of 60 Hz electric fields on membrane potential during the silent and burst phases in pancreatic β-cells using a mathematical model. Sinusoidal 60 Hz electric fields with amplitude ranging from 0.5 to 4 mV were applied on pancreatic β-cells model. The sinusoidal electric field changed burst duration, inter-burst intervals (silent phase) and spike sizes. The parameters above presented dose-dependent response with the voltage amplitude applied. In conclusion, theoretical analyses showed that a 60 Hz electric field with low amplitudes changes the membrane potential in pancreatic β-cells.

Networks in financial markets based on the mutual information rate

In the last few years there have been many efforts in econophysics studying how network theory can facilitate understanding of complex financial markets. These efforts consist mainly of the study of correlation-based hierarchical networks. This is somewhat surprising as the underlying assumptions of research looking at financial markets are that they are complex systems and thus behave in a nonlinear manner, which is confirmed by numerous studies, making the use of correlations which are inherently dealing with linear dependencies only baffling. In this paper we introduce a way to incorporate nonlinear dynamics and dependencies into hierarchical networks to study financial markets using mutual information and its dynamical extension: the mutual information rate. We show that this approach leads to different results than the correlation-based approach used in most studies, on the basis of 91 companies listed on the New York Stock Exchange 100 between 2003 and 2013, using minimal spanning trees and planar maximally filtered graphs.

Comparative assessment of nonlinear metrics to quantify organization-related events in surface electrocardiograms of atrial fibrillation

Atrial fibrillation (AF) is today the most common sustained arrhythmia, its treatment being not completely satisfactory. Electrical activity organization analysis within the atria could play a key role in the improvement of current AF therapies. The application of a nonlinear regularity index, such as sample entropy (SampEn), to the atrial activity (AA) fundamental waveform has proven to be a successful noninvasive AF organization estimator. However, the use of alternative nonlinear metrics within this context is a pending issue. The present work analyzes the ability of several nonlinear indices to assess regularity of patterns and, thus, organization, in the AA signal and its fundamental waveform, defined as the main atrial wave (MAW). Precisely, Fuzzy Entropy, Spectral Entropy, Lempel-Ziv Complexity and Hurst Exponents were studied, achieving more robust and accurate AF organization estimates than SampEn. Results also provided better AF organization estimates from the MAW than from the AA signal for all the tested nonlinear metrics, which agrees with previous works only focused on SampEn. Furthermore, some of these indices reported a discriminant ability close to 95% in the classification of AF organization-dependent events, thus outperforming the diagnostic accuracy of SampEn and other widely used noninvasive estimators, such as the dominant atrial frequency (DAF). As a conclusion, these nonlinear metrics could be considered as promising estimators of noninvasive AF organization and could be helpful in making appropriate decisions on the patients' management.