© 2005 The European Society of Cardiology. Published by Elsevier Ltd. All rights reserved.
Prediction of maintenance of sinus rhythm after electrical cardioversion of atrial fibrillation by non-deterministic modelling
ohara
i
baCardiology Department, Hospital Celje Slovenia; bLaboratory for Intelligent Manufacturing Systems, Faculty of Mechanical Engineering Maribor, Slovenia; cDepartment for Intensive Internal Medicine, General Hospital Celje Oblakova 5 3000 Celje, Slovenia
Manuscript submitted 30 March 2004. Revision received 14 July 2005. Accepted after revision 18 April 2005.
*Corresponding author. Tel.: +386 34233000; fax: +386 35481204.E-mail address: matej.podbregar{at}guest.arnes.si (M. Podbregar)
| Abstract |
|---|
|
|
|---|
AIMS: Atrial fibrillation (AF) is the most common rhythm disorder. Because of the high recurrence rate of AF after cardioversion and because of potential side effects of electrical cardioversion, it is clinically important to predict persistence of sinus rhythm after electrical cardioversion before it is attempted. The aim of our study was the development of a mathematical model by "genetic" programming (GP), a non-deterministic modelling technique, which would predict maintenance of sinus rhythm after electrical cardioversion of persistent AF.
PATIENTS AND METHODS: Ninety-seven patients with persistent AF lasting more than 48 h, undergoing the first attempt at transthoracic cardioversion were included in this prospective study. Persistence of AF before the cardioversion attempt, amiodarone treatment, left atrial dimension, mean, standard deviation and approximate entropy of ECG RR intervals were collected. The data of 53 patients were randomly selected from the database and used for GP modelling; the other 44 data sets were used for model testing.
RESULTS: In 23 patients sinus rhythm persisted at 3 months. In the other 21 patients sinus rhythm was not achieved or its duration was less than 3 months. The model developed by GP failed to predict maintenance of sinus rhythm at 3 months in one patient and in six patients falsely predicted maintenance of sinus rhythm. Positive and negative likelihood ratios of the model for testing data were 4.32 and 0.05, respectively. Using this model 15 of 21 (71.4%) cardioversions not resulting in sinus rhythm at 3 months would have been avoided, whereas 22 of 23 (95.6%) cardioversions resulting in sinus rhythm at 3 months would have been administered.
CONCLUSION: This model developed by GP, including clinical data, ECG data from the time-domain and nonlinear dynamics can predict maintenance of sinus rhythm. Further research is needed to explore its utility in the present or an expanded form.
Key Words: atrial fibrillation, electrical cardioversion, prediction
| Introduction |
|---|
|
|
|---|
Atrial fibrillation is the most commonly sustained cardiac arrhythmia in clinical practice; it is characterized by rapid ineffective atrial activity with irregular ventricular contractions [1]
AF may cause systemic thromboembolic complications, decrease exercise capacity, impair ventricular function, reduce quality of life and incur significant health care costs [2,
3]
. After adjusting for the underlying cardiac condition, AF is associated with a 1.51.9-fold increase in risk of mortality in both men and women across a wide spectrum of ages [4]
.
For many patients, maintenance of sinus rhythm is the main therapeutic goal. In patients with persistent AF repeated electrical cardioversions and prophylactic antiarrhythmic drugs are used to maintain sinus rhythm [5]
. The success rate of electrical cardioversions ranges from 64% to 96% [6]
. Although the success rate of electrical cardioversion is high, recurrence of AF is common, especially during the first 2 weeks following the procedure [7]
. The percentage of recurrences of AF within 2 min after electrical cardioversion is 10% [8]
. Electrical cardioversion has the potential of causing severe side effects, such as postshock bradycardia, malignant ventricular arrhythmias [9]
, arterial thromboembolism [10]
and complications related to anaesthesia [11]
. It has also minor complications such as local skin irritation or reversible muscle pain [12]
. Because of the high recurrence rate of AF and because of potential side effects of electrical cardioversion, it would be clinically useful to be able to predict persistence of sinus rhythm after electrical cardioversion before it is attempted.
There are essentially two kinds of predictive models: deterministic (mathematical models, empirical models, and computer simulation models) and non-deterministic (models developed by genetic methods, neural network models, and models based on chaos theory and soft logic theory), and each has its advantages and disadvantages [13]
. In general, when deterministic modelling is used the models obtained are the result of strict mathematical rules or they are set in advance. In that case, the goal is merely to discover a set of numerical coefficients for a model whose form has been pre-specified. However, nowadays more and more processes and systems are modelled and optimized using non-deterministic approaches. This is due to the high degree of complexity of the systems, and consequently, the inability to study them successfully with conventional methods only. In non-deterministic modelling of systems, there are no precise, strict mathematical rules. For example, in "genetic" programming (GP), no assumptions about the form, size, and complexity of models are made in advance. They are left to stochastic, self-organized, intelligent, and non-centralized evolutionary processes [14]
.
The aim of our study was to develop a mathematical model by GP, which would predict maintenance of sinus rhythm after electrical cardioversion of persistent AF.
| Methods |
|---|
|
|
|---|
Patients
Ninety seven patients with persistent AF lasting more than 48 h, during a 1-year period (January 2001 to January 2002), undergoing the first attempt at transthoracic electrical cardioversion were included in this prospective study. They were all treated at the cardiology department in a 860-bed Community General and Teaching Hospital. The study was approved by the Institutional Review Board. The patients gave informed consent.
Transthoracic echocardiography was performed in all patients before cardioversion.
Wet polymer gel pads (M3502A, Agilent Technologies, Andover, MA, USA) were applied to the standard antero-lateral position. The anterior pad was placed at the upper right sternal border and the lateral pad over the cardiac apex.
Patients received damped sine wave monophasic sequential shocks of 100, 200, 300 and 360 J (Code Master XL+, Hewlett Packard, USA). All shocks were delivered during exhalation. During the procedure the patients were sedated with midazolam. For the purpose of the study cardioversion was considered successful if atrial-P waves were unmistakably identified
2 min after the shock.
Appropriate anticoagulation for at least 3 weeks or transoesophageal echocardiography to exclude the presence of atrial thrombi before the procedure and anticoagulation with a therapeutic International Normalized Ratio (INR) of 2.0 to 3.0 for at least 4 weeks after the procedure was mandatory. All patients were required to have 12-lead ECG within 1 and 3 months of the cardioversion.
Patients were considered ineligible for the study if they were pregnant, were <18 years of age, were having AF of less than 48 h, had suffered from acute coronary syndrome within 3 weeks prior to cardioversion or if pharmacological cardioversion had been tried prior to electrical cardioversion.
Signal recording and analysis
A monophasic defibrillator was used for cardioversion and simultaneous on-line ECG recording. ECG signals of atrial fibrillation in length of 25 s before cardioversion were analyzed using MatLab R12 (MathWorks, USA).
Measurement of RR (beat to beat) intervals was made from the original time-domain ECG signal by calculating the difference between different consecutive R spikes using the same computer programme.
Mean and standard deviation of consecutive RR intervals were calculated.
Approximate entropy (ApEn), which is used in nonlinear analysis of RR, was calculated from the same RR intervals. ApEn measures the logarithmic likelihood that runs of patterns that are close to each other will remain close in the subsequent incremental comparison. A series containing many repetitive patterns has a relatively small ApEn; conversely, more random data produce higher values. Details of this method have been previously described [15]
.
"Genetic" programming
To avoid unnecessarily high classification performance of our model, data were divided into training and testing data sets. Data of 53 patients were randomized from the database and used for GP modelling. The theory of GP can be found in many books and articles that deal with evolutionary computation [16,
17]
(please see Appendix for more details).
In our study, computer programmes are in fact mathematical expressions (i.e. prediction models). The terminal genes (i.e. independent input variables) were persistence of AF before cardioversion attempt (1 = less than 3 month, 2 = between 3 and 6 months, 3 = between 6 and 12 months, 4 = more than 1 year), amiodarone treatment (yes = 1, no = 0), left atrial dimension (LA) (mm), mean (s) and SD (s) of RR intervals and ApEn. The dependant output variable was persistence of sinus rhythm at 3 months after cardioversion (1 = sinus rhythm persisted at 3 months, 0 = atrial fibrillation). The function genes were basic arithmetical operations (functions of addition, subtraction, division, and multiplication). With selected terminal and function genes the evolutionary attempts were made to construct, over several generations, as accurate a model as possible for prediction of sinus rhythm persistence at 3 months after cardioversion. The best model developed by GP from the training data was retested with the residual testing data of the 44 patients not included in the training set.
Statistics
Results were expressed as mean ± SD unless otherwise specified. Comparisons of dichotomous and continuous variables between different groups were calculated using
2 or nonparametric KruskalWallis test and Student's t test, respectively. SPSS 10.0 for Windows (SPSS Inc., USA) was used. A 2-tailed value of P < 0.05 was considered statistically significant. Classification of the performance of the decision model was determined by Bayesian analysis (true positive (TP), true negative (TN), false positive (FP), false negative (FN)) applying sensitivity (SENS) = TP/(TP + FN), specificity (SPEC) = TN/(FP + TN), positive predictive value = TP/(TP + FP), negative predictive value = TN/(FN + TN), positive likelihood ratio = SENS/(1SPEC), and negative likelihood ratio = (1SENS)/SPEC [18]
.
| Results |
|---|
|
|
|---|
Ninety-seven electrical cardioversion procedures were attempted. Sinus rhythm persisted in 89.9% of patients at 24 h after cardioversion, in 54.1% at 1 month and in 49.0% of patients at 3 months after cardioversion. The characteristics of patients are shown in Table 1. At 3 months sinus rhythm persisted more often in patients receiving amiodarone (22/48 vs. 7/49, P = 0.022). Maintenance of sinus rhythm was influenced by AF duration. At 3 months AF was registered more often compared with sinus rhythm in patients with AF duration of more than 12 months (24/49 vs. 9/48, P = 0.044).
|
The predictive model developed by GP has not excluded any of six variables included at the start of programming. The equation of the predictive model is presented in Fig. 1. The following rules were used: the result of equation more than zero predicted persistence of sinus rhythm at 3 months after the cardioversion attempt, otherwise the rhythm was predicted not to be sinus.
|
The testing data included 44 data sets before cardioversion. In 23 patients sinus rhythm persisted at 3 months. In the other 21 patients sinus rhythm was not achieved or it persisted for less than 3 months.
The model developed by GP failed to predict persistence of sinus rhythm at 3 months in one patient (false negative result) and in six patients falsely predicted (false positive result) persistence of sinus rhythm. The classification of performance of the model developed by GP on testing data is shown in Table 2.
|
| Discussion |
|---|
|
|
|---|
In this study the prediction power of a non-deterministic model developed by GP to predict persistence of sinus rhythm at 3 months after cardioversion of persistent atrial fibrillation was studied. High positive and low negative likelihood ratios of the GP model indicate that the model may have a very substantial impact on clinical decision-making through meaningful revision of performance of cardioversion [19]
In all patients, six variables from the time-domain (persistence of atrial fibrillation, mean RR interval before cardioversion, SD of RR interval before cardioversion), echocardiographic data (left atrial dimension), treatment (amiodarone) and nonlinear dynamics properties of RR intervals (ApEn) were collected. On the basis of these data alone, amiodarone and shorter duration of AF before cardioversion attempt treatment were associated with a higher rate of sinus rhythm maintenance as was previously shown [20]
.
All independent input variables included in our model developed by GP had already been studied as predictors of persistence of sinus rhythm after cardioversion or as predictors of the onset of paroxysmal atrial fibrillation [6,
21]
. Electrical cardioversion is the most effective method for termination of persistent atrial fibrillation. Rates of transthoracic cardioversion of chronic atrial fibrillation vary from 70% to more than 90%. Success rate is influenced by duration of atrial fibrillation, transthoracic impedance, initial energy selection and type of shock waveform [22
26]
. We have previously reported that the type of shock waveform does not influence the success rate of cardioversion in atrial fibrillation which has persisted for more than a year [27]
. It was also previously shown that amiodarone reduces procedures and cost related to atrial fibrillation and it reduces recurrences of atrial fibrillation after an electrical cardioversion [28
30]
. Left atrial function and size are also independent predictors of atrial fibrillation recurrence [31,
32]
. The autonomic nervous system may play an important role as a trigger for spontaneous onset of paroxysmal atrial fibrillation [33]
. Analysis of heart rate variability has become an important, non-invasive tool for assessing autonomic influence on the heart [34]
. It was recently shown that a decrease of the complexity of RR intervals and altered fractal properties (ApnE) in short-term RR interval dynamics precedes the spontaneous onset of atrial fibrillation in patients with no structural heart disease.
In our study the data were separated into training and testing data, so the "best-fit" threshold values were not used in prediction to obviate artificially good sensitivity and specificity.
The present study has at least three major limitations. First, the number of cardioversions was low. Second, a 25-s period of ECG signal just prior to cardioversion was analyzed rather than classical long term ECG signal analysis [35]
. The third weakness of our model is that it does not include data which potentially influence prediction of sinus rhythm maintenance after electrical conversion such as presence or absence of hypertension, level of C-reactive protein, presence or absence of heart failure etc. We have focused our data mainly on the ECG to develop computer programmes which could be simply integrated into defibrillators.
We have demonstrated that a non-deterministic model called genetic programming is able to develop a model, using clinical data, ECG data from the time-domain and nonlinear dynamics, which could predict maintenance of sinus rhythm after cardioversion of chronic atrial fibrillation with high positive and low negative likelihood ratios. Further research is needed to determine the utility of this model or an expanded version.
| Appendix |
|---|
|
|
|---|
Genetic programming
GP is probably the most general approach of evolutionary computational methods. In GP the structures subject to adaptation are the hierarchically organized computer programmes (organisms) the size and form of which dynamically change during simulated evolution. Computer programmes can be mathematical expressions, control strategies, plans, decision trees, state-transition rules, etc. Possible solutions in GP are all the possible computer programmes that can be composed in a recursive manner from a set of function genes and a set of terminal genes describing the problem area to be studied. The set of function genes can include e.g. arithmetical functions, Boolean functions, relation functions, program flow control functions, where the set of terminal genes, e.g., numerical constants, logical constants, variables. The aim of GP is to find the computer programme that best solves the problem.
The initial population is obtained with the creation of random computer programmes consisting of the available function and terminal genes. The next step is the evaluation of population (i.e. calculation of fitness for each organism). On the basis of fit, alternation of computer programmes with genetic operations follows. In GP the computer programs alter, in particular, with reproduction and crossover. The reproduction operation gives a higher probability of selection to more successful organisms. They are copied unchanged into the next generation. The crossover operation ensures the exchange of genetic material between computer programmes. Fig. 2 shows the operation of crossover of two computer programmes, which are in fact mathematical expressions. Two randomly selected parts of two parental organisms (in boldface) are interchanged. Thus two offspring are created.
|
After finishing the first cycle which includes: (1 ) the creation of initial population, (2) calculation of fit for each i (1vidual of the population, and (3) genetic alternation of contents of the computer programmes, an iterative repetition of points 2 and 3 follows. The evolution is terminated when a termination criterion is fulfilled. This can be a prescribed maximum number of generations or a sufficient quality of the solution.
| Acknowledgements |
|---|
|
|
|---|
We thank the patients who participated in the study and Drs Mojca Bervar, Matej Marin
ek, Dragan Kova
i
, Mojca Perkolj-Bi
ani
, Nanika
krabl-Mo
nik, Albina Dokler-Glavnik for treating the patients. | References |
|---|
|
|
|---|
[1] Benjamin EJ, Levy D, Vaziri SM, D'Agostino RB, Belanger AJ, Wolf PA. Independent risk factors for atrial fibrillation in a population-based cohort. The Framingham heart study. JAMA 1994; 271: 840844.
[2] Gronefeld G and Hohnloser SH. Towards a consensus in rate versus rhythm control for management of atrial fibrillation: insights from the PIAF trial. Card Electrophysiol Rev 2003; 7: 113117.[CrossRef][Medline]
[3] Gronefeld GC, Lilienthal J, Kuck KH, Hohnloser SH. Pharmacological intervention in atrial fibrillation (PIAF) study investigators. Impact of rate versus rhythm control on quality of life in patients with persistent atrial fibrillation. Results from a prospective randomized study. Eur Heart J 2003; 24: 14301436.
[4] Benjamin EJ, Wolf PA, D'Agostino RB, Silbershatz H, Kannel WB, Levy D. Impact of atrial fibrillation on the risk of death: the Framingham heart study. Circulation 1998; 98: 946952.
[5] Fuster V, Rydén LE, Asinger RW, et al. American College of Cardiology/American Heart Association Task Force on Practice Guidelines; European Society of Cardiology Committee for Practice Guidelines and Policy Conferences (Committee to Develop Guidelines for the Management of Patients With Atrial Fibrillation): North American Society of Pacing and Electrophysiology. ACC/AHA/ESC Guidelines for the Management of Patients With Atrial Fibrillation: Executive Summary A Report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines and the European Society of Cardiology Committee for Practice Guidelines and Policy Conferences (Committee to Develop Guidelines for the Management of Patients With Atrial Fibrillation) Developed in Collaboration With the North American Society of Pacing and Electrophysiology. Circulation 2001; 104: 21182150.
[6] Ricard P, Lévy S, Trigano J, et al. Prospective assessment of the minimum energy needed for external electrical cardioversion of atrial fibrillation. Am J Cardiol 1997; 79: 815816.[CrossRef][Web of Science][Medline]
[7] Tieleman RG, Van Gelder IC, Crijns HJ, et al. Early recurrences of atrial fibrillation after electrical cardioversion: a result of fibrillation-induced electrical re-modelling of the atria? J Am Coll Cardiol 1998; 31: 167173.
[8] Van Gelder IC, Tuinenburg AE, Schoonderwoerd BS, Tieleman RG, Crijns HJ. Pharmacologic versus direct-current electrical cardioversion of atrial flutter and fibrillation. Am J Cardiol 1999; 84: 147R151R.[Web of Science][Medline]
[9] Ebrahimi R and Rubin SA. Electrical cardioversion resulting in death from synchronization failure. Am J Cardiol 1994; 74: 100102.[CrossRef][Web of Science][Medline]
[10] Bjerkelund CJ and Orning OM. The efficacy of anticoagulant therapy in preventing embolism related to D.C. electrical conversion of atrial fibrillation. Am J Cardiol 1969; 25: 208216.
[11] Gale DW, Grissom TE, Mirenda JV. Titration of intravenous anesthetics for cardioversion: a comparison of propofol, methohexital, and midazolam. Crit Care Med 1993; 21: 15091513.[Web of Science][Medline]
[12] Groenefeld G, Ehrlich JR, Hohnloser SH. Comparison of outpatient vs inpatient direct current cardioversion of atrial fibrillation: safety, efficacy and cost savings. Eur Heart J 2003; 5:Suppl. H 1924.
[13] Dobnik-Dubovski P and Bezocnik M. Using genetic programming to predict the macroporosity of woven cotton fabrics. Text Res J 2002; 73: 187194.
[14] Bezocnik M and Balic J. A genetic-based approach to simulation of self-organizing assembly. Robot Comput Integrat Manufact 2001; 17: 113120.[CrossRef]
[15] Pincus SM. Approximate entropy as a measure of system complexity. Proc Natl Acad Sci U S A 1991; 88: 22972301.
[16] Koza JR. Genetic programming II 1994; Massachusetts The MIT Press.
[17] Banzhaf W, Nordin P, Keller ER, Francone FD. Genetic programming: an introduction 1998; San Francisco, CA Morgan Kaufmann.
[18] Gallagher EJ. Clinical utility of likelihood ratios. Ann Emerg Med 1998; 31: 391397.[CrossRef][Web of Science][Medline]
[19] Laupacis A, Sekar N, Stiell IG. Clinical prediction rules. A review and suggested modifications of methodological studies. JAMA 1997; 277: 488494.
[20] Roy D, Talajic M, Dorian P, et al. Amiodarone to prevent recurrence of atrial fibrillation. Canadian trial of atrial fibrillation investigators. N Engl J Med 2000; 342: 913920.
[21] Lévy S, Lauribe P, Dolla E., et al. A randomized comparision of external and internal cardioversion of chronic atrial fibrillation. Circulation 1992; 86: 14151420.
[22] Mittal S, Ayati S, Stein KM, et al. Transthoracic cardioversion of atrial fibrillation: comparison of rectilinear biphasic versus damped sine wave monophasic shocks. Circulation 2000; 101: 12821287.
[23] Kerber RE, Martins JB, Kienzle MG, et al. Energy, current and success in defibrillation and cardioversion: clinical studies using an automated impedance-based method of energy adjustment. Circulation 1988; 77: 10381046.
[24] Bjerregaard P, El-Shafei A, Janosik DL, et al. Double external direct-current shocks for refractory atrial fibrillation. Am J Cardiol 1999; 83: 972974.[CrossRef][Web of Science][Medline]
[25] Joglar JA, Hamdan MH, Ramaswamy K, et al. Initial energy for elective external cardioversion of persistent atrial fibrillation. Am J Cardiol 2000; 86: 348350.[CrossRef][Web of Science][Medline]
[26] Gallagher MM, Guo XH, Poloniecki JD, et al. Initial energy setting, outcome and efficiency in direct current cardioversion of atrial fibrillation and flutter. J Am Coll Cardiol 2001; 38: 14981504.
[27] Marinsek M, Larkin GL, Zohar P, et al. Efficacy and impact of monophasic versus biphasic countershocks for transthoracic cardioversion of persistent atrial fibrillation. Am J Cardiol 2003; 92: 988991.[CrossRef][Web of Science][Medline]
[28] Lumer GB, Roy D, Talajic M, et al. Amiodarone reduces procedures and costs related to atrial fibrillation in a controlled clinical trial. Eur Heart J 2002; 23: 10501056.
[29] De Simone A, De Pasquale M, De Matteis C, et al. Verapamil plus antiarrhythmic drugs reduce atrial fibrillation recurrences after an electrical cardioversion (VEPARAF study). Eur Heart J 2003; 24: 14251429.
[30] Galperin J, Elizari MV, Chiale PA, et al. Grupo de Estudio de Fibrilacion Auricular Con Amiodarona (GEFACA) Investigators. Pharmacologic reversion of persistent atrial fibrillation with amiodarone predicts long-term sinus rhythm maintenance. J Cardiovasc Pharmacol Ther 2003; 8: 179186.
[31] Mattioli AV, Bonatti S, Bonetti L, Mattioli G. Left atrial size after cardioversion for atrial fibrillation: effect of external direct current shock. J Am Soc Echocardiogr 2003; 16: 271276.[CrossRef][Web of Science][Medline]
[32] Ortiz De Murua JA, del Carmen Avila M, Ochoa C, de La Fuente L, Moreno De Vega JC, del Campo F, et al. Independent predictive factors of acute and first year success after electrical cardioversion in patients with chronic atrial fibrillation. Rev Esp Cardiol 2001; 54: 958964.[Web of Science][Medline]
[33] Coumel P. Paroxysmal atrial fibrillation: a disorder of autonomic tone? Eur Heart J 1994; 15: 916.
[34] Saul JP, Arai Y, Berger RD, Lilly LS, Colucci WS, Cohen RJ. Assessment of autonomic regulation in chronic congestive heart failure by heart rate spectral analysis. Am J Cardiol 1988; 61: 12921299.[CrossRef][Web of Science][Medline]
[35] Vikman S, Makikallio TH, Yli-Mayry S, Nurmi M, Airaksinen KE, Huikuri HV. Heart rate variability and recurrence of atrial fibrillation after electrical cardioversion. Ann Med 2003; 35: 3642.[CrossRef][Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

