Abstract
Ventricular tachycardia (VT) is a potentially fatal tachyarrhythmia, which causes a rapid heartbeat as a result of improper electrical activity of the heart. This is a potentially life-threatening arrhythmia because it can cause low blood pressure and may lead to ventricular fibrillation, asystole, and sudden cardiac death. To prevent VT, we developed an early prediction model that can predict this event one hour before its onset using an artificial neural network (ANN) generated using 14 parameters obtained from heart rate variability (HRV) and respiratory rate variability (RRV) analysis. De-identified raw data from the monitors of patients admitted to the cardiovascular intensive care unit at Asan Medical Center between September 2013 and April 2015 were collected. The dataset consisted of 52 recordings obtained one hour prior to VT events and 52 control recordings. Two-thirds of the extracted parameters were used to train the ANN, and the remaining third was used to evaluate performance of the learned ANN. The developed VT prediction model proved its performance by achieving a sensitivity of 0.88, specificity of 0.82, and AUC of 0.93.
Similar content being viewed by others
Introduction
Sudden cardiac death (SCD) causes more than 300,000 deaths annually in the United States1. Coronary artery disease, cardiomyopathy, structural heart problems, Brugada syndrome, and long QT syndrome are well known causes of SCD1,2,3,4. In addition, spontaneous ventricular tachyarrhythmia (VTA) is a main cause of SCD, contributing to about 80% of SCDs5. Ventricular tachycardia (VT) and ventricular fibrillation (VF) comprise VTA. VT is defined as a very rapid heartbeat (more than 100 times per minute), which does not allow enough time for the ventricles to fill with blood between beats. VT may terminate spontaneously after a few seconds; however, in some cases, VT can progress to more dangerous or fatal arrhythmia, VF. Accordingly, early prediction of VT will help in reducing mortality from SCD by allowing for preventive care of VTA.
Several studies have reported attempts at predicting VTAs by assessing the occurrence of syncope, left ventricular systolic dysfunction, QRS (Q, R, and S wave in electrocardiogram) duration, QT (Q and T wave) dispersion, Holter monitoring, signal averaged electrocardiograms (ECGs), heart rate variability (HRV), T wave alternans, electrophysiologic testing, B-type natriuretic peptides, and other parameters or method6,7,8,9,10. Among these studies, prediction of VTAs based on HRV analysis has recently emerged and shown potential for predicting VTA11,12,13.
Previous studies have focused on the prediction of VT using HRV analysis. In addition, most studies assessed the statistical value of each parameter calculated on or prior to the VT event and parameters of control data, which were collected from Holter recordings and implantable cardioverter defibrillators (ICDs)12,14,15. However, the results were not satisfactory in predicting fatal events like VT.
To make a better prediction model of VT, it is essential to utilize multiple parameters from various methods of HRV analysis and to generate a classifier that can deal with complex patterns composed of such parameters7. Artificial neural network (ANN) is a valuable tool for classification of a database with multiple parameters. ANN is a kind of machine learning algorithms, which can be trained using data with multiple parameters16. After training, the ANN calculates an output value according to the input parameters, and this output value can be used in pattern recognition or classification. ANN has not been widely used in medical analysis since the algorithm is not intuitive for physicians. However, utilization of ANN in medical research has recently emerged17,18,19. Our group previously reported an ANN-based prediction model for VTAs by utilizing parameters obtained from HRV analysis7. The model uses RR intervals within a 5-minute window and predicts VT events 10 seconds ahead of the events. The developed prediction model showed 73.3% sensitivity, 73.8% specificity and 75.6% accuracy, which is favorable compared to the performances of previously reported predictors13. However, predicting VT 10 seconds before it occurs is not sufficiently valuable in clinical practice. The forecast time was limited to 10 seconds in our previous study because of the length of the RR interval database. The database contains only 1024 RR intervals before VT events, which is about 6 to 12 minutes in time length. Short-term HRV analysis requires at least a 5-minute window, which greatly reduced the forecast time of VT events. To predict VT events earlier than in our previous work, gathering physiological signals including ECG from patients for a longer time period prior to VT events is essential.
In our current study, we propose a model to predict VT events 1 hour in advance by utilizing ANN and parameters from HRV and RRV analysis. Such a model may help to reduce mortality from VT events.
Results
Before training ANNs, HRV and RRV parameters of VT and control groups were compared to determine statistical differences for each parameter. Mean and standard deviation (SD) were calculated as well as corresponding p-values (Table 1). Standard deviation of NN intervals (SDNN), standard deviation 2 (SD2), and respiration period variability (RPdV) were found to be statistically significant (p <0.05). In contrast, these did not show statistical difference in previous studies7.
Initially, ANN with only 11 HRV parameters was developed. In our simulation, an ANN with 5 hidden neurons in one hidden layer showed the best performance. After training, the ANN showed 73.5% accuracy for the test set. The sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were 70.6% (12/17), 76.5% (13/17), 75.0% (12/16), and 72.2% (13/18), respectively. The results are summarized in Table 2. The receiver operating characteristic (ROC) curve of the ANN is shown in Fig. 1, with an area under curve (AUC) of 0.75.
We then developed an ANN with 3 RRV parameters. An ANN with 5 hidden neurons in one hidden layer showed the best performance. The trained ANN showed 82.4% accuracy for the test set. The sensitivity, specificity, PPV, and NPV were equal to 82.4% (Table 2). The AUC of the ROC curve was 0.83 (Fig. 1). The overall performance of the ANN using RRV parameters was found to be superior to that of the ANN using HRV parameters.
Finally, an ANN was generated using all parameters of the HRV and RRV analyses. This ANN with 13 hidden neurons in one hidden layer showed the best performance. The trained ANN showed 85.3% accuracy for the test set. The sensitivity, specificity, PPV, and NPV were 88.2% (15/17), 82.4% (14/17), 83.3% (15/18), and 87.5% (14/16), respectively (Table 2). The AUC of the ROC curve was 0.93 (Fig. 1). This model achieved the best performance among the three ANNs.
Discussion
As indicated in Table 1, the standard deviation of NN intervals (SDNN), standard deviation 2 (SD2), and respiration period variability (RPdV) were statistically significant (p < 0.05). Other parameters did not show a strong statistical significance between the VT and control dataset. However, one cannot compare the VT and control group based only on these variables. The ranges in SDNN, SD2, and PRdV values determined by the mean and SD values range of the two groups overlaps a little. Therefore, one cannot simply distinguish the two groups with only SDNN, SD2, and RPdV values. In such conditions, machine-learning based classifier is a suitable solution for developing predictors.
The results demonstrate that utilizing both ECG and respiration signals increases the performance of detecting VT one hour before its occurrence. Previously reported studies on the prediction of VT and other arrhythmias used statistical index of ECGs20,21 or only HRV and heart rate parameters12,13,14,15. Cappiello et al. showed good performance (accuracy of 98.44%) in prediction of VTA before about 356 ECG beats (about 3 to 5 minutes) but they used only 32 ECGs to develop and verifying their prediction method21. One study that focused on using HRV parameters included parametric changes in HRV, such as rise in heart rate, descent in LF, rise in LF/HF, rise in VLF, and reduction in HF observed prior to the VT event12. Other studies utilized point estimates of the correlation dimension14, slope reduction in heart rate turbulence15, and heart rate acceleration13. However, these models failed to show any significant performance. The classifier developed by Thong and Paitt showed 91% specificity, but the sensitivity was only 53 ~ 69% using 208 records from 90 subjects13. While the classifier developed by Skinner showed good performance (91% sensitivity and 85% specificity), only 14 VT subjects were included in their study14.
In contrast, the proposed ANN with HRV and RRV parameters achieved an AUC of 0.93 (88.2% sensitivity and 82.4% sensitivity) using 52 VT cases. The prediction performance of the developed predictor with only HRV showed an AUC of 0.75, which is similar to the VT prediction we previously reported. This implies that prediction of VT with only HRV parameters does not yield sufficient predictive performance. Physicians report that signs of panting or breathing difficulty frequently precede VT occurrence22. Furthermore, correlations between sleep apnea and cardiac arrhythmias have been described23. One paper reported that sleep-apnea increases the risk of atrial fibrillation24. Therefore, including RRV parameters may be important in predicting a VT event.
In addition, the model we have implemented shows great value in predicting VT one hour before its occurrence. The majority of the aforementioned studies reported predictions at or a few minutes before onset. As such, our model would allow for more time to respond to arrhythmia events. If VT leads to VF and ventricular circulation is not properly supplied to the body, organs can undergo damage, which decreases patient survival rate. Therefore, prompt medical action is necessary. One hour may be enough time to visit the hospital if there is no one around to help; as such, this model may be suitable for the general population at home as well as those in hospital.
The major limitation of this study is the small number of data (total: 104 recordings), which limits the statistical power of our analysis. Although we collected the data for about 2 years using 15 patient monitors, we saw a small numbers of actual VT cases. We are currently trying to expand the real-time vital signal collection system to other ICUs to acquire more VT cases. In addition, we will try to combine more signals, such as blood pressure or photo plethysmography. We believe that this could help in developing better predictor of VTs.
Nowadays, wearable devices, such as smart watches, smart necklaces, and smart bands that measure heart rate and other parameters are emerging. Measuring respiratory rate is a challenging issue for such wearable smart devices, especially smart watches, but this may be overcome in the near future. One major goal of such devices is personal healthcare. The developed algorithm can be easily loaded on such devices to predict and warn of fatal arrhythmias in advance.
Methods
Real-time acquisition of vital signals
Real-time physiological monitoring data used in this study was collected from 15 patient monitors (IntelliVue MP40, Phillips, Netherlands) at the cardiovascular intensive care unit (CCU) at Asan Medical Center (AMC). Every measuring signal from patient monitors was digitized and wirelessly transferred to a data server. To achieve this wireless transmission of real-time signals, a data export board (IntelliVue medical information bus/RS232 interface, Phillips, Netherlands) and a serial-to-Wi-Fi interface system (CSW-H80, Sollae System, Korea) were installed on every patient monitor in the CCU. With this configuration, the server can communicate wirelessly with each patient monitor and request alarms and vital signals that the monitor then measures, such as ECG and respiratory rate. Software for the communication with patient monitors, storing data, and data review was implemented using LabVIEW (National Instrument). The software was developed with a programming guide by Philips25. The constructed real-time collection system of vital signs is shown in Fig. 2. The institutional review board at the AMC approved the acquisition of data and all other experiments with the waiver of consent (No. 2013–0634). All experiments were performed in accordance with the approved guideline and regulations.
Database
The data acquisition server stored every vital signal from 15 patients in the CCU at AMC from September 2013 to April 2015. The total number of patients who visited the CCU within this period was 2,275. To build a VT database using the above raw data, the software obtains the event time and the monitor number when any patient monitor detects a VT. Each time the VT alarm was detected in the database, lead II ECG and respiration signals of 5-minute duration at the point of the VT event (from 5 minutes before the VT event to the VT event) and those of 1 hour before the VT event (from 65 minutes before the VT event to 6o minutes before VT event) were extracted. Data from patients with a pacemaker or ICD were excluded from the VT database because heart rhythm can be controlled by these devices. To remove false positives, ECG signals at onset of the VTs were independently reviewed by two experienced physicians (MSS, GBN). Signals collected at 1 hour before VTs were checked after removing false alarms. If there was a missing record or when a fatal error on the signals was detected, the case was also removed from the VT database. To build a control dataset, 5-minute-long ECG and respiration signals were collected when the patient was stable before and after 12 hours from the collecting point. Finally, the VT database consisted of 52 recordings prior to one hour from VT occurrence and 52 control recordings (total: 104 recordings). The data collection process is shown in Fig. 3.
Preprocessing and parameter extraction
Short-term HRV analysis was conducted on ECG signals from the database. Likewise, RRV analysis was also done by using respiratory signals from the database. To perform HRV analysis, R-peaks were detected first from the collected ECG using the detection algorithm developed by Sergey Chernenko26. Respiratory signals also underwent peak detection with the peak-detection algorithm by Eli Billauer27. Before extracting the HRV parameter, ectopic beats were removed as in a previous study7. The integrated pulse frequency modulation model28 was applied to handle the ectopic beats in the database. After revising ectopic beats, the time domain parameter and Poincaré nonlinear parameter were extracted. The frequency domain parameters were extracted by detrending the RR interval using time-varying finite impulse response high-pass filter after resampling at 7 Hz and cubic spline interpolation29. Power spectral density (PSD) of the resampled data was computed using Welch’s periodogram with a 512 points Hann window overlapped at 50%30. The spectral power was then calculated at a very low frequency (VLF, 0–0.04 Hz), low frequency (LF, 0.04–0.15 Hz), and high frequency (HF, 0.15–0.4 Hz). In extracting RRV parameters, the respiration signal was filtered with a band-pass filter of 0.1–0.5 Hz to remove noise. After that, the respiratory rate was calculated by measuring the distance between the time domain at positive peaks31. Three parameters of respiration period mean (RPdM), respiration period standard deviation (RPdSD), and respiration period variability (RPdV) were then calculated as described in a previous study32. Detailed description of other parameters used the study are summarized in Table 3.
Artificial neural network training
After extracting the 14 parameters from the database as in Table 3, two-thirds of the data (total 70 - VT: 35, control: 35) randomly selected from the complete dataset were utilized to train ANN for classification of VT events (training set). The remaining one-third of the data (total 34 - VT: 17, control: 17) was used to evaluate the performance of the trained ANN (test set).
A back propagation learning rule and a perceptron structure were used in training of ANN. The output value in response to input was set to 1 in the VT group and −1 in the control group. Training was stopped when the mean square error fell below 10−5. Finally, three ANNs (trained with parameters of HRV, RRV, and HRV + RRV, respectively) were determined, and the performance of each model was validated. Signal processing and ANN procedures were implemented directly with MATLAB 2012 (MathWorks).
Additional Information
How to cite this article: Lee, H. et al. Prediction of Ventricular Tachycardia One Hour before Occurrence Using Artificial Neural Networks. Sci. Rep. 6, 32390; doi: 10.1038/srep32390 (2016).
References
Zipes, D. P. & Wellens, H. J. J. Sudden Cardiac Death. Circulation 98, 2334–2351 (1998).
Borjesson, M. & Pelliccia, A. Incidence and aetiology of sudden cardiac death in young athletes: An international perspective. Br. J. Sports Med. 43, 644–648 (2009).
Corrado, D., Basso, C., Fau - Rizzoli, G., Rizzoli, G., Fau - Schiavon, M., Schiavon, M., Fau - Thiene, G. & Thiene, G. Does sports activity enhance the risk of sudden death in adolescents and young adults? J. Am. Coll. Cardiol. 42, 1959–1963 (2003).
Ellsworth, E. G. & Ackerman, M. J. The changing face of sudden cardiac death in the young. Heart Rhythm 2, 1283–1285 (2005).
De Luna, A. B., Coumel, P. & Leclercq, J. F. Ambulatory sudden cardiac death: mechanisms of production of fatal arrhythmia on the basis of data from 157 cases. Am. Heart J. 117, 151–159 (1989).
Hu, W. et al. Deceleration and acceleration capacities of heart rate associated with heart failure with high discriminating performance. Sci. Rep. 6, 23617, doi: 10.1038/srep23617 (2016).
Joo, S., Choi, K. J. & Huh, S. J. Prediction of spontaneous ventricular tachyarrhythmia by an artificial neural network using parameters gleaned from short-term heart rate variability. Expert Syst. Appl. 39, 3862–3866 (2012).
Lane, R. E., Cowie, M. R. & Chow, A. W. C. Prediction and prevention of sudden cardiac death in heart failure. Heart 91, 674–680 (2005).
Piccini, J. P. et al. Predictors of sudden cardiac death change with time after myocardial infarction: Results from the VALIANT trial. Eur. Heart J. 31, 211–221 (2010).
Tereshchenko, L. G., Fetics, B. J., Domitrovich, P. P., Lindsay, B. D. & Berger, R. D. Prediction of ventricular tachyarrhythmias by intracardiac repolarization variability analysis. Circ. Arrhythm. Electrophysiol. 2, 276–284 (2009).
Bilgin, S., Çolak, O. H., Polat, O. & Koklukaya, E. Estimation and evaluation of sub-bands on LF and HF base-bands in HRV for Ventricular Tachyarrhythmia patients. Expert Syst. Appl. 36, 10078–10084 (2009).
Reed, M. J., Robertson, C. E. & Addison, P. S. Heart rate variability measurements and the prediction of ventricular arrhythmias. QJM 98, 87–95 (2005).
Thong, T. & Raitt, M. H. Predicting imminent episodes of ventricular tachyarrhythmia using heart rate. PACE Pacing Clin. Electrophysiol. 30, 874–884 (2007).
Skinner, J. E., Pratt, C. M. & Vybiral, T. A reduction in the correlation dimension of heartbeat intervals precedes imminent ventricular fibrillation in human subjects. Am. Heart J. 125, 731–743 (1993).
Watanabe, M. A. Heart rate turbulence slope reduction in imminent ventricular tachyarrhythmia and its implications. J. Cardiovasc. Electrophysiol. 17, 735–740 (2006).
Solomatine, D. P. & Dulal, K. N. Model trees as an alternative to neural networks in rainfall-runoff modelling. Hydrolog. Sci. J. 48, 399–412 (2003).
Wei, J., Chen, T., Liu, G. & Yang, J. Higher-order Multivariable Polynomial Regression to Estimate Human Affective States. Sci. Rep. 6, 23384, doi: 10.1038/srep23384 (2016).
Lippmann, R. P. Pattern classification using neural networks. IEEE Commun. Mag. 27, 47–50, 59 (1989).
Patel, J. L. & Goyal, R. K. Applications of artificial neural networks in medical science. Curr. Clin. Pharmacol. 2, 217–226 (2007).
Koulaouzidis, G. et al. Prompt and accurate diagnosis of ventricular arrhythmias with a novel index based on phase space reconstruction of ECG. Int. J. Cardiol. 182 (2015).
Cappiello, G. et al. A statistical index for early diagnosis of ventricular arrhythmia from the trend analysis of ECG phase-portraits. Physiol. Meas. 36 (2014).
Schein, R. M., Hazday, N., Pena, M., Ruben, B. H. & Sprung, C. L. Clinical antecedents to in-hospital cardiopulmonary arrest. Chest 98, 1388–1392 (1990).
Mehra, R. et al. Nocturnal arrhythmias across a spectrum of obstructive and central sleep-disordered breathing in older men: Outcomes of sleep disorders in older men (MrOS sleep) study. Arch. Intern. Med. 169, 1147–1155 (2009).
Padeletti, M., Zacà, V., Mondillo, S. & Jelic, S. Sleep-disordered breathing increases the risk of arrhythmias. J. Cardiovasc. Med. 15, 411–416 (2014).
Philips. IntelliVue Patient Monitor X2, MP Series, MX Series: Data export interface programming guide (2008).
Chernenko, S. ECG processing: R-peaks detection, http://www.librow.com/cases/case-2 (2007).
Billauer, E. peakdet: Peak detection using MATLAB, http://billauer.co.il/peakdet.html (2012).
Solem, K., Laguna, P., Fau - Sornmo, L. & Sornmo, L. An efficient method for handling ectopic beats using the heart timing signal. IEEE Trans. Biomed. Eng. 53, 13–20 (2006).
Tarvainen, M. P., Ranta-Aho Po Fau - Karjalainen, P. A. & Karjalainen, P. A. An advanced detrending method with application to HRV analysis. IEEE Trans. Biomed. Eng. 49, 172–175 (2002).
Welch, P. D. The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging Over Short, Modified Periodograms. IEEE Trans. Audio Electroacoust. 15, 70–73 (1967).
Malik, M. et al. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Eur. Heart J. 17 (1996).
Wittenboer, G. K. W. & Dixhoorn, J. Respiratory variability and psychological well-being in schoolchildren. Behav. Modif. 27, 653–670 (2003).
Acknowledgements
This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2014R1A1A2053862) and by a grant (2013–7000) from the Asan Institute for Life Sciences, Seoul, Korea.
Author information
Authors and Affiliations
Contributions
H.L. built up classification models and tested the models. S.J conducted the model building and testing. S.-Y.S. collected data from patient monitors. M.S. reviewed medical records of enrolled patients and checked ECGs with G.-B.N. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Lee, H., Shin, SY., Seo, M. et al. Prediction of Ventricular Tachycardia One Hour before Occurrence Using Artificial Neural Networks. Sci Rep 6, 32390 (2016). https://doi.org/10.1038/srep32390
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep32390
This article is cited by
-
An Optimal Methodology for Early Prediction of Sudden Cardiac Death Using Advanced Heart Rate Variability Features of ECG Signal
Arabian Journal for Science and Engineering (2024)
-
Artificial Intelligence and Machine Learning in Electrophysiology—a Short Review
Current Treatment Options in Cardiovascular Medicine (2023)
-
Applications of Artificial Intelligence in Cardiac Electrophysiology and Clinical Diagnosis with Magnetic Resonance Imaging and Computational Modeling Techniques
Journal of Shanghai Jiaotong University (Science) (2023)
-
Machine learning techniques for arrhythmic risk stratification: a review of the literature
International Journal of Arrhythmia (2022)
-
Big Data in electrophysiology
Herzschrittmachertherapie + Elektrophysiologie (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.