Robust, real-time event detection from physiological signals acquired during long-term ambulatory monitoring still represents a major challenge for highly-artifacted signals. In this paper, we ...propose an original and generic multi-feature probabilistic detector (MFPD) and apply it to real-time QRS complex detection under noisy conditions. The MFPD method calculates a binary Bayesian probability for each derived feature and makes a centralized fusion, using the Kullback-Leibler divergence. The method is evaluated on two ECG databases: 1) the MIT-BIH arrhythmia database from Physionet containing clean ECG signals, 2) a benchmark noisy database created by adding noise recordings of the MIT-BIH noise stress test database, also from Physionet, to the MIT-BIH arrhythmia database. Results are compared with a well-known wavelet-based detector, and two recently published detectors: one based on spatiotemporal characteristic of the QRS complex and the second, as the MFDP, based on feature calculations from the University of New South Wales detector (UNSW). For both benchmark Physionet databases, the proposed MFPD method achieves the lowest standard deviation in sensitivity and positive predictivity (+P) despite its online algorithm architecture. While the statistics are comparable for low-to mildly artifactual ECG signals, the MFPD outperforms reference methods for artifacted ECG with low SNR levels reaching 87.48 ± 14.21% in sensitivity and 89.39 ± 14.67% in +P as compared to 88.30 ± 17.66% and 86.06 ± 19.67% respectively from UNSW, the best performing reference method. With demonstrations on the extensively studied QRS detection problem, we consider that the proposed generic structure of the multi-feature probabilistic detector should offer promising perspectives for long-term monitoring applications for highly-artifacted signals.
Heart rate variability (HRV) has been emerging in neonatal medicine. It may help for the early diagnosis of pathology and estimation of autonomous maturation. There is a lack of standardization and ...automation in the selection of the sequences to analyze and some features have not been explored in this specific population. The main objective of this study was to analyze the impact of the time length of the sequences on the estimation of linear and non-linear HRV features, including horizontal visibility graphs (HVG).
HRV features were repeatedly measured with linear and non-linear methods on 2-, 5-, 10-minute sequences selected from the longest 15-min sequence and recorded on a weekly basis in 39 infants less than 31 weeks at birth. The associations between HRV measurements were analyzed through principal component analysis and k-means clustering. The effects of the time lengths on HRV measurements and post-menstrual age (PMA) were analyzed by linear mixed effect model for repeated measures.
The domains of analysis were concordant for their descriptive parameters of short (rMSSD, SD1 and HF) and long-term (SD, SD2 and LF) variability. α1 was correlated with the LF/HF and SD2/SD1. DC and AC were correlated with short-term variability estimates and significantly increased with GA and PMA. Shortening the windows of analysis increased the random measurement error for all the features and increased the bias for all but short term features and HVGs.
The linear and non-linear measurements of HRV are correlated each other. Shortening the windows of analysis increased the random error for all the features and increased the bias for all but short term features and HVGs. Short-term HRV can be an index for evaluating the maturation in whatever sequence length.
Cardiac resynchronization therapy (CRT) is an implant-based therapy applied to patients with a specific heart failure (HF) profile. The identification of patients that may benefit from CRT is a ...challenging task and the application of current guidelines still induce a non-responder rate of about 30%. Several studies have shown that the assessment of left ventricular (LV) mechanics by speckle tracking echocardiography can provide useful information for CRT patient selection. A comprehensive evaluation of LV mechanics is normally performed using three different echocardioraphic views: 4, 3 or 2-chamber views. The aim of this study is to estimate the relative importance of strain-based features extracted from these three views, for the estimation of CRT response. Several features were extracted from the longitudinal strain curves of 130 patients and different methods of feature selection (out-of-bag random forest, wrapping and filtering) have been applied. Results show that more than 50% of the 20 most important features are calculated from the 4-chamber view. Although features from the 2- and 3-chamber views are less represented in the most important features, some of the former have been identified to provide complementary information. A thorough analysis and interpretation of the most informative features is also provided, as a first step towards the construction of a machine-learning chain for an improved selection of CRT candidates.
Although the therapeutic effects of Vagus Nerve Stimulation (VNS) have been recognized in pre-clinical and pilot clinical studies, the effect of different stimulation configurations on the ...cardiovascular response is still an open question, especially in the case of VNS delivered synchronously with cardiac activity. In this paper, we propose a formal mathematical methodology to analyze the acute cardiac response to different VNS configurations, jointly considering the chronotropic, dromotropic and inotropic cardiac effects. A latin hypercube sampling method was chosen to design a uniform experimental plan, composed of 75 different VNS configurations, with different values for the main parameters (current amplitude, number of delivered pulses, pulse width, interpulse period and the delay between the detected cardiac event and VNS onset). These VNS configurations were applied to 6 healthy, anesthetized sheep, while acquiring the associated cardiovascular response. Unobserved VNS configurations were estimated using a Gaussian process regression (GPR) model. In order to quantitatively analyze the effect of each parameter and their combinations on the cardiac response, the Sobol sensitivity method was applied to the obtained GPR model and inter-individual sensitivity markers were estimated using a bootstrap approach. Results highlight the dominant effect of pulse current, pulse width and number of pulses, which explain respectively 49.4%, 19.7% and 6.0% of the mean global cardiovascular variability provoked by VNS. More interestingly, results also quantify the effect of the interactions between VNS parameters. In particular, the interactions between current and pulse width provoke higher cardiac effects than the changes on the number of pulses alone (between 6 and 25% of the variability). Although the sensitivity of individual VNS parameters seems similar for chronotropic, dromotropic and inotropic responses, the interacting effects of VNS parameters provoke significantly different cardiac responses, showing the feasibility of a parameter-based functional selectivity. These results are of primary importance for the optimal, subject-specific definition of VNS parameters for a given therapy and may lead to new closed-loop methods allowing for the optimal adaptation of VNS therapy through time.
This paper proposes a model-based estimation of left ventricular (LV) pressure for the evaluation of constructive and wasted myocardial work of patients with aortic stenosis (AS). A model of the ...cardiovascular system is proposed, including descriptions of i) cardiac electrical activity, ii) elastance-based cardiac cavities, iii) systemic and pulmonary circulations and iv) heart valves. After a sensitivity analysis of model parameters, an identification strategy was implemented using a Monte-Carlo cross-validation approach. Parameter identification procedure consists in two steps for the estimation of LV pressures: step 1) from invasive, intraventricular measurements and step 2) from non-invasive data. The proposed approach was validated on data obtained from 12 patients with AS. The total relative errors between estimated and measured pressures were on average 11.9% and 12.27% and mean R2 were equal to 0.96 and 0.91, respectively for steps 1 and 2 of parameter identification strategy. Using LV pressures obtained from non-invasive measurements (step 2) and patient-specific simulations, Global Constructive (GCW), Wasted (GWW) myocardial Work and Global Work Efficiency (GWE) parameters were calculated. Correlations between measures and model-based estimations were 0.88, 0.80, 0.91 respectively for GCW, GWW and GWE. The main contributions concern the proposal of the parameter identification procedure, applied on an integrated cardiovascular model, able to reproduce LV pressure specifically to each AS patient, by non-invasive procedures, as well as a new method for the non-invasive estimation of constructive, wasted myocardial work and work efficiency in AS.
Cardiac resynchronization therapy (CRT) has been shown to improve cardiovascular function in specific patients suffering from heart failure. This procedure still needs to be optimized to overcome the ...high rate of implanted patients that do not respond to this therapy. We propose in this work a better characterization of the electro-mechanical (EM) coupling of each region of the left ventricle (LV) that could be useful to precise the best implantation site. A new descriptor is proposed with the extraction of local electro-mechanical delays. Their measurement is based on the fusion of anatomical, functional and electrical data acquired using computed tomography (CT), speckle tracking echocardiography (STE), and electro-anatomical mappings (EAM). We propose a workflow to place multimodal data in the same geometrical referential system and to extract local electro-mechanical descriptors. It implies the fusion of electrical and mechanical data on a 3D+ t anatomical model of the LV. It mainly consists in four steps: 1) the modeling of the endocardium using a dynamic surface estimated from CT images; 2) the semi-interactive registration of EAM data and CT images; 3) the automatic registration of STE data on the dynamic model, using a metric based on Fourier descriptors and dynamic time warping; 4) the temporal alignment between EAM and STE and the estimation of local electro-mechanical delays. The proposed process has been applied to real data corresponding to five patients undergoing CRT. Results show that local electro-mechanical delays provide meaningful information on the local characterization of the LV and may be useful for the optimal pacing site selection in CRT.
Brugada syndrome (BS) is a genetic pathological condition associated with a high risk for sudden cardiac death (SCD). Ventricular depolarization disorders have been suggested as a potential ...electrophysiological mechanism associated with high SCD risk on patients with BS. This paper aims to characterize the dynamic changes of ventricular depolarization observed during physical exercise in symptomatic and asymptomatic BS patients. To this end, cardiac ventricular depolarization features were automatically extracted from 12-lead ECG recordings acquired during standardized exercise stress test in 110 BS patients, of whom 25 were symptomatic. Conventional parameters were evaluated, including QRS duration, R and S wave amplitudes (Formula: see text, Formula: see text), as well as QRS morphological features, such as up-stroke and down-stroke slopes of the R and S waves (Formula: see text, Formula: see text and Formula: see text). The effects of physical exercise and recovery on the dynamics of these markers were assessed in both BS populations. Features showing significantly different dynamics between the studied groups were used alone and in combination with the clinical characteristics of the patients in a logistic regression analysis. Results show larger changes in the second half of the QRS complex through Formula: see text and Formula: see text measured in the right precordial leads for asymptomatic patients, especially during recovery, when the vagal tone is more pronounced. Multivariate analysis involving both types of features resulted in a reduced model of three relevant features (Formula: see text in lead V2, Sex and heart rate recovery, HRR), which achieved a suitable discrimination performance between groups; sensitivity = 80% and specificity = 75% (AUC = 83%). However, after controlling the model for possible confounding factors, only one feature (Formula: see text) remained meaningful. This adjusted model significantly improved the overall discrimination performance by up to: sensitivity = 84% and specificity = 100% (AUC = 94%). The study highlights the importance of physical exercise test to unmask differentiated behaviors between symptomatic and asymptomatic BS patients through depolarization dynamic analysis. This analysis together with the obtained model may help to identify asymptomatic patients at low or high risk of future cardiac events, but it should be confirmed by further prospective studies.
Vagus nerve stimulation (VNS) is an established adjunctive therapy for pharmacologically refractory epilepsy and depression and is currently in active clinical research for other applications. In ...current clinical studies, VNS is delivered in an open-loop approach, where VNS parameters are defined during a manual titration phase. However, the physiological response to a given VNS configuration shows significant inter and intra-patient variability and may significantly evolve through time. VNS closed-loop approaches, allowing for the optimization of the therapy in an adaptive manner, may be necessary to improve efficacy while reducing side effects. This paper proposes a generic, closed-loop control VNS system that is able to optimize a number of VNS parameters in an adaptive fashion, in order to keep a control variable within a specified range. Although the proposed control method is completely generic, an example application using the cardiac beat to beat interval (RR) as control variable will be developed in this paper. The proposed controller is based on a state transition model (STM) that can be configured using a partially or a fully-connected architecture, different model orders and different state-transition algorithms. The controller is applied to the adaptive regulation of heart rate and evaluated on 6 sheep, for 13 different targets, using partially-connected STM with 10 states. Also, partially and fully-connected STM defined by 30 states were applied to 7 other sheep for the same 10 targets. Results illustrate the interest of the proposed fully-connected STM and the feasibility of integrating this control system into an implantable neuromodulator.
As in other complex signal processing tasks, ECG beat delineation algorithms are usually constituted of a set of processing modules, each one characterized by a certain number of parameters (filter ...cutoff frequencies, threshold levels, time windows, etc.). It is well recognized that the adjustment of these parameters is a complex task that is traditionally performed empirically and manually, based on the experience of the designer. In this paper, we propose a new automated and quantitative method to optimize the parameters of such complex signal processing algorithms. To solve this multiobjective optimization problem, an evolutionary algorithm (EA) is proposed. This method for parameter optimization is applied to a wavelet-transform-based ECG delineator that has previously shown interesting performance. An evaluation of the final delineator, using the optimal parameters, has been performed on the QT database from Physionet and results are compared with previous algorithms reported in the literature. The optimized parameters provide a more accurate delineation, with a global improvement of 7.7%, over all the criteria evaluated, and over the best results found in the literature, which is a proof of interest in the approach.
Background and Aim. Variceal bleeding is the second most important precipitating factor related to the development of episodic hepatic encephalopathy; but to date there are no recommendations to ...prevent this complication. The aim of this study was to compare if primary prophylaxis with lactulose or L-ornithine L-aspartate or rifaximin, in cirrhotic patients with variceal bleeding, is better than placebo for avoiding the development of hepatic encephalopathy. Methods. A randomized, double-blind, placebo-controlled clinical trial (ClinicalTrials.gov identifier: NCT02158182) which included cirrhotic patients with variceal bleeding, without minimal or clinical hepatic encephalopathy at admission. Findings. 87 patients were randomized to one of four groups. The basal characteristics were similar between groups. Comparatively with placebo, the frequency with regard to the development of hepatic encephalopathy was as follows: lactulose (54.5% versus 27.3%; OR = 0.3, 95% CI 0.09-1.0; P = 0.06); L-ornithine L-aspartate (54.5% versus 22.7%, OR = 0.2, 95% CI 0.06-0.88; P = 0.03); rifaximin (54.5% versus 23.8%; OR = 0.3, 95% CI 0.07-0.9; P = 0.04). There was no significant difference between the three groups receiving any antiammonium drug (P = 0.94). In the group receiving lactulose, 59.1% had diarrhea, and 45.5% had abdominal discomfort, bloating, and flatulence. Two patients (10%) treated with lactulose and a patient (4.5%) in the placebo group developed spontaneous bacterial peritonitis due to E. coli; one of them died due to recurrent variceal bleeding. There were no other adverse effects. Conclusions. Antiammonium drugs, particularly L-ornithine L-aspartate and rifaximin, proved to be effective in preventing the development of hepatic encephalopathy in those cirrhotic patients with variceal bleeding.