Wavelet bidomain sample entropy analysis to predict spontaneous termination of atrial fibrillation

A NOVEL APPLICATION OF HIGHER ORDER STATISTICS OF R-R INTERVAL SIGNAL IN EMD DOMAIN FOR PREDICTING TERMINATION OF ATRIAL FIBRILLATION

Predicting termination of atrial fibrillation (AF), based on noninvasive techniques, can be invaluable in order to avoid useless therapeutic interventions and to minimize the risks for the patients. Currently, no reliable method exists to predict the termination of AF. We propose an algorithm for predicting termination of AF using higher order statistical moments of R-R interval signal calculated in both time and empirical mode decomposition (EMD) domains. In the proposed method, R-R interval signal is decomposed into a set of intrinsic mode functions (IMF) and higher order moments including skewness, and kurtosis, as well as mean and variance, are calculated from the first four IMFs. The appropriateness of these features in predicting the termination of AF is studied using atrial fibrillation termination database (AFTDB) which consists of three types of AF episodes: N-type (non-terminated AF episode), S-type (terminated 1'min after the end of the record), and T-type (terminated immediately after the end of the record). By using a support vector machine (SVM) classifier for classification of AF episodes, we obtained sensitivity, specificity, and positive predictivity 92.47%, 95.29%, and 92.80%, respectively. The important advantage of the proposed method compared to the other existing approaches is that our algorithm can simultaneously discriminate the three types of AF episodes with high accuracy. The results demonstrate that the EMD domain is a particularly well-suited domain for analyzing nonstationary and nonlinear R-R interval signal in AF termination prediction application.

AUTOMATED DETECTION OF ATRIAL FLUTTER AND FIBRILLATION USING ECG SIGNALS IN WAVELET FRAMEWORK

In this paper, an electrocardiogram (ECG)-based pattern analysis methodology is presented for the detection of artrial flutter and atrial fibrillation using fractal dimension (FD) of continuous wavelet transform (CWT) coefficients of raw ECG signals, sample entropy of heart beat interval time series, and mean heart beat interval features. Accurate diagnosis of atrial tachyarrhythmias is important, as they have different therapeutic options and clinical decisions. In view of this, we have made an attempt to develop a discrimination mechanism between artrial flutter and atrial fibrillation. The methodology consists of mean heart beat interval detection using Pan Tompkins algorithm, calculation of sample entropy of heart beat interval time series, computation of box counting FD from CWT coefficients of raw ECG, statistical significance test, and subsequent pattern classification using different classifiers. Different wavelet basis functions like Daubechies-4, Daubechies-6, Symlet-2, Symlet-4, Symlet-6, Symlet-8, Coiflet-2, Coiflet-5, Biorthogonal-1.3, Biorthogonal-3.1, and Mayer wavelet have been used to compute CWT coefficients. Features are evaluated using statistical analysis and subsequently two-class pattern classification is done using unsupervised (k-means, fuzzy c-means, and Gaussian mixture model) and supervised (error back propagation neural network and support vector machine) techniques. In order to reduce the bias in choosing the training and testing set, k-fold cross validation is used. The obtained results are compared and discussed. It is found that the supervised classifiers provide higher accuracy in comparison to the set of unsupervised classifiers.

Application of Wavelet Entropy to predict atrial fibrillation progression from the surface ECG

Atrial fibrillation (AF) is the most common supraventricular arrhythmia in clinical practice, thus, being the subject of intensive research both in medicine and engineering. Wavelet Entropy (WE) is a measure of the disorder degree of a specific phenomena in both time and frequency domains, allowing to reveal underlying dynamical processes out of sight for other methods. The present work introduces two different WE applications to the electrocardiogram (ECG) of patients in AF. The first application predicts the spontaneous termination of paroxysmal AF (PAF), whereas the second one deals with the electrical cardioversion (ECV) outcome in persistent AF patients. In both applications, WE was used with the objective of assessing the atrial fibrillatory (f) waves organization. Structural changes into the f waves reflect the atrial activity organization variation, and this fact can be used to predict AF progression. To this respect, results in the prediction of PAF termination regarding sensitivity, specificity, and accuracy were 95.38%, 91.67%, and 93.60%, respectively. On the other hand, for ECV outcome prediction, 85.24% sensitivity, 81.82% specificity, and 84.05% accuracy were obtained. These results turn WE as the highest single predictor of spontaneous PAF termination and ECV outcome, thus being a promising tool to characterize non-invasive AF signals.

Central tendency measure and wavelet transform combined in the non-invasive analysis of atrial fibrillation recordings

Background

Atrial fibrillation (AF) is the most common supraventricular arrhythmia in the clinical practice, being the subject of intensive research.

Methods

The present work introduces two different Wavelet Transform (WT) applications to electrocardiogram (ECG) recordings of patients in AF. The first one predicts spontaneous termination of paroxysmal AF (PAF), whereas the second one deals with the prediction of electrical cardioversion (ECV) outcome in persistent AF patients. In both cases, the central tendency measure (CTM) from the first differences scatter plot was applied to the AF wavelet decomposition. In this way, the wavelet coefficients vector CTM associated to the AF frequency scale was used to assess how atrial fibrillatory (f) waves variability can be related to AF events.

Results

Structural changes into the f waves can be assessed by combining WT and CTM to reflect atrial activity organization variation. This fact can be used to predict organization-related events in AF. To this respect, results in the prediction of PAF termination regarding sensitivity, specificity and accuracy were 100%, 91.67% and 96%, respectively. On the other hand, for ECV outcome prediction, 82.93% sensitivity, 90.91% specificity and 85.71% accuracy were obtained. Hence, CTM has reached the highest diagnostic ability as a single predictor published to date.

Conclusions

Results suggest that CTM can be considered as a promising tool to characterize non-invasive AF signals. In this sense, therapeutic interventions for the treatment of paroxysmal and persistent AF patients could be improved, thus, avoiding useless procedures and minimizing risks.

Noninvasive time and frequency predictors of long-standing atrial fibrillation early recurrence after electrical cardioversion

BACKGROUND

Several clinical factors have been studied to predict atrial fibrillation (AF) recurrence after electrical cardioversion (ECV) with limited predictive value.

METHODS

 A method able to predict robustly long-standing AF early recurrence by characterizing noninvasively the electrical atrial activity (AA) with parameters related to its time course and spectral features is presented. To this respect, 63 patients (20 men and 43 women; mean age 73.4 ± 9.0 years; under antiarrhythmic drug treatment with amiodarone) who were referred for ECV of persistent AF were studied. During a 4-week follow-up, AF recurrence was observed in 41 patients (65.1%).

RESULTS

 RR variability and the studied AA spectral features, including dominant atrial frequency (DAF), its first harmonic and their amplitude, provided poor statistical differences between groups. On the contrary, f waves power (fWP) and Sample Entropy (SampEn) of the AA behaved as very good predictors. Patients who relapsed to AF presented lower fWP (0.036 ± 0.019 vs 0.081 ± 0.029 n.u.(2) , P < 0.001) and higher SampEn (0.107 ± 0.022 vs 0.086 ± 0.033, P < 0.01). Furthermore, fWP presented the highest predictive accuracy of 82.5%, whereas SampEn provided a 79.4%. The remaining features revealed accuracies lower than 70%. A stepwise discriminant analysis (SDA) provided a model based on fWP and SampEn with 90.5% of accuracy.

CONCLUSIONS

 The fWP has proved to predict long-standing AF early recurrence after ECV and can be combined with SampEn to improve its diagnostic ability. Furthermore, a thorough analysis of the results allowed outlining possible associations between these two features and the concomitant status of atrial remodeling.

Non-invasive atrial fibrillation organization follow-up under successive attempts of electrical cardioversion

The development of non-invasive tools able to provide valuable information about the effectiveness of a shock in external electrical cardioversion (ECV) is clinically relevant to enhance these protocols in the treatment of atrial fibrillation (AF). The present contribution analyzes the ability of a non-linear regularity index, such as sample entropy (SampEn), to follow-up non-invasively AF organization under successive attempts of ECV and to predict the effectiveness of every single shock. To this respect, the atrial activity (AA) preceding each delivered shock was extracted by using a QRST cancellation method. Next, the main atrial wave (MAW), which can be considered as the fundamental waveform associated to the AA, was obtained by applying a selective filtering centered on the dominant atrial frequency (DAF). Finally, the MAW organization was estimated with SampEn and two thresholds (Th1 = 0.1223 and Th2 = 0.0832) were established to predict the ECV outcome. Results indicated that, prior to the first attempt, all the patients who needed only one shock to restore NSR were below Th1. In addition, most of them were above Th2 in case of AF relapsing during the first month. Regarding several shocks, all the patients who maintained NSR more than one month were below Th2 after the first shock. Moreover, all the patients who relapsed to AF during the first month were between Th1 and Th2 and, finally, all the patients with ineffective ECV were above Th1. After each unsuccessful shock, a SampEn relative decrease was observed for the patients who finally reverted to NSR, but the largest variation took place after the first attempt, thus indicating that this shock plays the most important role in the procedure. Indeed, by considering jointly the patients who needed only one shock and the patients who needed several shocks, 91.67% (22 out of 24) of ECVs resulting in NSR, 93.55% (29 out of 31) of ECVs relapsing to AF during the first month and 100% (10 out of 10) of ECVs in which NSR was not restored were correctly classified. As conclusion, the MAW organization analysis via SampEn can provide useful information that could improve the effectiveness of conventional external ECV protocols used in AF treatment.

Surface ECG organization time course analysis along onward episodes of paroxysmal atrial fibrillation

The complete understanding of the mechanisms leading to the initiation, maintenance and self-termination of atrial fibrillation (AF) still is an unsolved challenge for cardiac electrophysiology. Studies in which AF has been induced have shown that electrophysiological and structural remodeling of the atria during the arrhythmia could play an important role in the transition from paroxysmal to persistent AF. However, to this day, the time course of the atrial remodeling along onward episodes of non-induced paroxysmal AF has not been investigated yet. In this work, a non-invasive method, based on the regularity estimation of AF through sample entropy (SampEn), has been used to assess the organization evolution along onward episodes of paroxysmal AF. Given that AF organization has been associated to the number of existing wavelets wandering throughout the atrial tissue, SampEn could be considered as a concomitant estimator of atrial remodeling. The achieved results, in close agreement with previous findings obtained from invasive recordings, proved several relevant aspects of arial remodeling. Firstly, a progressive disorganization increase (SampEn increase) along onward episodes of AF has been observed for 63% of the analyzed patients, whereas a stable AF organization degree has been appreciated in the remaining 37%. Next, a positive correlation between episode duration and SampEn has been obtained (R=0.541, p<0.01). Finally, a remarkable influence of the fibrillation-free interval, preceding each episode, on the corresponding level of AF organization at the onset of the subsequent AF episode has been observed, with a correlation between these two indices of R=0.389 (p<0.01). As a consequence, it could be considered that atrial electrophysiological dynamics that occur along onward paroxysmal AF episodes are reflected and can be quantified from ECG recordings through non-invasive organization estimation.

Optimal parameters study for sample entropy-based atrial fibrillation organization analysis

Sample entropy (SampEn) is a nonlinear regularity index that requires the a priori selection of three parameters: the length of the sequences to be compared, m, the patterns similarity tolerance, r, and the number of samples under analysis, N. Appropriate values for m, r and N have been recommended and widely used in the literature for the application of SampEn to some physiological time series, such as heart rate, hormonal data, etc. However, no guidelines exist for the selection of that values in other cases. Therefore, an optimal parameters study should be required for the application of SampEn to not previously analyzed biomedical signals. In the present work, a thorough analysis on the optimal values for m, r and N is presented within the context of atrial fibrillation (AF) organization estimation, computed from surface electrocardiogram recordings. Recently, the evaluation of AF organization through SampEn, has revealed clinically useful information that could be used for a better treatment of this arrhythmia. The present study analyzed optimal SampEn parameter values within two different scenarios of AF organization estimation, such as the prediction of paroxysmal AF termination and the electrical cardioversion outcome in persistent AF. As a result, interesting recommendations about the selection of m, r and N, together with the relationship between N and the sampling rate (f(s)) were obtained. More precisely, (i) the proportion between N and f(s) should be higher than 1s and f(s)>or=256 Hz, (ii) overlapping between adjacent N-length windows does not improve AF organization estimation with respect to the analysis of non-overlapping windows, and (iii) values of m and r maximizing successful classification for the analyzed AF databases should be considered within a range wider than the proposed in the literature for heart rate analysis, i.e. m=1 and m=2 and r between 0.1 and 0.25 times the standard deviation of the data.

Noninvasive assessment of the complexity and stationarity of the atrial wavefront patterns during atrial fibrillation

A novel automated approach to quantitatively evaluate the degree of spatio-temporal organization in the atrial activity (AA) during atrial fibrillation (AF) from surface recordings, obtained from body surface potential maps (BSPM), is presented. AA organization is assessed by measuring the reflection of the spatial complexity and temporal stationarity of the wavefront patterns propagating inside the atria on the surface ECG, by means of principal component analysis (PCA). Complexity and stationarity are quantified through novel parameters describing the structure of the mixing matrices derived by the PCA of the different AA segments across the BSPM recording. A significant inverse correlation between complexity and stationarity is highlighted by this analysis. The discriminatory power of the parameters in identifying different groups in the set of patients under study is also analyzed. The obtained results present analogies with earlier invasive studies in terms of number of significant components necessary to describe 95% of the variance in the AA (four for more organized AF, and eight for more disorganized AF). These findings suggest that automated analysis of AF organization exploiting spatial diversity in surface recordings is indeed possible, potentially leading to an improvement in clinical decision making and AF treatment.

Sample entropy analysis of cervical neoplasia gene-expression signatures

Background

We introduce Approximate Entropy as a mathematical method of analysis for microarray data. Approximate entropy is applied here as a method to classify the complex gene expression patterns resultant of a clinical sample set. Since Entropy is a measure of disorder in a system, we believe that by choosing genes which display minimum entropy in normal controls and maximum entropy in the cancerous sample set we will be able to distinguish those genes which display the greatest variability in the cancerous set. Here we describe a method of utilizing Approximate Sample Entropy (ApSE) analysis to identify genes of interest with the highest probability of producing an accurate, predictive, classification model from our data set.

Results

In the development of a diagnostic gene-expression profile for cervical intraepithelial neoplasia (CIN) and squamous cell carcinoma of the cervix, we identified 208 genes which are unchanging in all normal tissue samples, yet exhibit a random pattern indicative of the genetic instability and heterogeneity of malignant cells. This may be measured in terms of the ApSE when compared to normal tissue. We have validated 10 of these genes on 10 Normal and 20 cancer and CIN3 samples. We report that the predictive value of the sample entropy calculation for these 10 genes of interest is promising (75% sensitivity, 80% specificity for prediction of cervical cancer over CIN3).

Conclusion

The success of the Approximate Sample Entropy approach in discerning alterations in complexity from biological system with such relatively small sample set, and extracting biologically relevant genes of interest hold great promise.

Non-invasive organization variation assessment in the onset and termination of paroxysmal atrial fibrillation

Atrial Fibrillation (AF) is the most common supraventricular tachyarrhythmia. Recently, it has been suggested that AF is partially organized on its onset and termination, thus being more suitable for antiarrhythmia and to avoid unnecessary therapy. Although several invasive and non-invasive AF organization estimators have been proposed, the organization time course in the first and last minutes of AF has not been quantified yet. The aim of this work is to study non-invasively the organization variation within the first and last minutes of paroxysmal AF. The organization was evaluated making use of sample entropy, which can robustly estimate electrical atrial activity organization from surface ECG recordings. This work proves an organization decrease in the first minutes of AF onset and an increase within the last minute before spontaneous AF termination. These results are in agreement with the conclusions reported by other authors who made use of invasive recordings.

Analysis of epicardial mapping electrogram of sustained atrial fibrillation based on Shannon entropy

The mechanisms of sustained atrial fibrillation (AF) are not well understood. And predicting the development of AF is a problem of great clinical interest. This paper proposed an AF analysis method by evaluating, based on Shannon entropy, the complexity of atrial activation in various AF stages. The first step was to preprocess and characterize electrograms. Then, Shannon entropy analysis and statistical analysis were applied to find the significance of interval entropy in sustained AF. Study results proved that interval entropy presented a degressive tendency in the process of sustained AF and some sites with high activation frequency but low entropy was possibly related to ectopic driver of AF.

Adaptive singular value cancelation of ventricular activity in single-lead atrial fibrillation electrocardiograms

The proper analysis and characterization of atrial fibrillation (AF) from surface electrocardiographic (ECG) recordings requires to cancel out the ventricular activity (VA), which is composed of the QRS complex and the T wave. Historically, for single-lead ECGs, the averaged beat subtraction (ABS) has been the most widely used technique. However, this method is very sensitive to QRST wave variations and, moreover, high-quality cancelation templates may be difficult to obtain when only short length and single-lead recordings are available. In order to overcome these limitations, a new QRST cancelation method based on adaptive singular value cancelation (ASVC) applied to each single beat is proposed. In addition, an exhaustive study about the optimal set of complexes for better cancelation of every beat is also presented for the first time. The whole study has been carried out with both simulated and real AF signals. For simulated AF, the cancelation performance was evaluated making use of a cross-correlation index and the normalized mean square error (nmse) between the estimated and the original atrial activity (AA). For real AF signals, two additional new parameters were proposed. First, the ventricular residue (VR) index estimated the presence of ventricular activity in the extracted AA. Second, the similarity (S) evaluated how the algorithm preserved the AA segments out of the QRST interval. Results indicated that for simulated AF signals, mean correlation, nmse, VR and S values were 0.945 +/- 0.024, 0.332 +/- 0.073, 1.552 +/- 0.386 and 0.986 +/- 0.012, respectively, for the ASVC method and 0.866 +/- 0.042, 0.424 +/- 0.120, 2.161 +/- 0.564 and 0.922 +/- 0.051 for ABS. In the case of real signals, the mean VR and S values were 1.725 +/- 0.826 and 0.983 +/- 0.038, respectively, for ASVC and 3.159 +/- 1.097 and 0.951 +/- 0.049 for ABS. Thus, ASVC provides a more accurate beat-to-beat ventricular QRST representation than traditional techniques. As a consequence, VA cancelation is optimized and the AA can be extracted more precisely. Finally, the study has proven that optimal VA cancelation is achieved when a number between 20 and 30 complexes is selected following a correlation-based strategy.