{"id":24826,"date":"2018-12-25T11:38:50","date_gmt":"2018-12-25T11:38:50","guid":{"rendered":"http:\/\/biomedpharmajournal.org\/?p=24826"},"modified":"2020-04-24T03:47:06","modified_gmt":"2020-04-24T03:47:06","slug":"comparative-analysis-of-heart-rate-variability-parameters-for-arrhythmia-and-atrial-fibrillation-using-anova","status":"publish","type":"post","link":"https:\/\/biomedpharmajournal.org\/staging\/vol11no4\/comparative-analysis-of-heart-rate-variability-parameters-for-arrhythmia-and-atrial-fibrillation-using-anova\/","title":{"rendered":"Comparative Analysis of Heart Rate Variability Parameters for Arrhythmia and Atrial Fibrillation using ANOVA"},"content":{"rendered":"<p><strong>Introduction<\/strong><\/p>\n<p>Electrocardiography (ECG) is a non-invasive method to extract electrical signals produced by Sino-Atrial (SA) node in the heart. By monitoring the ECG of a patient, prediction can be made which helps the doctor to take effective steps clinically.<sup>1<\/sup><sup>,<\/sup><sup>2<\/sup>\u00a0A significant decrease in the mortality rate of sudden heart attack cases can be achieved by regularly monitoring the ECG. ECG signal consists of several features delineating from <em>P<\/em>, <em>QRS<\/em> and <em>T<\/em> wave. Among all these features Heart Rate (HR) is an important feature to diagnose the cardiac.<sup>2<\/sup><sup>,<\/sup><sup>3<\/sup>\u00a0HR is calculated using the <em>RR<\/em> interval (the difference between two consecutive <em>R<\/em>-peaks). RR interval is the most prominent feature of the Heart Rate Variability (HRV) while other features are directly or indirectly dependent on it. HRV is a measure of fluctuation in between two heartbeats with respect to time.<\/p>\n<p>HRV is controlled by the primary part of our brain which is Autonomic Nervous System (ANS). It guides the working of HR, blood pressure, breathing and digestion despite from our desire. Nowadays, HRV is used to provide the information about our mental and physical health by observing one\u2019s lifestyle.<sup>4<\/sup><sup>,<\/sup><sup>5\u00a0<\/sup>HRV parameter assists clinicians to monitor sleep, stress, physical activity, emotion, thought, feelings etc. features on their patients. Low HRV signifies poor cardiac health that leads towards anxiety, high blood pressure and cardiac arrest, whereas high HRV represents a healthy heart. HRV components can be calculated in the time domain, geometrical domain and frequency domain.<sup>6<\/sup><\/p>\n<p>Various platforms are available nowadays for the collection of the database (offline and online).<sup>7<\/sup>\u00a0Firstly, pre-processing is done on the database to eliminate the undesired signal from the useful signal.<sup>8<\/sup><sup>,<\/sup><sup>9<\/sup>\u00a0Two types of unwanted signals are present which are Power Line Interference (PLI) and Baseline.<\/p>\n<p>Wander (BLW). To remove these noises, Finite Impulse Response (FIR) or Infinite Impulse Response (IIR) digital filters are used. IIR Butterworth filter is used by authors to diminish these noises.<sup>10<\/sup>\u00a0The same filter of 4th order is used by R. J. Ellis <em>et al<\/em>. for pre-processing.<sup>11<\/sup><\/p>\n<p>Butterworth filter is an IIR filter and it is not hardware implementable due to its unstable nature. A. M. A. Hassan used FIR filter to remove these unexpected signals due to ease of its hardware implementable nature.<sup>12<\/sup>\u00a0Secondly, feature extraction is done by F. Shaffer &amp; J. Ginsberg to extract the parameter of HRV for different duration and domains.<sup>13<\/sup>\u00a0There are various statistical tests which are used for the comparison in between and within groups like <em>t<\/em>-test, <em>f<\/em>-test, ANOVA etc. If the comparison is in between two groups then the <em>t<\/em>-test or <em>f<\/em>-test is used, while if more than two groups are there then ANOVA test is preferred.<sup>14<\/sup>\u00a0If we use the <em>t<\/em>-test or <em>f<\/em>-test on more than two groups, Type-I error will be introduced, which is unacceptable during statistical test performance. ANOVA test is a statistical test used for comparison or finding a relationship between or within features. Z. Qin <em>et al<\/em>. extract time and frequency domain features of HRV for the stress level evaluation.<sup>15<\/sup>\u00a0He also executed the ANOVA test for the statistical comparison among different level of stress. S. Agarwal and A. Wadhwani use ANOVA test on HRV features for the evaluation between meditative state and non-meditative state of a person.<sup>4<\/sup>\u00a0In this paper we have used three online databases: MIT\/BIH Normal Sinus Rhythm (NSR), MIT\/BIH Arrhythmia (AR) and MIT\/BIH Atrial Fibrillation (AF). From these databases of ECGs, removal of artifact is done by using digital FIR Kaiser Window filter because it has better frequency response than other conventional windows. Later from the filtered ECG signal, we have extracted time domain, geometrical domain and frequency domain HRV features. In the end, comparison of different databases on the basis of best features is done by using ANOVA test.The paper is structured as: Section 2 describes the different processing steps to be applied on different databases like MIT\/BIH NSR, MIT\/BIH AR and MIT\/BIH, Section 3 discusses the results of the different processing steps which is concluded in the last section.<\/p>\n<p><strong>Methodology <\/strong><\/p>\n<p>The process for the statistical analysis of feature processing consists primarily of pre-processing, feature extraction, feature selection and classification phase.<sup>16-18<\/sup>\u00a0In this paper, we are concentrating on the first three steps.<sup>19<\/sup>\u00a0Figure 1 illustrates the block diagram of HRV feature processing.<\/p>\n<p>&nbsp;<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24829\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig1-150x150.jpg\" alt=\"Figure 1: Block diagram of methodology of HRV Feature processing\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig1-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig1-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig1.jpg 481w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 1: Block diagram of methodology of HRV Feature processing<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig1.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p><strong>Database Collection<\/strong><\/p>\n<p>The online platform of Physionet is used for database collection. These databases are MIT\/BIH Normal Sinus Rhythm (NSR), MIT\/BIH Arrhythmia (AR) and MIT\/BIH Atrial Fibrillation (AF) with 18 records each.<sup>7<\/sup>\u00a0The sampling frequency of NSR, AR and AF database are 128 Hz, 360 Hz, and 250 Hz respectively. The database contains measurable data of the interval scale. Which helps in apply statistical parameter.<\/p>\n<p><strong>ECG Pre-processing<\/strong><\/p>\n<p>Several noises in ECG signals are introduced during the acquisition of the signal from the patient body. Among them, Power Line Interference (PLI) and Base Line Wander (BLW) are of main concern.<sup>20<\/sup>\u00a0These two unwanted signals create hindrance while extracting the useful information from the ECG wave. Pre-processing is done using a digital filter to diminish these undesirable signals.<sup>21<\/sup>\u00a0In this paper, FIR digital filters are used due to its simplicity in hardware implementation. Among various FIR filters, FIR Kaiser Window is used because it gives better frequency response as compared to other windowing technique. Useful bandwidth of ECG signal lies in between 0.5 Hz to 50 Hz. PLI is a high-frequency noise introduced due to the biomedical instrument used for the acquisition of ECG signal. PLI is eliminated using FIR Kaiser Low Pass Filter (LPF) having cut off frequency (\u2018<em>f<sub>c<\/sub><\/em>\u2019) of 50 Hz. Respiration during signal acquisition is the main cause of BLW containing low-frequency noise. Removal of BLW can be done be using FIR Kaiser High Pass Filter (HPF), having \u2018<em>f<sub>c<\/sub><\/em>\u2019 of 0.5 Hz. Later, we have obtained a denoised\/filtered signal which is used for HRV parameter extraction.<\/p>\n<p><strong><em>R<\/em><\/strong><strong>-peak Detection<\/strong><\/p>\n<p>ECG signal consists of mainly three waves <em>P<\/em>, <em>QRS<\/em> and <em>T<\/em>. Among these <em>QRS<\/em> wave is the prime focus, which can be further divided into three sub-waves <em>Q<\/em>, <em>R<\/em> and <em>S<\/em>. <em>R <\/em>wave is the most prominent feature of <em>QRS<\/em> due to its high positive peak. In this paper, the Pan-Tompkins Algorithm for <em>QRS<\/em> and <em>R<\/em> peak detection is used.<sup>22<\/sup>\u00a0For <em>QRS<\/em> detection, squaring and moving window integration is used. Since <em>R<\/em> wave is a positive peak, so squaring is done to eliminate all the negative peaks. Adaptive Thresholding is used for extraction of the <em>R<\/em>-peak from <em>QRS<\/em> wave after eliminating abnormal beats.<\/p>\n<p><strong>HRV Feature Extraction<\/strong><\/p>\n<p>Feature extraction from the <em>R<\/em>-peak of a filtered signal is the essential part for HRV feature processing.<sup>23<\/sup>\u00a0HRV features are divided into three domains: time domain, geometric domain and frequency domain. Spectrum for <em>RR<\/em> interval in time and frequency domain is demonstrated in Figure 2. Table 1 represents the explanation of HRV features of different domains with \u00a0their related formulas.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24830\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig2-150x150.jpg\" alt=\"Figure 2: Spectrum of RR interval (a)Time Domain(b) Frequency Domain\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig2-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig2.jpg 570w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td>&nbsp;<\/p>\n<p><strong>Figure 2: Spectrum of <em>RR <\/em>interval (<em>a<\/em>)Time Domain(<em>b<\/em>) Frequency Domain<\/strong><\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig2.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p><strong>HRV Feature Selection<\/strong><\/p>\n<p>A number of HRV features are present in various domains as defined but only some of them shows a noteworthy effect. That\u2019s why correlation among different HRV features of various domains expressed at two significance levels (95% and 99%) is shown in this paper. HRV parameters which show high correlation factor get selected while others are discarded. For the comparison of three groups of ECG databases, ANOVA test is used in different domains. ANOVA test demonstrates the variation among the mean of HRV parameters for different databases.<\/p>\n<p><strong>Result and Discussion<\/strong><\/p>\n<p>The input signal of ECG from different databases is used in this work which is demonstrated in figure 3. Fig. 3a represents the NSR ECG wave, Fig. 3b shows the AR ECG wave and Fig. 3c illustrates the AF ECG wave.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24831\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Tab1-150x150.jpg\" alt=\"Table 1: HRV measures of time, geometrical and frequency domain2\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Tab1-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Tab1-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Tab1.jpg 620w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Table 1: HRV measures of time, geometrical and frequency domain<sup>2<\/sup><\/strong><\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Tab1.jpg\" target=\"_blank\">Click here to View\u00a0table<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24836\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig3-150x150.jpg\" alt=\"Figure 3: ECG signal of (a) Normal Sinus Rhythm ECG (b) Arrhythmia ECG (c) Atrial Fibrillation\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig3-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig3-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig3.jpg 482w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 3: ECG signal of (<em>a<\/em>) Normal Sinus Rhythm ECG (<em>b<\/em>) Arrhythmia ECG (<em>c<\/em>) Atrial Fibrillation<\/strong><\/p>\n<p><strong><br \/>\n<\/strong><\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig3.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>In Pre-processing stage FIR Kaiser Window is used to denoise the signal from unwanted noises like PLI and BLW using LPF and HPF respectively. Fig. 4a shows the removal of PLI using LPF and Fig. 4b represents the removal of BLW using HPF.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24835\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig41-150x150.jpg\" alt=\"Figure 4: Filtered ECG signal using FIR Kaiser Window(a) LPF for PLI removal (b) HPF for BLW removal\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig41-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig41.jpg 585w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 4: Filtered ECG signal using FIR Kaiser Window(<em>a<\/em>) LPF for PLI removal (<em>b<\/em>) HPF for BLW removal<\/strong><\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig41.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>After getting filtered output, the Pan-Tompkins algorithm is used for real-time <em>QRS<\/em> detection. This algorithm is used for extracting the <em>R<\/em>-peak of <em>QRS <\/em>wave. <em>R<\/em>-peak is the most important feature to calculate various HRV parameters. Figure 5 demonstrates the <em>R<\/em>-peak from different <em>QRS<\/em> waves.<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td>\u00a0<img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-24837\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig51-150x150.jpg\" alt=\"Figure 5: R-peak detection\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig51-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig51.jpg 289w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/td>\n<td><strong>Figure 5: R-peak detection<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><a href=\"http:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2018\/12\/Vol11No4_Com_Kir_Fig51.jpg\" target=\"_blank\">Click here to View figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>Detection of <em>R<\/em>-peak helps in calculating various HRV features for different domains as described in Table 1. Time domain HRV feature describes the temporal variability of heart rate whereas the frequency domain of HRV feature defines the frequency distribution in a different frequency range. The geometric distribution gives the basic measurement to calculate the width of the <em>RR<\/em> interval histogram. Before doing any statistical test for comparing, the data must be normalized ( the \u2018<em>p<\/em>\u2019 value should be greater than 0.05). After a normality check, ANOVA test for comparison is used.<\/p>\n<p>In the ANOVA test, hypothesis testing is done on HRV features for different databases as mentioned earlier. In this paper, we have considered two hypothesis, null hypothesis (<em>H<\/em><sub>o<\/sub>) and alternative hypothesis (<em>H<\/em><sub>1<\/sub>).\u00a0 <em>H<\/em><sub>o <\/sub>means no variation and <em>H<\/em><sub>1<\/sub> signifies at least one of the mean is unequal.<sup>11<\/sup>\u00a0For ANOVA average case<\/p>\n<p><em>H<\/em><sub>o<\/sub>:<em>\u00b5<\/em><sub>1<\/sub>=<em>\u00b5<\/em><sub>2<\/sub>=<em>\u00b5<\/em><sub>3<\/sub><\/p>\n<p><em>H<\/em><sub>1<\/sub>=at least one of the mean is unequal<\/p>\n<p>where, <em>\u00b5<\/em><sub>1<\/sub> = Mean of NSR HRV feature<\/p>\n<p><em>\u00b5<\/em><sub>2<\/sub> = Mean of AR HRV feature<\/p>\n<p><em>\u00b5<\/em><sub>3<\/sub> = Mean of AF HRV feature<\/p>\n<p>Table 2 depicts the comparative results of the ANOVA test, which represents the significant variation among parameters of different datasets. The table illustrates the value in Mean \u00b1 Standard Deviation (SD) form. On the basis <em>p<\/em>-value significance level can be illustrated as.<\/p>\n<p><em>p<\/em>&gt; 0.05 = Weak Significance (WS)<\/p>\n<p>0.05&lt;<em>p<\/em> \u2264 0.10 = Moderately Significance (MS)<\/p>\n<p><em>p<\/em> \u2264 0.05 = Strongly Significance (SS)<\/p>\n<p><strong>Table 2: Comparison of HRV features using ANOVA<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"132\"><strong>Domain<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"109\"><strong>HRV Features<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\"><strong>\u00a0<\/strong><\/p>\n<p><strong>MIT\/BIH NSR<\/strong><\/p>\n<p><strong>Database<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"102\"><strong>\u00a0<\/strong><\/p>\n<p><strong>MIT\/BIH AR\u00a0 Database<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"104\"><strong>\u00a0<\/strong><\/p>\n<p><strong>MIT\/BIH AF<\/strong><\/p>\n<p><strong>Database<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"104\"><strong><em>\u00a0<\/em><\/strong><\/p>\n<p><strong><em>p<\/em><\/strong><strong>-value<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" rowspan=\"8\" width=\"132\"><strong>\u00a0<\/strong><\/p>\n<p><strong>Time<\/strong><\/p>\n<p><strong>Domain<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>mean RR<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">6.63E2\u00b16.72<\/td>\n<td style=\"text-align: center;\" width=\"102\">6.71E2\u00b114.40<\/td>\n<td style=\"text-align: center;\" width=\"104\">6.77E2\u00b125.96<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>SDNN<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">1.44E2\u00b13.89<\/td>\n<td style=\"text-align: center;\" width=\"102\">1.42E2\u00b117.20<\/td>\n<td style=\"text-align: center;\" width=\"104\">1.46E2\u00b19.17<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>SDSD<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">1.98E2\u00b15.53<\/td>\n<td style=\"text-align: center;\" width=\"102\">1.84E2\u00b120.76<\/td>\n<td style=\"text-align: center;\" width=\"104\">1.97E2\u00b113.47<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>RMSSD<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">1.98E2\u00b15.64<\/td>\n<td style=\"text-align: center;\" width=\"102\">1.84E2\u00b120.75<\/td>\n<td style=\"text-align: center;\" width=\"104\">1.97E2\u00b113.47<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>NN<\/em><\/strong><strong>50<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">7.88E2\u00b119.23<\/td>\n<td style=\"text-align: center;\" width=\"102\">5.37E2\u00b128.54<\/td>\n<td style=\"text-align: center;\" width=\"104\">8.09E2\u00b150.45<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p&lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>pNN<\/em><\/strong><strong>50<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">52.28\u00b10.89<\/td>\n<td style=\"text-align: center;\" width=\"102\">55.51\u00b12.65<\/td>\n<td style=\"text-align: center;\" width=\"104\">54.80\u00b13.12<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>Max-Min HR<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">61.11\u00b10.76<\/td>\n<td style=\"text-align: center;\" width=\"102\">61.78\u00b11.00<\/td>\n<td style=\"text-align: center;\" width=\"104\">61.28\u00b11.18<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>mean HR<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">97.05\u00b11.83<\/td>\n<td style=\"text-align: center;\" width=\"102\">93.78\u00b14.05<\/td>\n<td style=\"text-align: center;\" width=\"104\">94.89\u00b15.32<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>0.01&lt; p \u2264<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"132\"><strong>Geometrical Domain<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>HRVTI<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">22.24\u00b13.63<\/td>\n<td style=\"text-align: center;\" width=\"102\">26.92\u00b15.23<\/td>\n<td style=\"text-align: center;\" width=\"104\">21.49\u00b16.50<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.01<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" rowspan=\"11\" width=\"132\"><strong>Frequency<\/strong><\/p>\n<p><strong>Domain<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>aULF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">0.07\u00b10.008<\/td>\n<td style=\"text-align: center;\" width=\"102\">0.073\u00b10.007<\/td>\n<td style=\"text-align: center;\" width=\"104\">0.76\u00b10.008<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>aVLF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">1.34\u00b10.03<\/td>\n<td style=\"text-align: center;\" width=\"102\">1.33\u00b10.058<\/td>\n<td style=\"text-align: center;\" width=\"104\">1.37\u00b10.041<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p\u00a0 &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>aLF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">4.08\u00b10.60<\/td>\n<td style=\"text-align: center;\" width=\"102\">4.01\u00b10.11<\/td>\n<td style=\"text-align: center;\" width=\"104\">3.95\u00b10.66<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>aHF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">9.34\u00b10.082<\/td>\n<td style=\"text-align: center;\" width=\"102\">9.13\u00b10.25<\/td>\n<td style=\"text-align: center;\" width=\"104\">8.89\u00b11.90<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p\u00a0 &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>aTotal<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">14.76\u00b10.14<\/td>\n<td style=\"text-align: center;\" width=\"102\">14.47\u00b10.37<\/td>\n<td style=\"text-align: center;\" width=\"104\">14.81\u00b10.29<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &lt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>pVLF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">9.07\u00b10.18<\/td>\n<td style=\"text-align: center;\" width=\"102\">9.17\u00b10.28<\/td>\n<td style=\"text-align: center;\" width=\"104\">9.27\u00b10.30<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>pLF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">27.64\u00b10.27<\/td>\n<td style=\"text-align: center;\" width=\"102\">27.72\u00b10.47<\/td>\n<td style=\"text-align: center;\" width=\"104\">27.72\u00b10.36<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>pHF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">63.28\u00b10.21<\/td>\n<td style=\"text-align: center;\" width=\"102\">63.04\u00b10.58<\/td>\n<td style=\"text-align: center;\" width=\"104\">63.06\u00b10.57<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>Lfnorm<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">0.30\u00b10.002<\/td>\n<td style=\"text-align: center;\" width=\"102\">0.30\u00b10.005<\/td>\n<td style=\"text-align: center;\" width=\"104\">0.30\u00b10.004<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>Hfnorm<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">0.70\u00b10.002<\/td>\n<td style=\"text-align: center;\" width=\"102\">0.69\u00b10.005<\/td>\n<td style=\"text-align: center;\" width=\"104\">0.69\u00b10.004<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"109\"><strong><em>LF\/HF<\/em><\/strong><\/td>\n<td style=\"text-align: center;\" width=\"122\">0.44\u00b10.005<\/td>\n<td style=\"text-align: center;\" width=\"102\">0.44\u00b10.01<\/td>\n<td style=\"text-align: center;\" width=\"104\">0.44\u00b10.009<\/td>\n<td style=\"text-align: center;\" width=\"104\"><em>p &gt;<\/em>0.05<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<p>The feature values can be affected by the sampling frequency and waveform duration of the database taken while processing. In this paper, we have generated the results at 1000 Hz sampling frequency for a long-term database. After observing the results of ANOVA test for these specifications, it is interpreted that some HRV features like SDSD, RMSSD, NN50, pNN50, HRVTI, aULF, aVLF, aHF, and aTotal gives a strong indication to reject H<sub>o<\/sub>. It represents that these features show a significant variation in their mean value. While, all other parameters show weak evidence against to reject H<sub>o<\/sub>, except mean HR which represents moderate significant change. Mostly HRV feature represents a strong significant change which denotes that there is a substantial difference occurring among different databases. After analyzing the p-value of ANOVA which signifies the noticeable variation among different databases, we have selected the best HRV features. Best selected features are SDSD, RMSSD, NN50, pNN50, mean HR aULF, aVLF, aHF and aTotal which we will be used for further classification in the future.<\/p>\n<p><strong>Conclusion<\/strong><\/p>\n<p>This paper presents the statistical comparison among HRV features of different domains (time, geometrical and frequency) and databases (NSR, AR and AF). In the pre-processing stage FIR Kaiser digital filter is used to eliminate noise from different database. Then, twenty HRV features are extracted from different domain &amp; feature selection is done with correlation test. For the comparison between three databases, ANOVA test has been used to yield the best features. Only 9 out of 20 features show noteworthy difference between three databases. In future, these selected features will be used for classification. Later, hardware implementation of HRV feature processing will be done using VIVADO tool on FPGA.<\/p>\n<p><strong>References<\/strong><\/p>\n<ol>\n<li>Parashar N., Jain S., Sood M and Dogra J. Review of biomedical system for high performance applications. \u00a0<em>International Conference on Signal Processing. Computing and Control.\u00a0<\/em>2017;300-304<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/ISPCC.2017.8269693\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Shaffer F and Ginsberg J. An overview of heart rate variability metrics and norms.\u00a0<em>Frontiers in public health.<\/em> 2017;5:258.<br \/>\n<a href=\"https:\/\/doi.org\/10.3389\/fpubh.2017.00258\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Li Y., Tang X.,Xu Z and Yan H. A novel approach to phase space reconstruction of single lead ECG for QRS complex detection. <em>Biomedical Signal Processing and Control. <\/em>2018;39:405-415.<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.bspc.2017.06.007\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Agarwal S and Wadhwani A. Analysis of Heart Rate Variability During Meditative and Non-Meditative State Using Analysis of Variance. <em>International Journal of Advanced Biological and Biomedical Research. <\/em>2013;1:728-736.<\/li>\n<li>Parashar N., Jain S., Sood M. Semiautomatic Detection of Cardiac Diseases employing Dual Tree Complex Wavelet Transform. <em>Periodicals of Engineering and Natural Sciences (PEN). <\/em>2018;6:129-140.<\/li>\n<li>Dhiman A., Dubey S.,Jain S. Design of Lead II ECG Waveform and Classification Performance for Morphological features using Different Classifiers on Lead II.<em> Research Journal of Pharmaceutical. Biological and Chemical Sciences.\u00a0<\/em>2016;7:1226-1231.<\/li>\n<li>Goldberger A. L., Amaral A., Glass L.,Hausdorff J. M., Ivanov P. C., Mark R. G.<em>, et al.<\/em> PhysioBank, PhysioToolkit and PhysioNet: components of a new research resource for complex physiologic signals <em>Circulation.<\/em> 2000;101:e215-e220.<br \/>\n<a href=\"https:\/\/doi.org\/10.1161\/01.CIR.101.23.e215\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Acharya U. R., Fujita H., Adam M., Lih O. S., Sudarshan V. K.,Hong T. J.<em>, et al.<\/em>, Automated characterization and classification of coronary artery disease and myocardial infarction by decomposition of ECG signals: A comparative study. <em>Information Sciences.<\/em> 2017;377:17-29.<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.ins.2016.10.013\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Jain S. Implementation of marker proteins using standardised effect. <em>Journal of Global Pharma Technology.\u00a0<\/em>2017.<\/li>\n<li>Kambo T.,Avtarjaswal R. De-noising and statistical feature extraction of the ecg signal using wavelet analysis. <em>International Journal Of Electrical, Electronics And Data Communication.\u00a0<\/em>2016;4<\/li>\n<li>Ellis R. J., Zhu B., Koenig J., Thayer J. F and Wang Y. A careful look at ECG sampling frequency and R-peak interpolation on short-term measures of heart rate variability. <em>Physiological measurement.\u00a0<\/em>2015;36:1827.<br \/>\n<a href=\"https:\/\/doi.org\/10.1088\/0967-3334\/36\/9\/1827\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Hassan A. M. A. FPGA Realization for Baseline Wander Noise Cancellation of ECG Signals using Wavelet Transform. <em>International Journal of Computer Applications. <\/em>2017;168.<\/li>\n<li>Urooj M., Pillai K.,Tandon M., Venkateshan S and Saha N. Reference ranges for time domain parameters of heart rate variability in indian population and validation in hypertensive subjects and smokers.<em>Int J Pharm Pharm Sci.\u00a0<\/em>2011;3:36-39.<\/li>\n<li>Jain S. Regression analysis on different mitogenic pathways. <em>Network Biology.\u00a0<\/em>2016;6:40.<\/li>\n<li>Qin Z., Li M., Huang L and Zhao Y. Stress level evaluation using BP Neural network based on time-frequency analysis of HRV. <em>IEEE International Conference on Mechatronics and Automation.\u00a0<\/em>2017;1798-1803.<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/ICMA.2017.8016090\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Bhusri S and Jain S. Analysis of breast lesions using laws&#8217; mask texture features. <em>Fourth International Conference on Parallel, Distributed and Grid Computing.\u00a0<\/em>2016;56-60.<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/PDGC.2016.7913115\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Jain S. Classification of Protein Kinase B using discrete wavelet transform. <em>International Journal of Information Technology.\u00a0<\/em>2018;10:211-216.<br \/>\n<a href=\"https:\/\/doi.org\/10.1007\/s41870-018-0090-7\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Rana S., Jain S and Virmani J. Classification of focal kidney lesions using wavelet-based texture descriptors. <em>International Journal of Pharma and Bio Sciences.\u00a0<\/em>2016;7:646-652.<\/li>\n<li>Sharma S., Jain S and Bhusri S. Two Class Classification of Breast Lesions using Statistical and Transform Domain features. <em>Journal Of Global Pharma Technology Methodology.\u00a0<\/em>2017;9:18-24.<\/li>\n<li>Parashar N., Jain S., Sood M. Removal of electromyography noise from ECG for high performance biomedical systems. <em>Network Biology.\u00a0<\/em>2018;8:12.<\/li>\n<li>Singh G and Prakash N. R. FPGA Implementation of Higher Order FIR Filter. <em>International Journal of Electrical and Computer Engineering.\u00a0<\/em>2017;7:1874-1881.<br \/>\n<a href=\"https:\/\/doi.org\/10.11591\/ijece.v7i4.pp1874-1881\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Pan J and Tompkins W. J. A real-time QRS detection algorithm. <em>IEEE Trans. Biomed. Eng. <\/em>1985;32:230-236.<br \/>\n<a href=\"https:\/\/doi.org\/10.1109\/TBME.1985.325532\" target=\"_blank\">CrossRef<\/a><\/li>\n<li>Gupta R.,Singh S., Garg K and Jain S. Indigenous Design of Electronic Circuit for Electrocardiograph. International Journal of Innovative Research in Science, Engineering and Technology. 2014;3.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Introduction Electrocardiography (ECG) is a non-invasive method to extract electrical  [&#8230;]<\/p>\n","protected":false},"author":9,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[61],"tags":[],"class_list":["post-24826","post","type-post","status-publish","format-standard","hentry","category-vol11no4"],"_links":{"self":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/24826","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/users\/9"}],"replies":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/comments?post=24826"}],"version-history":[{"count":6,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/24826\/revisions"}],"predecessor-version":[{"id":32557,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/24826\/revisions\/32557"}],"wp:attachment":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/media?parent=24826"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/categories?post=24826"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/tags?post=24826"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}