- Research
- Open Access

# Atrial fibrillation detection by heart rate variability in Poincare plot

- Jinho Park
^{1}, - Sangwook Lee
^{2}and - Moongu Jeon
^{1}Email author

**8**:38

https://doi.org/10.1186/1475-925X-8-38

© Park et al; licensee BioMed Central Ltd. 2009

**Received:**17 September 2009**Accepted:**11 December 2009**Published:**11 December 2009

## Abstract

### Background

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.

### Methods

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.

### Results

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.

### Conclusions

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.

## Keywords

- Support Vector Machine
- Heart Rate Variability
- Step Increment
- Radial Basis Function Kernel
- Support Vector Machine Method

## Background

There is a growing tendency that atrial fibrillation (AFib) related disease affects quality of life. Its risk increases with age [1]; in fact, AFib is one of the most common types of arrhythmia in clinical practice [2]. 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 [3]. The more serious effect of AFib is formation of blood clots by congestion of blood in atria [2]. 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 [4].

Electrical remodeling is one of the features of AFib and it is related to decreased conduction velocity of electricity signals [2]. 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 [5]. 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 [6]. 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 [7]. Somerville et al. reported a screening result of AFib; Sensitivity and specificity were 100% and 77% respectively [8]. 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 [9].

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 [10]. 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 [11]. 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 [12]. 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 [13]. 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 [14]. 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 [15].

If the paroxysm of AFib occurs, variability of the inter-beat intervals increases from the onset to the end of AFib [16]; 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.

## Methods

### ECG data

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.

### Inter-beat intervals

*A*

_{ N }is an approximation coefficient vector, and

*D*

_{ i }, (

*i*= 1,...,

*N*) is a detail coefficient vector. We chose one detail coefficient vector

*D*

_{ i }by a criterion [19], and assigned zeros to the detail coefficient vectors

*D*

_{i-1},

*D*

_{i-2},⋯,

*D*

_{1}. Figure 1(b) shows a waveform obtained by applying inverse wavelet transform to the coefficient vectors

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).

### Poincaré plots

*I*

_{1},

*I*

_{2},

*I*

_{3},

*I*

_{4},

*I*

_{5},⋯,

*I*

_{ n }like Figure 3(a), we can make a Poincaré plot that is composed of the points (

*I*

_{1},

*I*

_{2}), (

*I*

_{2},

*I*

_{3}), (

*I*

_{3},

*I*

_{4}), (

*I*

_{4},

*I*

_{5}),⋯, (

*I*

_{n-1},

*I*

_{ n }). We connected the consecutive points with lines to observe dynamics of the inter-beat intervals.

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 *x*_{1}, *y*-coordinate of the point is mathematically related to
[20]. 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 *I*_{1}, *I*_{2}, *I*_{3}, *I*_{4}, *I*_{5}, *I*_{6}. Figure 3(b) describes the Poincaré plot made from these inter-beat intervals. This Poincaré plot has the points of (*I*_{1}, *I*_{2}), (*I*_{2}, *I*_{3}), (*I*_{3}, *I*_{4}), (*I*_{4}, *I*_{5}), (*I*_{5}, *I*_{6}), 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

*O*means the QRS complex detector found the time position corresponding to the ventricular activity. Figure 5(a) shows some PVCs exist. The inter-beat intervals change around the PVCs. This is represented in the Figure 5(b) as a wedge-shaped Poincaré plot. This type of Poincaré plot is also reported in Zemaityte et al.'s paper [21]. The difference between the Poincaré plot in this paper and the plot in Zemaityte et al.'s paper is whether the lines are drawn or not between the consecutive points in the plots.

#### Poincaré plot in case of AFib

### Feature selection

#### Mean stepping increment of inter-beat intervals

*I*

_{1},

*I*

_{2},

*I*

_{3},

*I*

_{4},

*I*

_{5},⋯,

*I*

_{ n }. The points in Poincaré plot will be (

*I*

_{1},

*I*

_{2}), (

*I*

_{2},

*I*

_{3}),⋯, (

*I*

_{n-1},

*I*

_{ n }) in order. If we designate two consecutive points as (

*I*

_{ j },

*I*

_{j+1}) and (

*I*

_{j+1},

*I*

_{j+2}), the distance between two points in the Poincaré plot will be . We calculated mean value of these quantities as . This implies rate of change of the inter-beat intervals in the Poincaré plot. To normalize this and make this quantity dimensionless, we divided it by mean inter-beat interval, . We defined next quantity as the mean stepping increment of the inter-beat intervals.

#### Dispersion of points around diagonal line in Poincaré plot

*I*

_{1},

*I*

_{2},

*I*

_{3},

*I*

_{4},

*I*

_{5},⋯,

*I*

_{ n }, the points of the Poincaré plot consist of (

*I*

_{1},

*I*

_{2}), (

*I*

_{2},

*I*

_{3}), (

*I*

_{3},

*I*

_{4}),⋯, (

*I*

_{n-1},

*I*

_{ n }). We tried to find a central point (

*x*,

*x*) minimizing sum of distance squares from this point to all the other points in the Poincaré plot. If we designate this sum as

*E*(

*x*), this will be represented as follows.

*x*. From , we found the central point (

*a*,

*a*) as follows.

*I*

_{ j },

*I*

_{j+1}) to the diagonal line

*y*=

*x*is represented as . Standard deviation of these terms is represented as follows.

#### 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 [22]. 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

*I*

_{1},

*I*

_{2}), (

*I*

_{2},

*I*

_{3}), (

*I*

_{3},

*I*

_{4}), (

*I*

_{4},

*I*

_{5}) where

*I*

_{1},

*I*

_{2},

*I*

_{4}and

*I*

_{5}are almost same but less than

*I*

_{3}.

We tried to correct this by identifying the triangle in the Poincaré plot as follows. Let us designate two consecutive points as (*I*_{
j
}, *I*_{j+1}) and (*I*_{j+1}, *I*_{j+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*, *I*_{j+1}- *a*) and (*I*_{j+1}- *a*, *I*_{j+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 [23]. We gave an option '-v 2' to this automatic tool.

## Results

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 [19] 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.

Test number | True Positive | True Negative | False Positive | False Negative | Sensitivity | Specificity |
---|---|---|---|---|---|---|

1 | 26 | 23 | 2 | 1 | 0.963 | 0.920 |

2 | 25 | 23 | 2 | 1 | 0.962 | 0.920 |

3 | 22 | 22 | 2 | 4 | 0.846 | 0.917 |

4 | 23 | 23 | 1 | 3 | 0.885 | 0.958 |

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.

## Conclusions

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.

## Declarations

### Acknowledgements

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.

## Authors’ Affiliations

## References

- 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]

## Copyright

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.