Publications

Heart diseases are the leading cause of death in developed countries. Ascertaining the etiology of cardiomyopathies is still a challenge. The objective of this study was to classify cardiomyopathy patients through cardio, respiratory and vascular variability analysis, considering the vascular activity as the input and output of the baroreflex response. Forty-one cardiomyopathy patients (CMP) classified as ischemic (ICM, 24 patients) and dilated (DCM, 17 patients) were analyzed. Thirty-nine elderly control subjects (CON) were used as reference. From the electrocardiographic, respiratory flow, and blood pressure signals, following temporal series were extracted: beat-to-beat intervals (BBI), total respiratory cycle time series (TT), and end– systolic (SBP) and diastolic (DBP) blood pressure amplitudes, respectively. Three-dimensional representation of the cardiorespiratory and vascular activities was characterized geometrically, by fitting a polygon that contains 95% of data, and by statistical descriptive indices. The best classifiers were used to build support vector machine models. The optimal model to classify ICM versus DCM patients achieved 92.7% accuracy, 94.1% sensitivity, and 91.7% specificity. When comparing CMP patients and CON subjects, the best model achieved 86.2% accuracy, 82.9% sensitivity, and 89.7% specificity. These results suggest a limited ability of cardiac and respiratory systems response to regulate the vascular variability in these patients.

JTD Keywords: Time series analysis, Support vector machines, Blood pressure, Sensitivity, Indexes, Electrocardiography, Kernel

The intermittency of the instantaneous concentration of a turbulent chemical plume is a fundamental cue for estimating the chemical source distance using chemical sensors. Such estimate is useful in applications such as environmental monitoring or localization of fugitive gas emissions by mobile robots or sensor networks. However, the inherent low-pass filtering of metal oxide (MOX) gas sensors-typically used in odor-guided robots and dense sensor networks due to their low cost, weight and size-hinders the quantification of concentration intermittency. In this paper, we design a digital differentiator to invert the low-pass dynamics of the sensor response, thus obtaining a much faster signal from which the concentration intermittency can be effectively computed. Using a fast photo-ionization detector as a reference instrument, we demonstrate that the filtered signal is a good approximation of the instantaneous concentration in a real turbulent plume. We then extract transient features from the filtered signal-the so-called “bouts”-to predict the chemical source distance, focusing on the optimization of the filter parameters and the noise threshold to make the predictions robust against changing wind conditions. This represents an advantage over previous bout-based models which require wind measurements-typically taken with expensive and bulky anemometers-to produce accurate predictions. The proposed methodology is demonstrated in a wind tunnel scenario where a MOX sensor is placed at various distances downwind of an emitting chemical source and the wind speed varies in the range 10-34 cm/s. The results demonstrate that models optimized with our methodology can provide accurate source distance predictions at different wind speeds.

JTD Keywords: Gas detectors, Chemical sensors, Signal processing, Machine learning, Time series analysis

Cardiovascular diseases are one of the most common causes of death in elderly patients. The etiology of cardiomyopathies is difficult to discern clinically. The objective of this study was to classify cardiomyopathy patients using coupling analysis, through their cardiovascular behavior and the baroreflex response. A total of thirty-eight cardiomyopathy patients (CMP) classified as ischemic (ICM, 25 patients) and dilated (DCM, 13 patients) were analyzed. Thirty elderly control subjects (CON) were used as reference. Their electrocardiographic (ECG) and blood pressure (BP) signals were studied. To characterize the cardiovascular activity, the following temporal series were extracted: beat-to-beat intervals (from the ECG signal), and end- systolic and diastolic blood pressure amplitudes (from the BP signal). Non-linear characterization techniques like high resolution joint symbolic dynamics, segmented Poincaré plot analysis, normalized shorttime partial directed coherence, and dual sequence method were used to characterize these times series. The best indices were used to build support vector machine models for classification. The optimal model for ICM versus DCM patients achieved 84.2% accuracy, 76.9% sensitivity, and 88% specificity. When CMP patients and CON subjects were compared, the best model achieved 95.5% accuracy, 97.3% sensitivity, and 93.3% specificity. These results suggest a disfunction in the baroreflex mechanism in cardiomyopathies patients.

JTD Keywords: Couplings, Time series analysis, Support vector machines, Electrocardiography, Baroreflex, Coherence, Sensitivity

In this paper, a non-linear HRV analysis was performed to assess fragmentation signatures observed in heartbeat time series after intermittent hypoxia (IH). Three markers quantifying short-term fragmentation levels, PIP, IALS and PSS, were evaluated on R-R interval series obtained in a rat model of recurrent apnea. Through airway obstructions, apnea episodes were periodically simulated in six anesthetized Sprague-Dawley rats. The number of apnea events per hour (AHI index) was varied during the first half of the experiment while apnea episodes lasted 15 s. For the second part, apnea episodes lasted 5, 10 or 15 s, but the AHI index was fixed. Recurrent apnea was repeated for 15-min time intervals in all cases, alternating with basal periods of the same duration. The fragmentation markers were evaluated in segments of 5 minutes, selected at the beginning and end of the experiment. The impact of the heart and breathing rates (HR and BR, respectively) on the parameter estimates was also investigated. The results obtained show a significant increase (from 5 to 10%, p <; 0.05) in fragmentation measures of heartbeat time series after IH, indicating a clear deterioration of the initial conditions. Moreover, there was a strong linear relationship (r > 0.9) between these markers and BR, as well as with the ratio given by HR/BR. Although fragmentation may be impacted by IH, we found that it is highly dependent on HR and BR values and thus, they should be considered during its calculation or used to normalize the fragmentation estimates.

JTD Keywords: Rats, Time series analysis, Radio access technologies, Protocols, Heart beat

Cardiac death risk is still a big problem by an important part of the population, especially in elderly patients. In this study, we propose to characterize and analyze the cardiovascular and cardiorespiratory systems using the Poincaré plot. A total of 46 cardiomyopathy patients and 36 healthy subjets were analyzed. Left ventricular ejection fraction (LVEF) was used to stratify patients with low risk (LR: LVEF > 35%, 16 patients), and high risk (HR: LVEF ≤ 35%, 30 patients) of heart attack. RR, SBP and TTot time series were extracted from the ECG, blood pressure and respiratory flow signals, respectively. Parameters that describe the scatterplott of Poincaré method, related to short- and long-term variabilities, acceleration and deceleration of the dynamic system, and the complex correlation index were extracted. The linear discriminant analysis (LDA) and the support vector machines (SVM) classification methods were used to analyze the results of the extracted parameters. The results showed that cardiac parameters were the best to discriminate between HR and LR groups, especially the complex correlation index (p = 0.009). Analising the interaction, the best result was obtained with the relation between the difference of the standard deviation of the cardiac and respiratory system (p = 0.003). When comparing HR vs LR groups, the best classification was obtained applying SVM method, using an ANOVA kernel, with an accuracy of 98.12%. An accuracy of 97.01% was obtained by comparing patients versus healthy, with a SVM classifier and Laplacian kernel. The morphology of Poincaré plot introduces parameters that allow the characterization of the cardiorespiratory system dynamics.

JTD Keywords: Time series analysis, Electrocardiography, Support vector machines, Kernel, Standards, Correlation, RF signals

Synchronization is a central key feature of neural information processing and communication between different brain areas. Disturbance of oscillatory brain rhythms and decreased synchronization have been associated with different disorders including schizophrenia. The aim of this study was to investigate whether synchronization (in relaxed conditions with no stimuli) between different brain areas within the delta, theta, alpha (alpha1, alpha2), beta (beta1, beta2), and gamma bands is altered in patients with a neurological disorder in order to generate significant cortical enhancements. To achieve this, we investigated schizophrenic patients (SZO; N=17, 37.5±10.4 years, 15 males) and compared them to healthy subjects (CON; N=21, 36.7±13.4 years, 15 males) applying the phase locking value (PLV). We found significant differences between SZO and CON in different brain areas of the theta, alpha1, beta2 and gamma bands. These areas are related to the central and parietal lobes for the theta band, the parietal lobe for the alpha1, the parietal and frontal for the beta2 and the frontal-central for the gamma band. The gamma band revealed the most significant differences between CON and SZO. PLV were 61.7% higher on average in SZO in most of the clusters when compared to CON. The related brain areas are directly related to cognition skills which are proved to be impaired in SZO. The results of this study suggest that synchronization in SZO is also altered when the patients were not asked to perform a task that requires their cognitive skills (i.e., no stimuli are applied - in contrast to other findings).

JTD Keywords: Synchronization, Electroencephalography, Electrodes, Brain, Time series analysis, Oscillators, Frequency synchronization

In this work, the cardiorespiratory pattern of patients undergoing extubation process is studied. First, the respiratory and cardiac signals were resampled, next the Symbolic Dynamics (SD) technique was implemented, followed of a dimensionality reduction applying Forward Selection (FS) and Moving Window with Variance Analysis (MWVA) methods. Finally, the Linear Discriminant Analysis (LDA) and Support Vector Machines (SVM) classifiers were used. The study analyzed 153 patients undergoing weaning process, classified into 3 groups: Successful Group (SG: 94 patients), Failed Group (FG: 39 patients), and patients who had been successful during the extubation and had to be reintubated before 48 hours, Reintubated Group (RG: 21 patients). According to the results, the best classification present an accuracy higher than 88.98 ± 0.013% in all proposed combinations.

JTD Keywords: Support vector machines, Standards, Time series analysis, Resonant frequency, Nonlinear dynamical systems, Ventilation

Cardiovascular diseases are the first cause of death in developed countries. Using electrocardiographic (ECG), blood pressure (BP) and respiratory flow signals, we obtained parameters for classifying cardiomyophaty patients. 42 patients with ischemic (ICM) and dilated (DCM) cardiomyophaties were studied. The left ventricular ejection fraction (LVEF) was used to stratify patients with low risk (LR: LVEF>35%, 14 patients) and high risk (HR: LVEF≤ 35%, 28 patients) of heart attack. RR, SBP and T_Tot time series were extracted from the ECG, BP and respiratory flow signals, respectively. The time series were transformed to a binary space and then analyzed using Joint Symbolic Dynamic with a word length of three, characterizing them by the probability of occurrence of the words. Extracted parameters were then reduced using correlation and statistical analysis. Principal component analysis and support vector machines methods were applied to characterize the cardiorespiratory and cardiovascular interactions in ICM and DCM cardiomyopaties, obtaining an accuracy of 85.7%.

JTD Keywords: Blood pressure, Electrocardiography, Joints, Kernel, Principal component analysis, Support vector machines, Time series analysis

The scope of our work focuses on investigating the potential use of the built-in accelerometer of the smartphones for the recording of the respiratory activity and deriving the respiratory rate. Five healthy subjects performed an inspiratory load protocol. The excursion of the right chest was recorded using the built-in triaxial accelerometer of a smartphone along the x, y and z axes and with an external uniaxial accelerometer. Simultaneously, the respiratory airflow and the inspiratory mouth pressure were recorded, as reference respiratory signals. The chest acceleration signal recorded in the z axis with the smartphone was denoised using a scheme based on the ensemble empirical mode decomposition, a noise data assisted method which decomposes nonstationary and nonlinear signals into intrinsic mode functions. To distinguish noisy oscillatory modes from the relevant modes we use the detrended fluctuation analysis. We reported a very strong correlation between the acceleration of the z axis of the smartphone and the reference accelerometer across the inspiratory load protocol (from 0.80 to 0.97). Furthermore, the evaluation of the respiratory rate showed a very strong correlation (0.98). A good agreement was observed between the respiratory rate estimated with the chest acceleration signal from the z axis of the smartphone and with the respiratory airflow signal: Bland-Altman limits of agreement between -1.44 and 1.46 breaths per minute with a mean bias of -0.01 breaths per minute. This preliminary study provides a valuable insight into the use of the smartphone and its built-in accelerometer for respiratory monitoring.

JTD Keywords: Acceleration, Accelerometers, Correlation, Empirical mode decomposition, Fluctuations, Protocols, Time series analysis

Due to the increasing elderly population and the extensive number of comorbidities that affect them, studies are required to determine future increments in admission to emergency departments. Some of these studies could focus on the relation between chronic diseases and breathing pattern in elderly patients. Variations in the fractal properties of respiratory signals can be associated with several diseases. To determine the relationship between these variations and breathing patterns, and to quantify the fractal properties of respiratory flow signals, we estimated the Hurst exponent (H). Detrended fluctuation analysis (DFA) and discrete wavelet transform-based estimation (DWTE) methods were applied. The estimation methods were analyzed using simulated data series generated by fractional Gaussian noise. 43 elderly patients (19 patients with a non-periodic breathing pattern - nPB, and 24 patients with a periodic breathing pattern - PB) were studied. The results were evaluated according to the length of data and the number of averaged data series used to obtain a good estimation. The DWTE method estimated the respiratory flow signals better than the DFA method, and obtained Hurst values clustered by group. We found significant differences in the H exponent (p = 0.002) between PB and nPB patients, which showed different behavior in the fractal properties.

JTD Keywords: Discrete wavelet transforms, Diseases, Estimation, Fractals, Modulation, Senior citizens, Time series analysis

One of the most challenging problems in intensive care is still the process of discontinuing mechanical ventilation, called weaning process. Both an unnecessary delay in the discontinuation process and a weaning trial that is undertaken too early are undesirable. In this study, we analyzed respiratory pattern variability using the respiratory volume signal of patients submitted to two different levels of pressure support ventilation (PSV), prior to withdrawal of the mechanical ventilation. In order to characterize the respiratory pattern, we analyzed the following time series: inspiratory time, expiratory time, breath duration, tidal volume, fractional inspiratory time, mean inspiratory flow and rapid shallow breathing. Several autoregressive modeling techniques were considered: autoregressive models (AR), autoregressive moving average models (ARMA), and autoregressive models with exogenous input (ARX). The following classification methods were used: logistic regression (LR), linear discriminant analysis (LDA) and support vector machines (SVM). 20 patients on weaning trials from mechanical ventilation were analyzed. The patients, submitted to two different levels of PSV, were classified as low PSV and high PSV. The variability of the respiratory patterns of these patients were analyzed. The most relevant parameters were extracted using the classifiers methods. The best results were obtained with the interquartile range and the final prediction errors of AR, ARMA and ARX models. An accuracy of 95% (93% sensitivity and 90% specificity) was obtained when the interquartile range of the expiratory time and the breath duration time series were used a LDA model. All classifiers showed a good compromise between sensitivity and specificity.

JTD Keywords: autoregressive moving average processes, feature extraction, medical signal processing, patient care, pneumodynamics, signal classification, support vector machines, time series, ARX, autoregressive modeling techniques, autoregressive models with exogenous input, autoregressive moving average model, breath duration time series, classification method, classifier method, discontinuing mechanical ventilation, expiratory time, feature extraction, final prediction errors, fractional inspiratory time, intensive care, interquartile range, linear discriminant analysis, logistic regression analysis, mean inspiratory flow, patient respiratory volume signal, pressure support level, pressure support ventilation, rapid shallow breathing, respiratory pattern variability characterization, support vector machines, tidal volume, weaning trial, Analytical models, Autoregressive processes, Biological system modeling, Estimation, Support vector machines, Time series analysis, Ventilation

The process of weaning from mechanical ventilation is one of the challenges in intensive care units. 66 patients under extubation process (T-tube test) were studied: 33 patients with successful trials and 33 patients who failed to maintain spontaneous breathing and were reconnected. Each patient was characterized using 7 time series from respiratory signals, and for each serie was evaluated the discrete wavelet transform. It trains a neural network for discriminating between patients from the two groups.

JTD Keywords: discrete wavelet transforms, neural nets, patient treatment, pneumodynamics, time series, ventilation, T-tube test, discrete wavelet transform, extubation process, intensive care units, mechanical ventilation, moment of disconnection, neural network, patients, respiratory signals, spontaneous breathing, time series, weaning, Mechanical Ventilation, Neural Networks, Time series from respiratory signals, Wavelet Transform

In this paper, we present a new Multi-Agent System for reduction of the bullwhip effect in fuzzy supply chains. First, we show that a supply chain that uses an optimal ordering policy without data sharing among echelons still suffers from the bullwhip effect. Then, we propose the multi-agent solution to manage and reduce the bullwhip effect. The proposed multi-agent system includes four different types of agents in which each agent has its own list of actions. The proposed Multi-agent System applies a new Tabu Search algorithm for fuzzy rule generation, and a new data filtering algorithm for extraction of the bullwhip-free data from supply chain data warehouse. We validate the multi-agent system under different conditions and discuss how the system responds to different factors. The results show that the proposed multi-agent system reduces the bullwhip effect significantly in a rational time.

JTD Keywords: Bullwhip effect, Bullwhip-free data, Decentralized decision making, Fuzzy rule base, Fuzzy supply chain, Fuzzy time series, Multi-agent system, Supply chain management

Weaning trials process of patients in intensive care units is a complex clinical procedure. 153 patients under extubation process (T-tube test) were studied: 94 patients with successful trials (group S), 38 patients who failed to maintain spontaneous breathing and were reconnected (group F), and 21 patients with successful test but that had to be reintubated before 48 hours (group R). The respiratory pattern of each patient was characterized through the following time series: inspiratory time (T_I), expiratory time (T_E), breathing cycle duration (T_Tot), tidal volume (V_T), inspiratory fraction (T_I/T_Tot), half inspired flow (V_T/T_I), and rapid shallow index (f/V_T), where f is respiratory rate. Using techniques as autoregressive models (AR), autoregressive moving average models (ARMA) and autoregressive models with exogenous input (ARX), the most relevant parameters of the respiratory pattern were obtained. We proposed the evaluation of these parameters using classifiers as logistic regression (LR), linear discriminant analysis (LDA), support vector machines (SVM) and classification and regression tree (CART) to discriminate between patients from groups S, F and R. An accuracy of 93% (98% sensitivity and 82% specificity) has been obtained using CART classification.

JTD Keywords: Accuracy, Indexes, Logistics, Regression tree analysis, Support vector machines, Time series analysis, Autoregressive moving average processes, Medical signal processing, Pattern classification, Pneumodynamics, Regression analysis, Sensitivity, Signal classification, Support vector machines, Time series, SVM, T-tube testing, Autoregressive models-with-exogenous input, Autoregressive moving average models, Breathing cycle duration, Classification-and-regression tree, Expiratory time, Extubation process, Half inspired flow, Inspiratory fraction, Inspiratory time, Intensive care units, Linear discriminant analysis, Logistic regression, Rapid shallow index, Respiratory pattern parameter performance, Sensitivity, Spontaneous breathing, Support vector machines, Tidal volume, Time 48 hr, Time series, Weaning process classifiers