An open source benchmarked toolbox for cardiovascular waveform and interval analysis

OBJECTIVE This work aims to validate a set of data processing methods for variability metrics, which hold promise as potential indicators for autonomic function, prediction of adverse cardiovascular outcomes, psychophysiological status, and general wellness. Although the investigation of heart rate variability (HRV) has been prevalent for several decades, the methods used for preprocessing, windowing, and choosing appropriate parameters lacks consensus among academic and clinical investigators. Moreover, many of the important steps are omitted from publications, preventing reproducibility. APPROACH To address this, we have compiled a comprehensive and open-source modular toolbox for calculating HRV metrics and other related variability indices, on both raw cardiovascular time series and RR intervals. The software, known as the PhysioNet Cardiovascular Signal Toolbox, is implemented in the MATLAB programming language, with standard (open) input and output formats, and requires no external libraries. The functioning of our software is compared with other widely used and referenced HRV toolboxes to identify important differences. MAIN RESULTS Our findings demonstrate how modest differences in the approach to HRV analysis can lead to divergent results, a factor that might have contributed to the lack of repeatability of studies and clinical applicability of HRV metrics. SIGNIFICANCE Existing HRV toolboxes do not include standardized preprocessing, signal quality indices (for noisy segment removal), and abnormal rhythm detection and are therefore likely to lead to significant errors in the presence of moderate to high noise or arrhythmias. We therefore describe the inclusion of validated tools to address these issues. We also make recommendations for default values and testing/reporting.

[1]  J. Richman,et al.  Physiological time-series analysis using approximate entropy and sample entropy. , 2000, American journal of physiology. Heart and circulatory physiology.

[2]  Joachim Behar,et al.  A Comparison of Single Channel Fetal ECG Extraction Methods , 2014, Annals of Biomedical Engineering.

[3]  P. Sleight,et al.  Effects of controlled breathing, mental activity and mental stress with or without verbalization on heart rate variability. , 2000, Journal of the American College of Cardiology.

[4]  Patrick R. Norris,et al.  Reduced heart rate multiscale entropy predicts death in critical illness: a study of physiologic complexity in 285 trauma patients. , 2008, Journal of critical care.

[5]  Patrick E. McSharry,et al.  Advanced Methods And Tools for ECG Data Analysis , 2006 .

[6]  Gari D Clifford,et al.  Multimodal heart beat detection using signal quality indices , 2015, Physiological measurement.

[7]  J. Brouwer,et al.  The influence of recording length on time and frequency domain analysis of heart rate variability , 1998, Computers in Cardiology 1998. Vol. 25 (Cat. No.98CH36292).

[8]  Marek Malik,et al.  Phase-rectified signal averaging for the detection of quasi-periodicities and the prediction of cardiovascular risk. , 2007, Chaos.

[9]  Gulmira Kudaiberdieva,et al.  Bispectral analysis of heart rate variability signal , 2002, 2002 11th European Signal Processing Conference.

[10]  N. Lomb Least-squares frequency analysis of unequally spaced data , 1976 .

[11]  G D Clifford,et al.  Phase-rectified signal averaging as a sensitive index of autonomic changes with aging. , 2010, Journal of applied physiology.

[12]  Jing Yan,et al.  Do the deceleration/acceleration capacities of heart rate reflect cardiac sympathetic or vagal activity? A model study , 2016, Medical & Biological Engineering & Computing.

[13]  G. Breithardt,et al.  Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. , 1996 .

[14]  A. Camm,et al.  Heart-rate turbulence after ventricular premature beats as a predictor of mortality after acute myocardial infarction , 1999, The Lancet.

[15]  Adrian R. Hill,et al.  Inter-observer and intra-observer variability of the Oxford clinical cataract classification and grading system , 2004, International Ophthalmology.

[16]  Joachim Behar,et al.  Crowd-Sourced Annotation of ECG Signals Using Contextual Information , 2013, Annals of Biomedical Engineering.

[17]  H. Mølgaard,et al.  Evaluation of the Reynolds Pathfinder II system for 24 h heart rate variability analysis. , 1991, European heart journal.

[18]  Eduardo Villacorta,et al.  Variabilidad inter e intraobservador en la valoración ecocardiográfica de la función del ventrículo derecho , 2010 .

[19]  Wojciech Zareba,et al.  Risk Stratification Using Heart Rate Turbulence and Ventricular Arrhythmia in MADIT II: Usefulness and Limitations of a 10‐Minute Holter Recording , 2004, Annals of noninvasive electrocardiology : the official journal of the International Society for Holter and Noninvasive Electrocardiology, Inc.

[20]  J. Scargle Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data , 1982 .

[21]  Shamim Nemati,et al.  Benchmarking heart rate variability toolboxes. , 2017, Journal of electrocardiology.

[22]  G. Schmidt,et al.  Heart rate turbulence. , 2003, Journal of electrocardiology.

[23]  Patrick E. McSharry,et al.  Method for generating an artificial RR tachogram of a typical healthy human over 24-hours , 2002, Computers in Cardiology.

[24]  Qiao Li,et al.  Signal processing and feature selection preprocessing for classification in noisy healthcare data , 2016 .

[25]  George B. Moody,et al.  A robust open-source algorithm to detect onset and duration of QRS complexes , 2003, Computers in Cardiology, 2003.

[26]  Lionel Tarassenko,et al.  Quantifying errors in spectral estimates of HRV due to beat replacement and resampling , 2005, IEEE Transactions on Biomedical Engineering.

[27]  A. Kadish,et al.  Uncertainty principle of signal-averaged electrocardiography. , 2000, Circulation.

[28]  James C Lacefield,et al.  Three-dimensional visualization and thickness estimation of aortic valve cusps using high-frequency ultrasound. , 2004, Physiological measurement.

[29]  M. Saeed,et al.  Estimating cardiac output from arterial blood pressurewaveforms: a critical evaluation using the MIMIC II database , 2005, Computers in Cardiology, 2005.

[30]  Georg Schmidt,et al.  Phase-rectified signal averaging method to predict intermediate perinatal outcomes in infants with very preterm fetal growth restriction – a secondary analysis of TRUFFLE-trial , 2015 .

[31]  Keith R. Godfrey,et al.  A new application for heart rate variability: diagnosing the sleep apnoea syndrome , 1998, Computers in Cardiology 1998. Vol. 25 (Cat. No.98CH36292).

[32]  Claude Gronfier,et al.  Dynamic heart rate variability: a tool for exploring sympathovagal balance continuously during sleep in men. , 1998, American journal of physiology. Heart and circulatory physiology.

[33]  R G Mark,et al.  Robust heart rate estimation from multiple asynchronous noisy sources using signal quality indices and a Kalman filter , 2008, Physiological measurement.

[34]  Patrick E. McSharry,et al.  A dynamical model for generating synthetic electrocardiogram signals , 2003, IEEE Transactions on Biomedical Engineering.

[35]  Fernando Andreotti,et al.  R-peak estimation using multimodal lead switching , 2014, Computing in Cardiology 2014.

[36]  James Xin Sun,et al.  Cardiac output estimation using arterial blood pressure waveforms , 2006 .

[37]  Federico Lombardi,et al.  Heart rate turbulence: standards of measurement, physiological interpretation, and clinical use: International Society for Holter and Noninvasive Electrophysiology Consensus. , 2008, Journal of the American College of Cardiology.

[38]  E. Villacorta,et al.  Inter- and intra-observer variability in the echocardiographic evaluation of right ventricular function. , 2010, Revista espanola de cardiologia.

[39]  H. Stanley,et al.  Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. , 1995, Chaos.

[40]  G. Clifford,et al.  Signal quality and data fusion for false alarm reduction in the intensive care unit. , 2012, Journal of electrocardiology.

[41]  Junichiro Hayano,et al.  Multi-scale heart rate dynamics detected by phase-rectified signal averaging predicts mortality after acute myocardial infarction. , 2013, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[42]  Jeffrey M. Hausdorff,et al.  Physionet: Components of a New Research Resource for Complex Physiologic Signals". Circu-lation Vol , 2000 .

[43]  L. Tarassenko,et al.  Characterizing artefact in the normal human 24-hour RR time series to aid identification and artificial replication of circadian variations in human beat to beat heart rate using a simple threshold , 2002, Computers in Cardiology.

[44]  P. McSharry,et al.  Quantifying self-similarity in cardiac inter-beat interval time series , 2005, Computers in Cardiology, 2005.

[45]  Madalena Costa,et al.  Multiscale entropy analysis of complex physiologic time series. , 2002, Physical review letters.

[46]  G. B. Young,et al.  Sleep in critically ill patients requiring mechanical ventilation. , 2000, Chest.

[47]  Qiao Li,et al.  Ventricular Fibrillation and Tachycardia Classification Using a Machine Learning Approach , 2014, IEEE Transactions on Biomedical Engineering.

[48]  Madalena Costa,et al.  Multiscale entropy analysis of biological signals. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[49]  Mika P. Tarvainen,et al.  Kubios HRV - Heart rate variability analysis software , 2014, Comput. Methods Programs Biomed..

[50]  Qichen Li,et al.  Ventricular ectopic beat detection using a wavelet transform and a convolutional neural network , 2019, Physiological measurement.

[51]  L Tarassenko,et al.  Segmenting cardiac-related data using sleep stages increases separation between normal subjects and apnœic patients , 2004, Physiological measurement.

[52]  D. Eckberg Sympathovagal balance: a critical appraisal. , 1997, Circulation.

[53]  M. Malik,et al.  Phase-rectified signal averaging detects quasi-periodicities in non-stationary data , 2006 .

[54]  G.B. Moody,et al.  The impact of the MIT-BIH Arrhythmia Database , 2001, IEEE Engineering in Medicine and Biology Magazine.

[55]  Gari D. Clifford,et al.  Signal processing methods for heart rate variability , 2002 .

[56]  Gari D Clifford,et al.  Impact of the presence of noise on RR interval-based atrial fibrillation detection. , 2015, Journal of electrocardiology.

[57]  Patrick R. Norris,et al.  HEART RATE MULTISCALE ENTROPY AT THREE HOURS PREDICTS HOSPITAL MORTALITY IN 3,154 TRAUMA PATIENTS , 2008, Shock.

[58]  Yi-Lwun Ho,et al.  Multi-scale symbolic entropy analysis provides prognostic prediction in patients receiving extracorporeal life support , 2014, Critical Care.

[59]  Keith Willson,et al.  Relationship between detrended fluctuation analysis and spectral analysis of heart-rate variability. , 2002, Physiological measurement.