Cardiac state diagnosis using higher order spectra of heart rate variability

Transformations between rotational and translational invariants formulated in reciprocal spaces

Correlation functions play an important role in the theoretical underpinnings of many disparate areas of the physical sciences: in particular, scattering theory. More recently, they have become useful in the classification of objects in areas such as computer vision and our area of cryoEM. Our primary classification scheme in the cryoEM image processing system, EMAN2, is now based on third order invariants formulated in Fourier space. This allows a factor of 8 speed up in the two classification procedures inherent in our software pipeline, because it allows for classification without the need for computationally costly alignment procedures. In this work, we address several formal and practical aspects of such multispectral invariants. We show that we can formulate such invariants in the representation in which the original signal is most compact. We explicitly construct transformations between invariants in different orientations for arbitrary order of correlation functions and dimension. We demonstrate that third order invariants distinguish 2D mirrored patterns (unlike the radial power spectrum), which is a fundamental aspects of its classification efficacy. We show the limitations of 3rd order invariants also, by giving an example of a wide family of patterns with identical (vanishing) set of 3rd order invariants. For sufficiently rich patterns, the third order invariants should distinguish typical images, textures and patterns.

Nonlinear analysis of heart rate variability signals in meditative state: a review and perspective

INTRODUCTION

In recent times, an upsurge in the investigation related to the effects of meditation in reconditioning various cardiovascular and psychological disorders is seen. In majority of these studies, heart rate variability (HRV) signal is used, probably for its ease of acquisition and low cost. Although understanding the dynamical complexity of HRV is not an easy task, the advances in nonlinear analysis has significantly helped in analyzing the impact of meditation of heart regulations. In this review, we intend to present the various nonlinear approaches, scientific findings and their limitations to develop deeper insights to carry out further research on this topic.

RESULTS

Literature have shown that research focus on nonlinear domain is mainly concentrated on assessing predictability, fractality, and entropy-based dynamical complexity of HRV signal. Although there were some conflicting results, most of the studies observed a reduced dynamical complexity, reduced fractal dimension, and decimated long-range correlation behavior during meditation. However, techniques, such as multiscale entropy (MSE) and multifractal analysis (MFA) of HRV can be more effective in analyzing non-stationary HRV signal, which were hardly used in the existing research works on meditation.

CONCLUSIONS

After going through the literature, it is realized that there is a requirement of a more rigorous research to get consistent and new findings about the changes in HRV dynamics due to the practice of meditation. The lack of adequate standard open access database is a concern in drawing statistically reliable results. Albeit, data augmentation technique is an alternative option to deal with this problem, data from adequate number of subjects can be more effective. Multiscale entropy analysis is scantily employed in studying the effect of meditation, which probably need more attention along with multifractal analysis.

METHODS

Scientific databases, namely PubMed, Google Scholar, Web of Science, Scopus were searched to obtain the literature on "HRV analysis during meditation by nonlinear methods". Following an exclusion criteria, 26 articles were selected to carry out this scientific analysis.

Autonomic nervous system assessment using heart rate variability

The role of the autonomic nervous system in the onset of supraventricular and ventricular arrhythmias is well established. It can be analysed by the spontaneous behaviour of the heart rate with ambulatory ECG recordings, through heart rate variability measurements. Input of heart rate variability parameters into artificial intelligence models to make predictions regarding the detection or forecast of rhythm disorders is becoming routine and neuromodulation techniques are now increasingly used for their treatment. All this warrants a reappraisal of the use of heart rate variability for autonomic nervous system assessment.Measurements performed over long periods such as 24H-variance, total power, deceleration capacity, and turbulence are suitable for estimating the individual basal autonomic status. Spectral measurements performed over short periods provide information on the dynamics of systems that disrupt this basal balance and may be part of the triggers of arrhythmias, as well as premature atrial or ventricular beats. All heart rate variability measurements essentially reflect the modulations of the parasympathetic nervous system which are superimposed on the impulses of the adrenergic system. Although heart rate variability parameters have been shown to be useful for risk stratification in patients with myocardial infarction and patients with heart failure, they are not part of the criteria for prophylactic implantation of an intracardiac defibrillator, because of their high variability and the improved treatment of myocardial infarction. Graphical methods such as Poincaré plots allow quick screening of atrial fibrillation and are set to play an important role in the e-cardiology networks. Although mathematical and computational techniques allow manipulation of the ECG signal to extract information and permit their use in predictive models for individual cardiac risk stratification, their explicability remains difficult and making inferences about the activity of the ANS from these models must remain cautious.

Short-Term HRV Detection and Human Fatigue State Analysis Based on Optical Fiber Sensing Technology

Mental fatigue is a key cause of chronic diseases and traffic accidents, which is difficult to be quantitatively evaluated. In order to non-intrusively detect fatigue state, an optical fiber sensing system is proposed, which is non-invasive and does not require direct contact with skin. The fiber sensor was fabricated through phase mask exposure method and packaged by sensitivity-enhanced structure, which can suppress transverse force and increase signal amplitude by 5%. A fatigue-inducing experiment was carried out, and the heartbeat signals of 20 subjects under different fatigue states were collected by the proposed sensing system. A series of heart rate variability indicators were calculated from the sensing signals, and their statistical significance for fatigue was analyzed. The experiment results showed that the values of SDNN and LF/HF increased significantly with subjects' fatigue level. This study shows that the proposed fiber optic sensing system has practical value in fatigue state monitoring.

Bispectral Analysis of Heart Rate Variability to Characterize and Help Diagnose Pediatric Sleep Apnea

Pediatric obstructive sleep apnea (OSA) is a breathing disorder that alters heart rate variability (HRV) dynamics during sleep. HRV in children is commonly assessed through conventional spectral analysis. However, bispectral analysis provides both linearity and stationarity information and has not been applied to the assessment of HRV in pediatric OSA. Here, this work aimed to assess HRV using bispectral analysis in children with OSA for signal characterization and diagnostic purposes in two large pediatric databases (0-13 years). The first database (training set) was composed of 981 overnight ECG recordings obtained during polysomnography. The second database (test set) was a subset of the Childhood Adenotonsillectomy Trial database (757 children). We characterized three bispectral regions based on the classic HRV frequency ranges (very low frequency: 0-0.04 Hz; low frequency: 0.04-0.15 Hz; and high frequency: 0.15-0.40 Hz), as well as three OSA-specific frequency ranges obtained in recent studies (BW1: 0.001-0.005 Hz; BW2: 0.028-0.074 Hz; BWRes: a subject-adaptive respiratory region). In each region, up to 14 bispectral features were computed. The fast correlation-based filter was applied to the features obtained from the classic and OSA-specific regions, showing complementary information regarding OSA alterations in HRV. This information was then used to train multi-layer perceptron (MLP) neural networks aimed at automatically detecting pediatric OSA using three clinically defined severity classifiers. Both classic and OSA-specific MLP models showed high and similar accuracy (Acc) and areas under the receiver operating characteristic curve (AUCs) for moderate (classic regions: Acc = 81.0%, AUC = 0.774; OSA-specific regions: Acc = 81.0%, AUC = 0.791) and severe (classic regions: Acc = 91.7%, AUC = 0.847; OSA-specific regions: Acc = 89.3%, AUC = 0.841) OSA levels. Thus, the current findings highlight the usefulness of bispectral analysis on HRV to characterize and diagnose pediatric OSA.

Automated detection of conduct disorder and attention deficit hyperactivity disorder using decomposition and nonlinear techniques with EEG signals

BACKGROUND AND OBJECTIVES

Attention deficit hyperactivity disorder (ADHD) is often presented with conduct disorder (CD). There is currently no objective laboratory test or diagnostic method to discern between ADHD and CD, and diagnosis is further made difficult as ADHD is a common neuro-developmental disorder often presenting with other co-morbid difficulties; and in particular with conduct disorder which has a high degree of associated behavioural challenges. A novel automated system (AS) is proposed as a convenient supplementary tool to support clinicians in their diagnostic decisions. To the best of our knowledge, we are the first group to develop an automated classification system to classify ADHD, CD and ADHD+CD classes using brain signals.

METHODS

The empirical mode decomposition (EMD) and discrete wavelet transform (DWT) methods were employed to decompose the electroencephalogram (EEG) signals. Autoregressive modelling coefficients and relative wavelet energy were then computed on the signals. Various nonlinear features were extracted from the decomposed coefficients. Adaptive synthetic sampling (ADASYN) was then employed to balance the dataset. The significant features were selected using sequential forward selection method. The highly discriminatory features were subsequently fed to an array of classifiers.

RESULTS

The highest accuracy of 97.88% was achieved with the K-Nearest Neighbour (KNN) classifier. The proposed system was developed using ten-fold validation strategy on EEG data from 123 children. To the best of our knowledge this is the first study to develop an AS for the classification of ADHD, CD and ADHD+CD classes using EEG signals.

POTENTIAL APPLICATION

Our AS can potentially be used as a web-based application with cloud system to aid the clinical diagnosis of ADHD and/or CD, thus supporting faster and accurate treatment for the children. It is important to note that testing with larger data is required before the AS can be employed for clinical applications.

Study of the impact of clicks and murmurs on cardiac sounds S1 and S2 through bispectral analysis

This paper presents a study of the impact of clicks, and murmurs on cardiac sound S1, and S2, and the measure of severity degree through synchronization degree between frequencies, using bispectral analysis. The algorithm is applied on three groups of Phonocardiogram (PCG) signal: group A represents PCG signals having a morphology similar to that of the normal PCG signal without click or murmur, group B represents PCG signals with a click (reduced murmur), and group C represent PCG signals with murmurs.

The proposed algorithm permits us to evaluate and quantify the relationship between the two sounds S1 and S2 on one hand and between the two sounds, click and murmur on the other hand. The obtained results show that the clicks and murmurs can affect both the heart sounds, and vice versa. This study shows that the heart works in perfect harmony and that the frequencies of sounds S1, S2, clicks, and murmurs are not accidentally generated; but they are generated by the same generator system. It might also suggest that one of the obtained frequencies causes the others.

The proposed algorithm permits us also to determine the synchronization degree. It shows high values in group C; indicating high severity degrees, low values for group B, and zero in group A.

The algorithm is compared to Short-Time Fourier Transform (STFT) and continuous wavelet transform (CWT) analysis. Although the STFT can provide correctly the time, it can’t distinguish between the internal components of sounds S1 and S2, which are successfully determined by CWT, which, in turn, cannot find the relationship between them.

The algorithm was also evaluated and compared to the energetic ratio. the obtained results show very satisfactory results and very good discrimination between the three groups.

We can conclude that the three algorithms (STFT, CWT, and bispectral analysis) are complementary to facilitate a good approach and to better understand the cardiac sounds

Detecting Congestive Heart Failure by Extracting Multimodal Features and Employing Machine Learning Techniques

The adaptability of heart to external and internal stimuli is reflected by the heart rate variability (HRV). Reduced HRV can be a predictor of negative cardiovascular outcomes. Based on the nonlinear, nonstationary, and highly complex dynamics of the controlling mechanism of the cardiovascular system, linear HRV measures have limited capability to accurately analyze the underlying dynamics. In this study, we propose an automated system to analyze HRV signals by extracting multimodal features to capture temporal, spectral, and complex dynamics. Robust machine learning techniques, such as support vector machine (SVM) with its kernel (linear, Gaussian, radial base function, and polynomial), decision tree (DT), k-nearest neighbor (KNN), and ensemble classifiers, were employed to evaluate the detection performance. Performance was evaluated in terms of specificity, sensitivity, positive predictive value (PPV), negative predictive value (NPV), and area under the receiver operating characteristic curve (AUC). The highest performance was obtained using SVM linear kernel (TA = 93.1%, AUC = 0.97, 95% CI [lower bound = 0.04, upper bound = 0.89]), followed by ensemble subspace discriminant (TA = 91.4%, AUC = 0.96, 95% CI [lower bound 0.07, upper bound = 0.81]) and SVM medium Gaussian kernel (TA = 90.5%, AUC = 0.95, 95% CI [lower bound = 0.07, upper bound = 0.86]). The results reveal that the proposed approach can provide an effective and computationally efficient tool for automatic detection of congestive heart failure patients.

Comparison of methods for the assessment of nonlinearity in short-term heart rate variability under different physiopathological states

Despite the widespread diffusion of nonlinear methods for heart rate variability (HRV) analysis, the presence and the extent to which nonlinear dynamics contribute to short-term HRV are still controversial. This work aims at testing the hypothesis that different types of nonlinearity can be observed in HRV depending on the method adopted and on the physiopathological state. Two entropy-based measures of time series complexity (normalized complexity index, NCI) and regularity (information storage, IS), and a measure quantifying deviations from linear correlations in a time series (Gaussian linear contrast, GLC), are applied to short HRV recordings obtained in young (Y) and old (O) healthy subjects and in myocardial infarction (MI) patients monitored in the resting supine position and in the upright position reached through head-up tilt. The method of surrogate data is employed to detect the presence and quantify the contribution of nonlinear dynamics to HRV. We find that the three measures differ both in their variations across groups and conditions and in the percentage and strength of nonlinear HRV dynamics. NCI and IS displayed opposite variations, suggesting more complex dynamics in O and MI compared to Y and less complex dynamics during tilt. The strength of nonlinear dynamics is reduced by tilt using all measures in Y, while only GLC detects a significant strengthening of such dynamics in MI. A large percentage of detected nonlinear dynamics is revealed only by the IS measure in the Y group at rest, with a decrease in O and MI and during T, while NCI and GLC detect lower percentages in all groups and conditions. While these results suggest that distinct dynamic structures may lie beneath short-term HRV in different physiological states and pathological conditions, the strong dependence on the measure adopted and on their implementation suggests that physiological interpretations should be provided with caution.

Spectral Analysis of Heart Rate Variability: Time Window Matters

Spectral analysis of heart rate variability (HRV) is a valuable tool for the assessment of cardiovascular autonomic function. Fast Fourier transform and autoregressive based spectral analysis are two most commonly used approaches for HRV analysis, while new techniques such as trigonometric regressive spectral (TRS) and wavelet transform have been developed. Short-term (on ECG of several minutes) and long-term (typically on ECG of 1–24 h) HRV analyses have different advantages and disadvantages. This article reviews the characteristics of spectral HRV studies using different lengths of time windows. Short-term HRV analysis is a convenient method for the estimation of autonomic status, and can track dynamic changes of cardiac autonomic function within minutes. Long-term HRV analysis is a stable tool for assessing autonomic function, describe the autonomic function change over hours or even longer time spans, and can reliably predict prognosis. The choice of appropriate time window is essential for research of autonomic function using spectral HRV analysis.

Bispectrum and Histogram Features for the Identification of Atrial Fibrillation Based on Electrocardiogram

Atrial Fibrillation (AF) is probably the most common serious abnormal heart rhythm. It affects about 2% to 3% of the population in Europe and North America. In this study, we proposed an effective Atrial Fibrillation (AF) identification system based on RR interval (RRI) analysis. Two preprocessing methods were employed to remove the motion artifacts and ectopic beats. Three categories of RRI features, including base, bispectrum, and histogram features, were proposed to enhance the performance of the identifier. The roles of different feature categories were evaluated. The combination of the three categories of features were demonstrated to compensate with one another to construct an effective feature set for AF identification. When compared to other representative AF identifiers in the literature, the proposed method outperforms them with superior recognition rates by using much larger number of testing data.

Higher-order spectral analysis of spontaneous speech signals in Alzheimer's disease

An early and accurate diagnosis of Alzheimer's disease (AD) has been progressively attracting more attention in recent years. One of the main problems of AD is the loss of language skills. This paper presents a computational framework for classifying AD patients compared to healthy control subjects using information from spontaneous speech signals. Spontaneous speech data are obtained from 30 AD patients and 30 healthy controls. Because of the nonlinear and dynamic nature of speech signals, higher order spectral features (specifically bispectrum) were used for analysis. Four classifiers (k-Nearest Neighbor, Support Vector Machine, Naïve Bayes and Decision tree) were used to classify subjects into three different levels of AD and healthy group based on their performance in terms of the HOS-based features. Ten-fold cross-validation method was used to test the reliability of the classifier results. The results showed that the proposed method had a good potential in AD diagnosis. The proposed method was also able to diagnose the earliest stage of AD with high accuracy. The method has the great advantage of being non-invasive, cost-effective, and associated with no side effects. Therefore, the proposed method can be a spontaneous speech directed test for pre-clinical evaluation of AD diagnosis.

Bispectral analysis of spontaneous EEG activity from patients with moderate dementia due to Alzheimer's disease

Dementia due to Alzheimer's disease (AD) is a common disorder with a great impact on the patients' quality of life. The aim of this pilot study was to characterize spontaneous electroencephalography (EEG) activity in dementia due to AD using bispectral analysis. Five minutes of EEG activity were recorded from 17 patients with moderate dementia due to AD and 19 age-matched controls. Bispectrum results revealed that AD patients are characterized by an increase of phase coupling at low frequencies in comparison with controls. Additionally, some bispectral features calculated from the bispectrum showed significant differences between both groups (p <; 0.05, Mann-Whitney U test with Bonferroni's correction). Finally, a stepwise logistic regression analysis with a leave-one-out cross-validation procedure was used for classification purposes. An accuracy of 86.11% (sensitivity = 88.24%; specificity =84.21%) was achieved. This study suggests the usefulness of bispectral analysis to provide further insights into the underlying brain dynamics associated with AD.

ASSESSMENT OF CARDIAC HEART FAILURE AND CARDIAC ARTERY DISEASE BY THE HIGHER ORDER SPECTRA

Cardiac diseases are major reason of death in the world populace and the numeral of cases is upsurging every year. Due to cardiac artery disease (CAD), the strength of heart muscles becomes weak and heart pumping is disturbed which may eventually lead to abnormal heart beat and heart failure. Therefore, the beginning stage detection of CAD and cardiac heart failure (CHF) are of prime importance. In this work, we have used a non-invasive diagnosis method as higher order spectra (HOS) for assessment of cardiac diseases. The method indicates whether or not a cardiac heart disease is present, by assessing the cardiac health of subjects using extracted features from heart rate variability (HRV) signals. This assessment is based on 10 spectra nonlinear features. These features were extracted from HRV signals by using the HOS method. For this study, the R-R interval data (i.e. HRV signals) were taken from the standard database of cardiac heart failure (CHF), CAD patients, healthy young (YNG) and Self recorded of healthy young (SELF_YNG) subjects. Statistical assessments were performed on the group of database sets as YNG-CAD, YNG-CHF, SELF_YNG-CAD and Self_YNG-CHF subjects. A Wilcoxon rank sum test ([Formula: see text]-value) was used to statistically compare the features extracted by HOS for group of data sets. It indicates whether or not the same features of individual classes of HRV data sets are dissimilar. The results depicted that the all features are very significant ([Formula: see text]) except the phase entropy (PHE) feature which is not significant for CAD-CHF, SELF_YNG-CAD and SELF_YNG-CHF group of subjects. While in the case of YNG-CAD group of subjects, features like first-order spectral moment of amplitudes of diagonal elements (H3), PHE and logarithmic amplitudes of diagonal elements (H2) are significant ([Formula: see text]) and excluding these features, the remaining features are very significant except MM and H1 which are not significant. The results also depicted that the mean value of sum of logarithmic amplitude (H1), H2, normalized entropy (P1), normalized squared entropy (P2) and PHE features of healthy YNG subjects are having higher values than that of CAD and CHF patients. While weighted center of bi-spectrum (WCOB2) and FLAT spectrum features are lower than CAD and CHF patients compared to YNG subjects. In case of CAD and CHF patients, all the features of CAD patients are having higher values compared to CHF except P1, P2 and WCOB1.

Utility of bispectrum in the screening of pediatric sleep apnea-hypopnea syndrome using oximetry recordings

BACKGROUND AND OBJECTIVE

The aim of this study was to assess the utility of bispectrum-based oximetry approaches as a complementary tool to traditional techniques in the screening of pediatric sleep apnea-hypopnea syndrome (SAHS).

METHODS

298 blood oxygen saturation (SpO₂) signals from children ranging 0-13 years of age were recorded during overnight polysomnography (PSG). These recordings were divided into three severity groups according to the PSG-derived apnea hypopnea index (AHI): AHI < 5 events per hour (e/h), 5 ≤ AHI < 10 e/h, AHI ≥ 10 e/h. For each pediatric subject, anthropometric variables, 3% oxygen desaturation index (ODI3) and spectral features from power spectral density (PSD) and bispectrum were obtained. Then, the fast correlation-based filter (FCBF) was applied to select a subset of relevant features that may be complementary, excluding those that are redundant. The selected features fed a multiclass multi-layer perceptron (MLP) neural network to build a model to estimate the SAHS severity degrees.

RESULTS

An optimum subset with features from all the proposed methodological approaches was obtained: variables from bispectrum, as well as PSD, ODI3, Age, and Sex. In the 3-class classification task, the MLP model trained with these features achieved an accuracy of 76.0% and a Cohen's kappa of 0.56 in an independent test set. Additionally, high accuracies were reached using the AHI cutoffs for diagnosis of moderate (AHI = 5 e/h) and severe (AHI = 10 e/h) SAHS: 81.3% and 85.3%, respectively. These results outperformed the diagnostic ability of a MLP model built without using bispectral features.

CONCLUSIONS

Our results suggest that bispectrum provides additional information to anthropometric variables, ODI3 and PSD regarding characterization of changes in the SpO₂ signal caused by respiratory events. Thus, oximetry bispectrum can be a useful tool to provide complementary information for screening of moderate-to-severe pediatric SAHS.

Optimized multi-level elongated quinary patterns for the assessment of thyroid nodules in ultrasound images

Ultrasound imaging is one of the most common visualizing tools used by radiologists to identify the location of thyroid nodules. However, visual assessment of nodules is difficult and often affected by inter- and intra-observer variabilities. Thus, a computer-aided diagnosis (CAD) system can be helpful to cross-verify the severity of nodules. This paper proposes a new CAD system to characterize thyroid nodules using optimized multi-level elongated quinary patterns. In this study, higher order spectral (HOS) entropy features extracted from these patterns appropriately distinguished benign and malignant nodules under particle swarm optimization (PSO) and support vector machine (SVM) frameworks. Our CAD algorithm achieved a maximum accuracy of 97.71% and 97.01% in private and public datasets respectively. The evaluation of this CAD system on both private and public datasets confirmed its effectiveness as a secondary tool in assisting radiological findings.

Reference-Free Assessment of Speech Intelligibility Using Bispectrum of an Auditory Neurogram

Sensorineural hearing loss occurs due to damage to the inner and outer hair cells of the peripheral auditory system. Hearing loss can cause decreases in audibility, dynamic range, frequency and temporal resolution of the auditory system, and all of these effects are known to affect speech intelligibility. In this study, a new reference-free speech intelligibility metric is proposed using 2-D neurograms constructed from the output of a computational model of the auditory periphery. The responses of the auditory-nerve fibers with a wide range of characteristic frequencies were simulated to construct neurograms. The features of the neurograms were extracted using third-order statistics referred to as bispectrum. The phase coupling of neurogram bispectrum provides a unique insight for the presence (or deficit) of supra-threshold nonlinearities beyond audibility for listeners with normal hearing (or hearing loss). The speech intelligibility scores predicted by the proposed method were compared to the behavioral scores for listeners with normal hearing and hearing loss both in quiet and under noisy background conditions. The results were also compared to the performance of some existing methods. The predicted results showed a good fit with a small error suggesting that the subjective scores can be estimated reliably using the proposed neural-response-based metric. The proposed metric also had a wide dynamic range, and the predicted scores were well-separated as a function of hearing loss. The proposed metric successfully captures the effects of hearing loss and supra-threshold nonlinearities on speech intelligibility. This metric could be applied to evaluate the performance of various speech-processing algorithms designed for hearing aids and cochlear implants.

APPLICATION OF EMPIRICAL MODE DECOMPOSITION–BASED FEATURES FOR ANALYSIS OF NORMAL AND CAD HEART RATE SIGNALS

Coronary Artery Disease (CAD) is a heart disease caused due to insufficient supply of nutrients and oxygen to the heart muscles. Hence, reduced supply of nutrients and oxygen causes heart attack or stroke and may cause death. Also significant number of people are suffering from CAD around the world so timely diagnosis of CAD can save the life of patients. In this work, we have proposed computer assisted diagnosis of CAD using Heart Rate (HR) signals obtained from Electrocardiogram (ECG) signals. We have used the Empirical Mode Decomposition (EMD) technique to process the HR signals. The features namely: Second-Order Difference Plot (SODP) area, Analytic Signal Representation (ASR) area, Amplitude Modulation (AM) bandwidth, Frequency Modulation (FM) bandwidth and Fourier–Bessel expansion (FBE)- based mean frequency computed from the Intrinsic Mode Functions (IMFs) are extracted to discriminate normal and CAD subjects. Thereafter, Kruskal–Wallis statistical test is performed on these features. The features having p-value less than 0.05 are considered to be significant. Our results show that three features namely: AM bandwidth, FM bandwidth and FBE-based mean frequency are more suitable than ASR area and SODP area features for discrimination of normal and CAD subjects.

An integrated index for automated detection of infarcted myocardium from cross-sectional echocardiograms using texton-based features (Part 1)

Cross-sectional view echocardiography is an efficient non-invasive diagnostic tool for characterizing Myocardial Infarction (MI) and stages of expansion leading to heart failure. An automated computer-aided technique of cross-sectional echocardiography feature assessment can aid clinicians in early and more reliable detection of MI patients before subsequent catastrophic post-MI medical conditions. Therefore, this paper proposes a novel Myocardial Infarction Index (MII) to discriminate infarcted and normal myocardium using features extracted from apical cross-sectional views of echocardiograms. The cross-sectional view of normal and MI echocardiography images are represented as textons using Maximum Responses (MR8) filter banks. Fractal Dimension (FD), Higher-Order Statistics (HOS), Hu's moments, Gabor Transform features, Fuzzy Entropy (FEnt), Energy, Local binary Pattern (LBP), Renyi's Entropy (REnt), Shannon's Entropy (ShEnt), and Kapur's Entropy (KEnt) features are extracted from textons. These features are ranked using t-test and fuzzy Max-Relevancy and Min-Redundancy (mRMR) ranking methods. Then, combinations of highly ranked features are used in the formulation and development of an integrated MII. This calculated novel MII is used to accurately and quickly detect infarcted myocardium by using one numerical value. Also, the highly ranked features are subjected to classification using different classifiers for the characterization of normal and MI LV ultrasound images using a minimum number of features. Our current technique is able to characterize MI with an average accuracy of 94.37%, sensitivity of 91.25% and specificity of 97.50% with 8 apical four chambers view features extracted from only single frame per patient making this a more reliable and accurate classification.

Computer-Aided Diagnosis of Depression Using EEG Signals

The complex, nonlinear and non-stationary electroencephalogram (EEG) signals are very tedious to interpret visually and highly difficult to extract the significant features from them. The linear and nonlinear methods are effective in identifying the changes in EEG signals for the detection of depression. Linear methods do not exhibit the complex dynamical variations in the EEG signals. Hence, chaos theory and nonlinear dynamic methods are widely used in extracting the EEG signal features for computer-aided diagnosis (CAD) of depression. Hence, this article presents the recent efforts on CAD of depression using EEG signals with a focus on using nonlinear methods. Such a CAD system is simple to use and may be used by the clinicians as a tool to confirm their diagnosis. It should be of a particular value to enable the early detection of depression.

Empirical mode decomposition analysis of near-infrared spectroscopy muscular signals to assess the effect of physical activity in type 2 diabetic patients

Type 2 diabetes is a metabolic disorder that may cause major problems to several physiological systems. Exercise has proven to be very effective in the prevention, management and improvement of this pathology in patients. Muscle metabolism is often studied with near-infrared spectroscopy (NIRS), a noninvasive technique that can measure changes in the concentration of oxygenated (O2Hb) and reduced hemoglobin (HHb) of tissues. These NIRS signals are highly non-stationary, non-Gaussian and nonlinear in nature. The empirical mode decomposition (EMD) is used as a nonlinear adaptive model to extract information present in the NIRS signals. NIRS signals acquired from the tibialis anterior muscle of controls and type 2 diabetic patients are processed by EMD to yield three intrinsic mode functions (IMF). The sample entropy (SE), fractal dimension (FD), and Hurst exponent (HE) are computed from these IMFs. Subjects are monitored at the beginning of the study and after one year of a physical training programme. Following the exercise programme, we observed an increase in the SE and FD and a decrease in the HE in all diabetic subjects. Our results show the influence of physical exercise program in improving muscle performance and muscle drive by the central nervous system in the patients. A multivariate analysis of variance performed at the end of the training programme also indicated that the NIRS metabolic patterns of controls and diabetic subjects are more similar than at the beginning of the study. Hence, the proposed EMD technique applied to NIRS signals may be very useful to gain a non-invasive understanding of the neuromuscular and vascular impairment in diabetic subjects.

Multiple feature extraction and classification of electroencephalograph signal for Alzheimers' with spectrum and bispectrum

In this paper, we have combined experimental neurophysiologic recording and statistical analysis to investigate the nonlinear characteristic and the cognitive function of the brain. Spectrum and bispectrum analyses are proposed to extract multiple effective features of electroencephalograph (EEG) signals from Alzheimer's disease (AD) patients and further applied to distinguish AD patients from the normal controls. Spectral analysis based on autoregressive Burg method is first used to quantify the power distribution of EEG series in the frequency domain. Compared to the control group, the relative power spectral density of AD group is significantly higher in the theta frequency band, while lower in the alpha frequency bands. In addition, median frequency of spectrum is decreased, and spectral entropy ratio of these two frequency bands undergoes drastic changes at the P3 electrode in the central-parietal brain region, implying that the electrophysiological behavior in AD brain is much slower and less irregular. In order to explore the nonlinear high order information, bispectral analysis which measures the complexity of phase-coupling is further applied to P3 electrode in the whole frequency band. It is demonstrated that less bispectral peaks appear and the amplitudes of peaks fall, suggesting a decrease of non-Gaussianity and nonlinearity of EEG in ADs. Notably, the application of this method to five brain regions shows higher concentration of the weighted center of bispectrum and lower complexity reflecting phase-coupling by bispectral entropy. Based on spectrum and bispectrum analyses, six efficient features are extracted and then applied to discriminate AD from the normal in the five brain regions. The classification results indicate that all these features could differentiate AD patients from the normal controls with a maximum accuracy of 90.2%. Particularly, different brain regions are sensitive to different features. Moreover, the optimal combination of features obtained by discriminant analysis may improve the classification accuracy. These results demonstrate the great promise for scape EEG spectral and bispectral features as a potential effective method for detection of AD, which may facilitate our understanding of the pathological mechanism of the disease.

Inference of human affective states from psychophysiological measurements extracted under ecologically valid conditions

Compared to standard laboratory protocols, the measurement of psychophysiological signals in real world experiments poses technical and methodological challenges due to external factors that cannot be directly controlled. To address this problem, we propose a hybrid approach based on an immersive and human accessible space called the eXperience Induction Machine (XIM), that incorporates the advantages of a laboratory within a life-like setting. The XIM integrates unobtrusive wearable sensors for the acquisition of psychophysiological signals suitable for ambulatory emotion research. In this paper, we present results from two different studies conducted to validate the XIM as a general-purpose sensing infrastructure for the study of human affective states under ecologically valid conditions. In the first investigation, we recorded and classified signals from subjects exposed to pictorial stimuli corresponding to a range of arousal levels, while they were free to walk and gesticulate. In the second study, we designed an experiment that follows the classical conditioning paradigm, a well-known procedure in the behavioral sciences, with the additional feature that participants were free to move in the physical space, as opposed to similar studies measuring physiological signals in constrained laboratory settings. Our results indicate that, by using our sensing infrastructure, it is indeed possible to infer human event-elicited affective states through measurements of psychophysiological signals under ecological conditions.

Revealing real-time emotional responses: a personalized assessment based on heartbeat dynamics

Emotion recognition through computational modeling and analysis of physiological signals has been widely investigated in the last decade. Most of the proposed emotion recognition systems require relatively long-time series of multivariate records and do not provide accurate real-time characterizations using short-time series. To overcome these limitations, we propose a novel personalized probabilistic framework able to characterize the emotional state of a subject through the analysis of heartbeat dynamics exclusively. The study includes thirty subjects presented with a set of standardized images gathered from the international affective picture system, alternating levels of arousal and valence. Due to the intrinsic nonlinearity and nonstationarity of the RR interval series, a specific point-process model was devised for instantaneous identification considering autoregressive nonlinearities up to the third-order according to the Wiener-Volterra representation, thus tracking very fast stimulus-response changes. Features from the instantaneous spectrum and bispectrum, as well as the dominant Lyapunov exponent, were extracted and considered as input features to a support vector machine for classification. Results, estimating emotions each 10 seconds, achieve an overall accuracy in recognizing four emotional states based on the circumplex model of affect of 79.29%, with 79.15% on the valence axis, and 83.55% on the arousal axis.

DEPRESSION DIAGNOSIS SUPPORT SYSTEM BASED ON EEG SIGNAL ENTROPIES

Electroencephalography (EEG) is a measure which represents the functional activity of the brain. We show that a detailed analysis of EEG measurements provides highly discriminant features which indicate the mental state of patients with clinical depression. Our feature extraction method revolves around a novel processing structure that combines wavelet packet decomposition (WPD) and non-linear algorithms. WPD was used to select appropriate EEG frequency bands. The resulting signals were processed with the non-linear measures of approximate entropy (ApEn), sample entropy (SampEn), renyi entropy (REN) and bispectral phase entropy ( P h). The features were selected using t-test and only discriminative features were fed to various classifiers, namely probabilistic neural network (PNN), support vector machine (SVM), decision tree (DT), k-nearest neighbor algorithm (k-NN), naive bayes classification (NBC), Gaussian mixture model (GMM) and Fuzzy Sugeno Classifier (FSC). Our classification results show that, with a classification accuracy of 99.5%, the PNN classifier performed better than the rest of classifiers in discriminating between normal and depression EEG signals. Hence, the proposed decision support system can be used to diagnose, and monitor the treatment of patients suffering from depression.

Electrocardiogram pattern recognition and analysis based on artificial neural networks and support vector machines: a review

Computer systems for Electrocardiogram (ECG) analysis support the clinician in tedious tasks (e.g., Holter ECG monitored in Intensive Care Units) or in prompt detection of dangerous events (e.g., ventricular fibrillation). Together with clinical applications (arrhythmia detection and heart rate variability analysis), ECG is currently being investigated in biometrics (human identification), an emerging area receiving increasing attention. Methodologies for clinical applications can have both differences and similarities with respect to biometrics. This paper reviews methods of ECG processing from a pattern recognition perspective. In particular, we focus on features commonly used for heartbeat classification. Considering the vast literature in the field and the limited space of this review, we dedicated a detailed discussion only to a few classifiers (Artificial Neural Networks and Support Vector Machines) because of their popularity; however, other techniques such as Hidden Markov Models and Kalman Filtering will be also mentioned.

Linear and nonlinear analysis of normal and CAD-affected heart rate signals

Coronary artery disease (CAD) is one of the dangerous cardiac disease, often may lead to sudden cardiac death. It is difficult to diagnose CAD by manual inspection of electrocardiogram (ECG) signals. To automate this detection task, in this study, we extracted the heart rate (HR) from the ECG signals and used them as base signal for further analysis. We then analyzed the HR signals of both normal and CAD subjects using (i) time domain, (ii) frequency domain and (iii) nonlinear techniques. The following are the nonlinear methods that were used in this work: Poincare plots, Recurrence Quantification Analysis (RQA) parameters, Shannon entropy, Approximate Entropy (ApEn), Sample Entropy (SampEn), Higher Order Spectra (HOS) methods, Detrended Fluctuation Analysis (DFA), Empirical Mode Decomposition (EMD), Cumulants, and Correlation Dimension. As a result of the analysis, we present unique recurrence, Poincare and HOS plots for normal and CAD subjects. We have also observed significant variations in the range of these features with respect to normal and CAD classes, and have presented the same in this paper. We found that the RQA parameters were higher for CAD subjects indicating more rhythm. Since the activity of CAD subjects is less, similar signal patterns repeat more frequently compared to the normal subjects. The entropy based parameters, ApEn and SampEn, are lower for CAD subjects indicating lower entropy (less activity due to impairment) for CAD. Almost all HOS parameters showed higher values for the CAD group, indicating the presence of higher frequency content in the CAD signals. Thus, our study provides a deep insight into how such nonlinear features could be exploited to effectively and reliably detect the presence of CAD.

Entropy analysis of muscular near-infrared spectroscopy (NIRS) signals during exercise programme of type 2 diabetic patients: quantitative assessment of muscle metabolic pattern

Diabetes mellitus (DM) is a metabolic disorder that is widely rampant throughout the world population these days. The uncontrolled DM may lead to complications of eye, heart, kidney and nerves. The most common type of diabetes is the type 2 diabetes or insulin-resistant DM. Near-infrared spectroscopy (NIRS) technology is widely used in non-invasive monitoring of physiological signals. Three types of NIRS signals are used in this work: (i) variation in the oxygenated haemoglobin (O2Hb) concentration, (ii) deoxygenated haemoglobin (HHb), and (iii) ratio of oxygenated over the sum of the oxygenated and deoxygenated haemoglobin which is defined as: tissue oxygenation index (TOI) to analyze the effect of exercise on diabetes subjects. The NIRS signal has the characteristics of non-linearity and non-stationarity. Hence, the very small changes in this time series can be efficiently extracted using higher order statistics (HOS) method. Hence, in this work, we have used sample and HOS entropies to analyze these NIRS signals. These computer aided techniques will assist the clinicians to diagnose and monitor the health accurately and easily without any inter or intra observer variability. Results showed that after a one-year of physical exercise programme, all diabetic subjects increased the sample entropy of the NIRS signals, thus revealing a better muscle performance and an improved recruitment by the central nervous system. Moreover, after one year of physical therapy, diabetic subjects showed a NIRS muscular metabolic pattern that was not distinguished from that of controls. We believe that sample and bispectral entropy analysis is need when the aim is to compare the inner structure of the NIRS signals during muscle contraction, particularly when dealing with neuromuscular impairments.

Automated classification of patients with coronary artery disease using grayscale features from left ventricle echocardiographic images

Coronary Artery Disease (CAD), caused by the buildup of plaque on the inside of the coronary arteries, has a high mortality rate. To efficiently detect this condition from echocardiography images, with lesser inter-observer variability and visual interpretation errors, computer based data mining techniques may be exploited. We have developed and presented one such technique in this paper for the classification of normal and CAD affected cases. A multitude of grayscale features (fractal dimension, entropies based on the higher order spectra, features based on image texture and local binary patterns, and wavelet based features) were extracted from echocardiography images belonging to a huge database of 400 normal cases and 400 CAD patients. Only the features that had good discriminating capability were selected using t-test. Several combinations of the resultant significant features were used to evaluate many supervised classifiers to find the combination that presents a good accuracy. We observed that the Gaussian Mixture Model (GMM) classifier trained with a feature subset made up of nine significant features presented the highest accuracy, sensitivity, specificity, and positive predictive value of 100%. We have also developed a novel, highly discriminative HeartIndex, which is a single number that is calculated from the combination of the features, in order to objectively classify the images from either of the two classes. Such an index allows for an easier implementation of the technique for automated CAD detection in the computers in hospitals and clinics.

Automated diagnosis of normal and alcoholic EEG signals

Electroencephalogram (EEG) signals, which record the electrical activity in the brain, are useful for assessing the mental state of a person. Since these signals are nonlinear and non-stationary in nature, it is very difficult to decipher the useful information from them using conventional statistical and frequency domain methods. Hence, the application of nonlinear time series analysis to EEG signals could be useful to study the dynamical nature and variability of the brain signals. In this paper, we propose a Computer Aided Diagnostic (CAD) technique for the automated identification of normal and alcoholic EEG signals using nonlinear features. We first extract nonlinear features such as Approximate Entropy (ApEn), Largest Lyapunov Exponent (LLE), Sample Entropy (SampEn), and four other Higher Order Spectra (HOS) features, and then use them to train Support Vector Machine (SVM) classifier of varying kernel functions: 1st, 2nd, and 3rd order polynomials and a Radial basis function (RBF) kernel. Our results indicate that these nonlinear measures are good discriminators of normal and alcoholic EEG signals. The SVM classifier with a polynomial kernel of order 1 could distinguish the two classes with an accuracy of 91.7%, sensitivity of 90% and specificity of 93.3%. As a pre-analysis step, the EEG signals were tested for nonlinearity using surrogate data analysis and we found that there was a significant difference in the LLE measure of the actual data and the surrogate data.

NONLINEAR INDICES OF HEART RATE VARIABILITY FOR DIFFERENTIATING ARRHYTHMIAS

Heart rate variability (HRV) is the temporal variation between sequences of consecutive heartbeats. Chaos and fractal-based measurements have been widely used for quantifying the HRV for cardiac risk stratification purposes. In this paper, five different sets of HRVs, viz., normal sinus rhythm (NSR), congestive heart failure (CHF), cardiac arrhythmia suppression trial (CAST), supra ventricular tachyarrhythmia (SVTA) and atrial fibrillation (AF), have been analysed using nonlinear parameters to fix the ranges of each parameter. Data were downloaded from the PhysioNet database with 15 sets in each case. The parameters used for analysis were Poincare plot measures: SD1, SD2 and SD12, largest Lyapunov exponent (LLE), correlation dimension (CD); recurrence plot measures: recurrence rate (REC), determinism (DET), mean diagonal length (L mean ), maximal diagonal length (L max ) and entropy (ENTR); detrended fluctuation analysis measures: scaling exponent (α) and fractal dimension (FD); sample entropy (SampEn); and approximate entropy (ApEn). Analysis of variance (ANOVA) was done for confirming the differences in parameter values between various cases. All parameters except LLE showed a significant statistical difference for different cases.

APPLICATION OF HIGHER ORDER CUMULANT FEATURES FOR CARDIAC HEALTH DIAGNOSIS USING ECG SIGNALS

Electrocardiogram (ECG) is the electrical activity of the heart indicated by P, Q-R-S and T wave. The minute changes in the amplitude and duration of ECG depicts a particular type of cardiac abnormality. It is very difficult to decipher the hidden information present in this nonlinear and nonstationary signal. An automatic diagnostic system that characterizes cardiac activities in ECG signals would provide more insight into these phenomena thereby revealing important clinical information. Various methods have been proposed to detect cardiac abnormalities in ECG recordings. Application of higher order spectra (HOS) features is a seemingly promising approach because it can capture the nonlinear and dynamic nature of the ECG signals. In this paper, we have automatically classified five types of beats using HOS features (higher order cumulants) using two different approaches. The five types of ECG beats are normal (N), right bundle branch block (RBBB), left bundle branch block (LBBB), atrial premature contraction (APC) and ventricular premature contraction (VPC). In the first approach, cumulant features of segmented ECG signal were used for classification; whereas in the second approach cumulants of discrete wavelet transform (DWT) coefficients were used as features for classifiers. In both approaches, the cumulant features were subjected to data reduction using principal component analysis (PCA) and classified using three layer feed-forward neural network (NN) and least square — support vector machine (LS-SVM) classifiers. In this study, we obtained the highest average accuracy of 94.52%, sensitivity of 98.61% and specificity of 98.41% using first approach with NN classifier. The developed system is ready clinically to run on large datasets.

Prostate tissue characterization/classification in 144 patient population using wavelet and higher order spectra features from transrectal ultrasound images

In this work, we have proposed an on-line computer-aided diagnostic system called "UroImage" that classifies a Transrectal Ultrasound (TRUS) image into cancerous or non-cancerous with the help of non-linear Higher Order Spectra (HOS) features and Discrete Wavelet Transform (DWT) coefficients. The UroImage system consists of an on-line system where five significant features (one DWT-based feature and four HOS-based features) are extracted from the test image. These on-line features are transformed by the classifier parameters obtained using the training dataset to determine the class. We trained and tested six classifiers. The dataset used for evaluation had 144 TRUS images which were split into training and testing sets. Three-fold and ten-fold cross-validation protocols were adopted for training and estimating the accuracy of the classifiers. The ground truth used for training was obtained using the biopsy results. Among the six classifiers, using 10-fold cross-validation technique, Support Vector Machine and Fuzzy Sugeno classifiers presented the best classification accuracy of 97.9% with equally high values for sensitivity, specificity and positive predictive value. Our proposed automated system, which achieved more than 95% values for all the performance measures, can be an adjunct tool to provide an initial diagnosis for the identification of patients with prostate cancer. The technique, however, is limited by the limitations of 2D ultrasound guided biopsy, and we intend to improve our technique by using 3D TRUS images in the future.

AUTOMATED DIAGNOSIS OF EPILEPSY USING CWT, HOS AND TEXTURE PARAMETERS

Epilepsy is a chronic brain disorder which manifests as recurrent seizures. Electroencephalogram (EEG) signals are generally analyzed to study the characteristics of epileptic seizures. In this work, we propose a method for the automated classification of EEG signals into normal, interictal and ictal classes using Continuous Wavelet Transform (CWT), Higher Order Spectra (HOS) and textures. First the CWT plot was obtained for the EEG signals and then the HOS and texture features were extracted from these plots. Then the statistically significant features were fed to four classifiers namely Decision Tree (DT), K-Nearest Neighbor (KNN), Probabilistic Neural Network (PNN) and Support Vector Machine (SVM) to select the best classifier. We observed that the SVM classifier with Radial Basis Function (RBF) kernel function yielded the best results with an average accuracy of 96%, average sensitivity of 96.9% and average specificity of 97% for 23.6 s duration of EEG data. Our proposed technique can be used as an automatic seizure monitoring software. It can also assist the doctors to cross check the efficacy of their prescribed drugs.

DETECTION OF HUMAN STRESS USING SHORT-TERM ECG AND HRV SIGNALS

This paper introduces a method for resolving the problem of human stress detection through short-term (less than 5 min) electrocardiogram (ECG) and heart rate variability (HRV) signals. The explored methodology helps to improve the stress detection rate and reliability through multiple evidences originated in same sensor. In this work, stress-inducing protocol, data acquisition, preprocessing, feature extraction and classification are the major steps involved to detect the stress. In total, 60 subjects (30 males and 30 females) participated in the Stroop color word-based stress-inducing task and ECG signal was acquired simultaneously. The wavelet denoising algorithm was applied to remove high frequency, baseline wander and power line noises. Discrete wavelet transform (DWT)-based heart rate (HR) detection algorithm is used for deriving HRV signal from the preprocessed ECG signal. The ectopic beat removal method is employed to eliminate the ectopic beat and noise peaks in the HRV signal. In order to detect the stress, the issue of uneven sampling with the HRV signal has been successfully rectified using the Lomb-Scargle periodogram (LSP). The application of LSP in short-term HRV signals (32 s), uneven sampling issue, and power spectral information issue has been rectified and the trustworthiness of the short-term HRV signal has been proved by hypothesis as well as experimental results. Theoretical analysis suggested that a minimum 25 s of online or offline ECG data is required to analyze the autonomous nervous system (ANS) activity related to stress. In addition to the HRV signal, ECG-based stress assessment has been proposed to detect the stress through optimum features using fast Fourier transform (FFT). Various features extracted from the ECG and HRV signal have been classified into normal and stress using PNN and kNN classifiers with different smoothing factor and k values. The experimental results indicate that the proposed methodology for short-term ECG and HRV signal can achieve the overall average classification accuracy of 91.66% and 94.66% in the subject-independent mode.