- Open Access
Atrial fibrillation detection by heart rate variability in Poincare plot
© Park et al; licensee BioMed Central Ltd. 2009
- Received: 17 September 2009
- Accepted: 11 December 2009
- Published: 11 December 2009
Atrial fibrillation (AFib) is one of the prominent causes of stroke, and its risk increases with age. We need to detect AFib correctly as early as possible to avoid medical disaster because it is likely to proceed into a more serious form in short time. If we can make a portable AFib monitoring system, it will be helpful to many old people because we cannot predict when a patient will have a spasm of AFib.
We analyzed heart beat variability from inter-beat intervals obtained by a wavelet-based detector. We made a Poincare plot using the inter-beat intervals. By analyzing the plot, we extracted three feature measures characterizing AFib and non-AFib: the number of clusters, mean stepping increment of inter-beat intervals, and dispersion of the points around a diagonal line in the plot. We divided distribution of the number of clusters into two and calculated mean value of the lower part by k-means clustering method. We classified data whose number of clusters is more than one and less than this mean value as non-AFib data. In the other case, we tried to discriminate AFib from non-AFib using support vector machine with the other feature measures: the mean stepping increment and dispersion of the points in the Poincare plot.
We found that Poincare plot from non-AFib data showed some pattern, while the plot from AFib data showed irregularly irregular shape. In case of non-AFib data, the definite pattern in the plot manifested itself with some limited number of clusters or closely packed one cluster. In case of AFib data, the number of clusters in the plot was one or too many. We evaluated the accuracy using leave-one-out cross-validation. Mean sensitivity and mean specificity were 91.4% and 92.9% respectively.
Because pulse beats of ventricles are less likely to be influenced by baseline wandering and noise, we used the inter-beat intervals to diagnose AFib. We visually displayed regularity of the inter-beat intervals by way of Poincare plot. We tried to design an automated algorithm which did not require any human intervention and any specific threshold, and could be installed in a portable AFib monitoring system.
- Support Vector Machine
- Heart Rate Variability
- Step Increment
- Radial Basis Function Kernel
- Support Vector Machine Method
There is a growing tendency that atrial fibrillation (AFib) related disease affects quality of life. Its risk increases with age ; in fact, AFib is one of the most common types of arrhythmia in clinical practice . Blood circulation of AFib patients is not smooth; therefore, AFib patients feel dizzy and uncomfortable while they exercise. AFib can be one of the deadliest symptoms to patients with preexcitation; in this case, it often induces tachycardia of ventricles or atrioventricular fibrillation . The more serious effect of AFib is formation of blood clots by congestion of blood in atria . If these blood clots come out of the atria and occlude a vessel somewhere in the brain, dire stroke can come about.
AFib can be classified into three grades: paroxysmal, persistent and permanent AFib. The paroxysmal AFib can be a preceding omen of the persistent AFib. Takahashi, Seki and Imatak observed that their patients with the paroxysmal AFib were highly affected with the more serious form of AFib; 25.3% of paroxysmal AFib patients developed into the more serious form of AFib in one year .
Electrical remodeling is one of the features of AFib and it is related to decreased conduction velocity of electricity signals . When the heart experiences the electrical remodeling, an area of slow conduction takes place in atria because of insufficient recovery of excitability. Slow conduction shortens wavelengths of the wandering electricity signals; thus, this area of slow conduction increases the number of re-current wave fronts of depolarization in atria and contributes to the sustaining of AFib . The wandering wave fronts around the atria fork themselves or collide with one another; accordingly, these maintain the turbulence process of electric conduction in atria . The re-entrant wave fronts induce inappropriate heart pumping; consequently, they deteriorate solidity of cardiac hemodynamics.
There are several studies about screening AFib by palpating an electrocardiogram (ECG) manually. Sudlow et al. tried to screen AFib by two methods: digoxin prescriptions and pulse palpation of ECG. Sensitivity and specificity using digoxin prescriptions were somewhat low. Sensitivity and specificity using pulse palpation were (93%, 71%) in case of women elder than 75, (100%, 86%) in case of 65-74 aged women, (95%, 71%) in case of men elder than 75, (100%, 79%) in case of 65-74 aged men, (sensitivity, specificity) respectively . Somerville et al. reported a screening result of AFib; Sensitivity and specificity were 100% and 77% respectively . Mant et al. showed the screening result by general practitioners and practice nurses observing 12 lead ECG. Sensitivity and specificity were 79.8% and 91.6% by general practitioners, 77.1% and 85.1% by practice nurses .
There are lots of studies about detecting AFib. Xu et al. chose five feature parameters which were input regularity, input atrial rate, energy distribution, time interval corresponding to zero amplitude signal, and number of points reaching baseline. They used Bayesian discriminator to classify the input data as one of sinus rhythm, AFib or atrial flutter . Petrucci et al. used two histograms which were calculated from the inter-beat intervals. One histogram consisted of differences between two successive inter-beat intervals and the other histogram consisted of normalized deviations from mean value of the inter-beat intervals. They calculated distribution widths from these histograms to discriminate AFib from non-AFib . Kikillus et al. made a Poincaré plot from inter-beat intervals and estimated density of points in each segment of Poincaré plot. They calculated an indicator of AFib from standard deviation of temporal differences of the consecutive inter-beat intervals . Thuraisingham used wavelet method to obtain a filtered time series from the input ECG. He calculated the standard deviation of the time series and the standard deviation of successive differences, and the length of the ellipse that characterized the Poincaré plot. He used these indicators to discriminate AFib from non-AFib . Shouldice et al. made feature vectors from inter-beat intervals, and then applied Fisher's linear discriminant method to estimate the likelihood of a block of inter-beat intervals containing the paroxysmal AFib . Kikillus et al. tried to detect AFib using a method of neural network. They calculated 25 parameters of time domain, frequency and non-linear domain, with which they applied two neural networks to decide whether the input ECG implied AFib .
If the paroxysm of AFib occurs, variability of the inter-beat intervals increases from the onset to the end of AFib ; hence, we analyzed the pulse-beat patterns to detect AFib. If the input data is contaminated with noise, it's difficult to discriminate fibrillatory wave from the noise; on the other hand, the pulse beat patterns in ECG are less likely to be influenced by baseline wandering and noise because they have clear appearances. We tried to design an algorithm which could be installed in a portable heart monitoring system since we cannot predict when the paroxysm of AFib will come about. It should endure noise well to diagnose AFib in a mobile situation. In this regard, we focused on dynamics of the inter-beat intervals to detect the onset of AFib.
We used two databases, Computers in Cardiology challenge 2001 and 2004 (CinC 2001, 2004) of physionet [17, 18]. The CinC 2001 database includes both AFib and non-AFib data files. These files were made from 24 hour ECG by cutting appropriate segments, and came from 48 different people. The files whose names begin with 'n' contain the ECG data from people who do not have any AFib. However, those people had several diseases only except AFib or else they were normal. Even numbered files whose names begin with 'p' and end with 'c' contain AFib ECG data. We dissected each data into one minute quantity to analyze easily and took only first one minute amount. The CinC 2004 database includes only AFib data files. Each file of the database had one minute amount of data. Sampling frequency of each file was 128 Hz. Each ECG data has two simultaneous components which record two different leads of ECG. We chose the first component. ECG files 'n27' and 'n27c' from the CinC 2001 database had too much noise, so we could not detect heart beat well; hence, we omitted two files. We used 25 AFib and 98 non-AFib data files from the CinC 2001 database and 80 AFib data files from the CinC 2004 database. Almost two non-AFib data files were obtained from one person and only one AFib data file was obtained from each patient. There was no explanation about the number of patients in the CinC 2004 database.
where 0 means a zero vector. Figure 1(c) is a result obtained by subtracting the waveform of Figure 1(b) from the waveform of Figure 1(a). We can see the waveform of Figure 1(c) is leveled out and the details of the waveform were preserved with respect to the waveform of Figure 1(a).
The Poincaré plot applicable to discrete data is closely related to a conventional phase plane of continuous data. If an x-coordinate of a point in Poincaré plot is x1, y-coordinate of the point is mathematically related to . If the x axis of phase plane is x, the y axis corresponds to .
Figure 3 describes the procedure of building a Poincaré plot. Figure 3(a) indicates an ECG data containing a premature ventricular contraction (PVC); in addition, it represents the inter-beat intervals I1, I2, I3, I4, I5, I6. Figure 3(b) describes the Poincaré plot made from these inter-beat intervals. This Poincaré plot has the points of (I1, I2), (I2, I3), (I3, I4), (I4, I5), (I5, I6), and we drew the lines between the consecutive points to observe the dynamics more easily. The points revolve clockwise and make a wedge-shaped diagram. This is because the inter-beat intervals changed around the PVC.
Typical patterns of Poincaré plots
Poincaré plot in case of AFib
Mean stepping increment of inter-beat intervals
Dispersion of points around diagonal line in Poincaré plot
Number of clusters in Poincaré plot
To determine the number of clusters in Poincaré plot, we developed a clustering method based on spectral graph theory . The Poincaré plots in Figure 4(b) and 5(b) show that non-AFib data sets have a limited number of clusters. On the other hand, Figure 6(b) shows that AFib data sets can have many clusters or just one conglomerate lump.
Correcting faults in QRS complex detection
We tried to correct this by identifying the triangle in the Poincaré plot as follows. Let us designate two consecutive points as (I j , Ij+1) and (Ij+1, Ij+2). We calculated the coordinates of the central point (a, a) on diagonal line in Poincaré plot in the above section. We can think of the central point (a, a) as a new origin; then the coordinates of two consecutive points will be (I j - a, Ij+1- a) and (Ij+1- a, Ij+2- a) with regard to this new origin. We identified a mistake done by the QRS complex detector when the distances , are similar and the coordinates of the middle point are positive and this middle point is located near the diagonal line and the smaller of , is substantially large.
Support vector machine classifier
The parameters C and γ were selected by an automatic tool provided by a support vector machine program . We gave an option '-v 2' to this automatic tool.
We tried to detect AFib using the above three feature measures. First we extracted inter-beat intervals from an ECG data by the wavelet method  and made Poincaré plot using these inter-beat intervals. If there were a limited number of clusters in Poincaré plot like the Figure 5(b), we concluded the Poincaré plot implied non-AFib because the small plural number of clusters meant the ECG data had some pattern.
We know that the Poincaré plot of the Figure 4(b) has the points accumulated around a central point and the Poincaré plot of the Figure 6(b) has the points scattered all over the place. The Poincaré plot of the Figure 4(b) has one cluster, and the estimated number of clusters in the Figure 6(b) is one or too many. If there was an error of calculating the inter-beat intervals from AFib data, we often got the result that the number of clusters in the Poincaré plot was too many. We had to deal with this error situation.
The leave-one-out 4-fold cross-validation.
Mean sensitivity and mean specificity were obtained by Mean sensitivity = Total True Positive/(Total True Positive + Total False Negative), Mean specificity = Total True Negative/(Total False Positive + Total True Negative). Mean sensitivity and mean specificity were 91.4% and 92.9% respectively.
We tried to propose an automated detection method which did not require any human intervention. Usually an ECG data is recorded for some period of time; then, an expert reads the data to find an abnormality in it. However, this process is tedious and inconvenient for patients since they should be in some place equipped with an electrocardiograph. Because pulse beats of ventricles are less likely to be influenced by baseline wandering and noise, we used the inter-beat intervals to diagnose AFib. It will be useful to make a real time portable monitoring electrocardiograph because we cannot predict when the paroxysm of AFib will come about. Our algorithm requires only one lead of ECG to acquire inter-beat intervals. We tried to design our algorithm not using any specific threshold.
Heart rate variability is closely related to homeostasis of the autonomous nervous system. The dynamics of inter-beat intervals come to change after the onset of AFib. Therefore we used the Poincaré plot because it was a useful tool in studying dynamics of ECG data. We found that people without any AFib showed some patterns in the Poincaré plots and these patterns were regular. The plots of AFib patients, however, were very irregular and changed too much in the course of time.
We described feature selection process and classifier using clustering and support vector machine method. We captured inter-beat intervals from an input ECG data, and drew Poincaré plot using them. By analyzing the Poincaré plot, we obtained three feature measures: the number of clusters, mean stepping increment of inter-beat intervals, and dispersion of the points in the plot. These three feature measures were dimensionless quantities.
There are some limitations in this paper. First, one feature, the dispersion of points in Poincaré plot, was not powerful to discriminate AFib from non-AFib as you can see in the Figure 9. Second, the number of data used in this paper was somewhat small. Third, the approach taken by us is susceptible to performance of QRS complex detector. If the QRS complex detector made a mistake and we failed to correct it, the accuracy of this algorithm would greatly fall off.
This work was supported by the systems biology infrastructure establishment grant provided by Gwangju Institute of Science & Technology, Korea and the Center for Distributed Sensor Network at GIST.
- Furberg C, Psaty B, Manolio T, Gardin J, Smith V, Rautaharju P: Prevalence of atrial fibrillation in elderly subjects. The American Journal of Cardiology 1994,74(3):236–241. 10.1016/0002-9149(94)90363-8View ArticleGoogle Scholar
- Bollmann A, Lombardi F: Electrocardiology of Atrial Fibrillation. Current Knowledge and Future Challenges. IEEE Engineering in Medicine and Biology Magazine 2006,25(6):15–23. 10.1109/EMB-M.2006.250504View ArticleGoogle Scholar
- Al-Rawahi N, Green M: Diagnosis of Supraventricular Tachycardia. Journal of The Association of Physicians of India 2007, (Suppl 55):21–24.Google Scholar
- Takahashi N, Seki A, Imataka K, Fujii J: Clinical Features of Paroxysmal Atrial Fibrillation. An Observation of 94 Patients [abstract]. Japanese Heart Journal 1981, 22: 143–149.View ArticleGoogle Scholar
- Konings KT, Smeets JL, Penn OC, Wellens HJ, Allessie MA: Configuration of Unipolar Atrial Electrograms During Electrically Induced Atrial Fibrillation in Humans. Circulation 1997,95(5):1231–1241.View ArticleGoogle Scholar
- Jacquemet V, Van Oosterom A, Vesin JM, Kappenberger L: Analysis of Electrocardiograms During Atrial Fibrillation. A Biophysical Model Approach. IEEE Engineering In Medicine And Biology Magazine 2006,25(6):79–88. 10.1109/EMB-M.2006.250511View ArticleGoogle Scholar
- Sudlow M, Rodgers H, Kenny RA, Thomson R: Identification of patients with atrial fibrillation in general practice: a study of screening methods. BMJ 1998, 317: 327–328.View ArticleGoogle Scholar
- Somerville S, Somerville J, Croft P, Lewis M: Atrial fibrillation: a comparison of methods to identify cases in general practice. Br J Gen Pract 2000, 50: 727–729.Google Scholar
- Mant J, Fitzmaurice DA, Hobbs FDR, Jowett S, Murray ET, Holder R, Davies M, Lip GYH: Accuracy of diagnosing atrial fibrillation on electrocardiogram by primary care practitioners and interpretative diagnostic software: analysis of data from screening for atrial fibrillation in the elderly (SAFE) trial. BMJ 2007, 335: 380. 10.1136/bmj.39227.551713.AEView ArticleGoogle Scholar
- Xu W, Tse HF, Chan F, Fung P, Lee K, Lau CP: New Bayesian Discriminator for Detection of Atrial Tachyarrhythmias. Circulation 2002, 105: 1472–1479. 10.1161/01.CIR.0000012349.14270.54View ArticleGoogle Scholar
- Petrucci E, Balian V, Filippini G, Mainardi L: Atrial Fibrillation Detection Algorithms for Very Long Term ECG Monitoring. Computers in Cardiology 2005, 2005: 623–626. full_textGoogle Scholar
- Kikillus N, Hammer G, Wieland S, Bolz A: Algorithm for Identifying Patients with Paroxysmal Atrial Fibrillation without Appearance on the ECG. Conf Proc IEEE Eng Med Biol Soc 2007, 2007: 275–278.Google Scholar
- Thuraisingham R: An electrocardiogram marker to detect paroxysmal atrial fibrillation. Journal of Electrocardiology 2007, 40: 344–347. 10.1016/j.jelectrocard.2006.10.061View ArticleGoogle Scholar
- Shouldice R, Heneghan C, de Chazal P: Automated Detection of Paroxysmal Atrial Fibrillation from Inter-Heartbeat Intervals. Conf Proc IEEE Eng Med Biol Soc 2007, 2007: 686–689.Google Scholar
- Kikillus N, Schweikert M, Bolz A: Identifying Patients Suffering From Atrial Fibrillation During Atrial Fibrillation and Non-Atrial Fibrillation Episodes. 4th European Conference of the International Federation for Medical and Biological Engineering 2008, 1349–1352.Google Scholar
- Petrutiu S, Ng J, Nijm GM, Al-Angari HM, Swiryn S, Sahakian AV: Atrial Fibrillation and Waveform Characterization. A Time Domain Perspective in the Surface ECG. IEEE Engineering in Medicine and Biology Magazine 2006,25(6):24–30. 10.1109/EMB-M.2006.250505View ArticleGoogle Scholar
- The Paroxysmal Atrial Fibrillation Prediction Challenge Database, Computers in Cardiology Challenge 2001[http://www.physionet.org/physiobank/database/afpdb]
- The Atrial Fibrillation Termination Challenge Database, Computers in Cardiology Challenge 2004[http://www.physionet.org/physiobank/database/aftdb]
- Park J: Analysis of Heart Beat Irregularity for Atrial Fibrillation by Wavelet Method of QRS Complex Detection. In Master's thesis. Korea University of Science and Technology, Medical Physics; 2008.Google Scholar
- Kasmacher H, Wiese S, Lahl M: Monitoring the Complexity of Ventricular Response in Atrial Fibrillation. Discrete Dynamics In Nature And Society 2000, 4: 63–75. 10.1155/S1026022600000066View ArticleGoogle Scholar
- Zemaityte D, Varoneckas G, Ozeraitis E, Podlipskyte A, Valyte G, Zakarevicius L: Heart Rate Poincare Plots and Their Hemodynamic Correlates: Discrimination Between Sinus and Ectopic Rhythms. Biomedicine 2001,1(2):80–89.Google Scholar
- Chung F: Spectral graph theory. American Mathematical Society; 1997.MATHGoogle Scholar
- Chih-Chung Chang and Chih-Jen Lin, LIBSVM: a library for support vector machines, 2001[http://www.csie.ntu.edu.tw/~cjlin/libsvm]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.