A Markov-Switching Model Approach to Heart Sound Segmentation and Classification.

OBJECTIVE: We consider challenges in accurate segmentation of heart sound signals recorded under noisy clinical environments for subsequent classification of pathological events. Existing state-of-the-art solutions to heart sound segmentation use probabilistic models such as hidden Markov models (HMMs), which, however, are limited by its observation independence assumption and rely on pre-extraction of noise-robust features. METHODS: We propose a Markov-switching autoregressive (MSAR) process to model the raw heart sound signals directly, which allows efficient segmentation of the cyclical heart sound states according to the distinct dependence structure in each state. To enhance robustness, we extend the MSAR model to a switching linear dynamic system (SLDS) that jointly model both the switching AR dynamics of underlying heart sound signals and the noise effects. We introduce a novel algorithm via fusion of switching Kalman filter and the duration-dependent Viterbi algorithm, which incorporates the duration of heart sound states to improve state decoding. RESULTS: Evaluated on Physionet/CinC Challenge 2016 dataset, the proposed MSAR-SLDS approach significantly outperforms the hidden semi-Markov model (HSMM) in heart sound segmentation based on raw signals and comparable to a feature-based HSMM. The segmented labels were then used to train Gaussian-mixture HMM classifier for identification of abnormal beats, achieving high average precision of 86.1% on the same dataset including very noisy recordings. CONCLUSION: The proposed approach shows noticeable performance in heart sound segmentation and classification on a large noisy dataset. SIGNIFICANCE: It is potentially useful in developing automated heart monitoring systems for pre-screening of heart pathologies.

Estimating Dynamic Connectivity States in fMRI Using Regime-Switching Factor Models.

We consider the challenges in estimating the state-related changes in brain connectivity networks with a large number of nodes. Existing studies use the sliding-window analysis or time-varying coefficient models, which are unable to capture both smooth and abrupt changes simultaneously, and rely on ad-hoc approaches to the high-dimensional estimation. To overcome these limitations, we propose a Markov-switching dynamic factor model, which allows the dynamic connectivity states in functional magnetic resonance imaging (fMRI) data to be driven by lower-dimensional latent factors. We specify a regime-switching vector autoregressive (SVAR) factor process to quantity the time-varying directed connectivity. The model enables a reliable, data-adaptive estimation of change-points of connectivity regimes and the massive dependencies associated with each regime. We develop a three-step estimation procedure: 1) extracting the factors using principal component analysis, 2) identifying connectivity regimes in a low-dimensional subspace based on the factor-based SVAR model, and 3) constructing high-dimensional state connectivity metrics based on the subspace estimates. Simulation results show that our estimator outperforms -means clustering of time-windowed coefficients, providing more accurate estimate of time-evolving connectivity. It achieves percentage of reduction in mean squared error by 60% when the network dimension is comparable to the sample size. When applied to the resting-state fMRI data, our method successfully identifies modular organization in the resting-statenetworksin consistencywith other studies. It further reveals changes in brain states with variations across subjects and distinct large-scale directed connectivity patterns across states.

A Unified Estimation Framework for State-Related Changes in Effective Brain Connectivity.

OBJECTIVE: This paper addresses the critical problem of estimating time-evolving effective brain connectivity. Current approaches based on sliding window analysis or time-varying coefficient models do not simultaneously capture both slow and abrupt changes in causal interactions between different brain regions. METHODS: To overcome these limitations, we develop a unified framework based on a switching vector autoregressive (SVAR) model. Here, the dynamic connectivity regimes are uniquely characterized by distinct vector autoregressive (VAR) processes and allowed to switch between quasi-stationary brain states. The state evolution and the associated directed dependencies are defined by a Markov process and the SVAR parameters. We develop a three-stage estimation algorithm for the SVAR model: 1) feature extraction using time-varying VAR (TV-VAR) coefficients, 2) preliminary regime identification via clustering of the TV-VAR coefficients, 3) refined regime segmentation by Kalman smoothing and parameter estimation via expectation-maximization algorithm under a state-space formulation, using initial estimates from the previous two stages. RESULTS: The proposed framework is adaptive to state-related changes and gives reliable estimates of effective connectivity. Simulation results show that our method provides accurate regime change-point detection and connectivity estimates. In real applications to brain signals, the approach was able to capture directed connectivity state changes in functional magnetic resonance imaging data linked with changes in stimulus conditions, and in epileptic electroencephalograms, differentiating ictal from nonictal periods. CONCLUSION: The proposed framework accurately identifies state-dependent changes in brain network and provides estimates of connectivity strength and directionality. SIGNIFICANCE: The proposed approach is useful in neuroscience studies that investigate the dynamics of underlying brain states.

Electroencephalographic Motor Imagery Brain Connectivity Analysis for BCI: A Review.

Recent research has reached a consensus on the feasibility of motor imagery brain-computer interface (MI-BCI) for different applications, especially in stroke rehabilitation. Most MI-BCI systems rely on temporal, spectral, and spatial features of single channels to distinguish different MI patterns. However, no successful communication has been established for a completely locked-in subject. To provide more useful and informative features, it has been recommended to take into account the relationships among electroencephalographic (EEG) sensor/source signals in the form of brain connectivity as an efficient tool of neuroscience. In this review, we briefly report the challenges and limitations of conventional MI-BCIs. Brain connectivity analysis, particularly functional and effective, has been described as one of the most promising approaches for improving MI-BCI performance. An extensive literature on EEG-based MI brain connectivity analysis of healthy subjects is reviewed. We subsequently discuss the brain connectomes during left and right hand, feet, and tongue MI movements. Moreover, key components involved in brain connectivity analysis that considerably affect the results are explained. Finally, possible technical shortcomings that may have influenced the results in previous research are addressed and suggestions are provided.

Is First-Order Vector Autoregressive Model Optimal for fMRI Data?

We consider the problem of selecting the optimal orders of vector autoregressive (VAR) models for fMRI data. Many previous studies used model order of one and ignored that it may vary considerably across data sets depending on different data dimensions, subjects, tasks, and experimental designs. In addition, the classical information criteria (IC) used (e.g., the Akaike IC (AIC)) are biased and inappropriate for the high-dimensional fMRI data typically with a small sample size. We examine the mixed results on the optimal VAR orders for fMRI, especially the validity of the order-one hypothesis, by a comprehensive evaluation using different model selection criteria over three typical data types--a resting state, an event-related design, and a block design data set--with varying time series dimensions obtained from distinct functional brain networks. We use a more balanced criterion, Kullback's IC (KIC) based on Kullback's symmetric divergence combining two directed divergences. We also consider the bias-corrected versions (AICc and KICc) to improve VAR model selection in small samples. Simulation results show better small-sample selection performance of the proposed criteria over the classical ones. Both bias-corrected ICs provide more accurate and consistent model order choices than their biased counterparts, which suffer from overfitting, with KICc performing the best. Results on real data show that orders greater than one were selected by all criteria across all data sets for the small to moderate dimensions, particularly from small, specific networks such as the resting-state default mode network and the task-related motor networks, whereas low orders close to one but not necessarily one were chosen for the large dimensions of full-brain networks.

Morphological study of the newly designed cementless femoral stem.

A morphology study was essential to the development of the cementless femoral stem because accurate dimensions for both the periosteal and endosteal canal ensure primary fixation stability for the stem, bone interface, and prevent stress shielding at the calcar region. This paper focused on a three-dimensional femoral model for Asian patients that applied preoperative planning and femoral stem design. We measured various femoral parameters such as the femoral head offset, collodiaphyseal angle, bowing angle, anteversion, and medullary canal diameters from the osteotomy level to 150 mm below the osteotomy level to determine the position of the isthmus. Other indices and ratios for the endosteal canal, metaphyseal, and flares were computed and examined. The results showed that Asian femurs are smaller than Western femurs, except in the metaphyseal region. The canal flare index (CFI) was poorly correlated (r < 0.50) to the metaphyseal canal flare index (MCFI), but correlated well (r = 0.66) with the corticomedullary index (CMI). The diversity of the femoral size, particularly in the metaphyseal region, allows for proper femoral stem design for Asian patients, improves osseointegration, and prolongs the life of the implant.

Primary stability recognition of the newly designed cementless femoral stem using digital signal processing.

Stress shielding and micromotion are two major issues which determine the success of newly designed cementless femoral stems. The correlation of experimental validation with finite element analysis (FEA) is commonly used to evaluate the stress distribution and fixation stability of the stem within the femoral canal. This paper focused on the applications of feature extraction and pattern recognition using support vector machine (SVM) to determine the primary stability of the implant. We measured strain with triaxial rosette at the metaphyseal region and micromotion with linear variable direct transducer proximally and distally using composite femora. The root mean squares technique is used to feed the classifier which provides maximum likelihood estimation of amplitude, and radial basis function is used as the kernel parameter which mapped the datasets into separable hyperplanes. The results showed 100% pattern recognition accuracy using SVM for both strain and micromotion. This indicates that DSP could be applied in determining the femoral stem primary stability with high pattern recognition accuracy in biomechanical testing.

Design process of cementless femoral stem using a nonlinear three dimensional finite element analysis.

BACKGROUND: Minimal available information concerning hip morphology is the motivation for several researchers to study the difference between Asian and Western populations. Current use of a universal hip stem of variable size is not the best option for all femur types. This present study proposed a new design process of the cementless femoral stem using a three dimensional model which provided more information and accurate analysis compared to conventional methods. METHODS: This complete design cycle began with morphological analysis, followed by femoral stem design, fit and fill analysis, and nonlinear finite element analysis (FEA). Various femur parameters for periosteal and endosteal canal diameters are measured from the osteotomy level to 150 mm below to determine the isthmus position. RESULTS: The results showed better total fit (53.7%) and fill (76.7%) canal, with more load distributed proximally to prevent stress shielding at calcar region. The stem demonstrated lower displacement and micromotion (less than 40 µm) promoting osseointegration between the stem-bone and providing primary fixation stability. CONCLUSION: This new design process could be used as a preclinical assessment tool and will shorten the design cycle by identifying the major steps which must be taken while designing the femoral stem.

Fabrication of low-cost, cementless femoral stem 316L stainless steel using investment casting technique.

Total hip arthroplasty is a flourishing orthopedic surgery, generating billions of dollars of revenue. The cost associated with the fabrication of implants has been increasing year by year, and this phenomenon has burdened the patient with extra charges. Consequently, this study will focus on designing an accurate implant via implementing the reverse engineering of three-dimensional morphological study based on a particular population. By using finite element analysis, this study will assist to predict the outcome and could become a useful tool for preclinical testing of newly designed implants. A prototype is then fabricated using 316L stainless steel by applying investment casting techniques that reduce manufacturing cost without jeopardizing implant quality. The finite element analysis showed that the maximum von Mises stress was 66.88 MPa proximally with a safety factor of 2.39 against endosteal fracture, and micromotion was 4.73 µm, which promotes osseointegration. This method offers a fabrication process of cementless femoral stems with lower cost, subsequently helping patients, particularly those from nondeveloped countries.

Discriminative tandem features for HMM-based EEG classification.

We investigate the use of discriminative feature extractors in tandem configuration with generative EEG classification system. Existing studies on dynamic EEG classification typically use hidden Markov models (HMMs) which lack discriminative capability. In this paper, a linear and a non-linear classifier are discriminatively trained to produce complementary input features to the conventional HMM system. Two sets of tandem features are derived from linear discriminant analysis (LDA) projection output and multilayer perceptron (MLP) class-posterior probability, before appended to the standard autoregressive (AR) features. Evaluation on a two-class motor-imagery classification task shows that both the proposed tandem features yield consistent gains over the AR baseline, resulting in significant relative improvement of 6.2% and 11.2% for the LDA and MLP features respectively. We also explore portability of these features across different subjects.

Linear dynamic models for classification of single-trial EEG.

This paper investigates the use of linear dynamic models (LDMs) to improve classification of single-trial EEG signals. Existing dynamic classification of EEG uses discrete-state hidden Markov models (HMMs) based on piecewise-stationary assumption, which is inadequate for modeling the highly non-stationary dynamics underlying EEG. The continuous hidden states of LDMs could better describe this continuously changing characteristic of EEG, and thus improve the classification performance. We consider two examples of LDM: a simple local level model (LLM) and a time-varying autoregressive (TVAR) state-space model. AR parameters and band power are used as features. Parameter estimation of the LDMs is performed by using expectation-maximization (EM) algorithm. We also investigate different covariance modeling of Gaussian noises in LDMs for EEG classification. The experimental results on two-class motor-imagery classification show that both types of LDMs outperform the HMM baseline, with the best relative accuracy improvement of 14.8% by LLM with full covariance for Gaussian noises. It may due to that LDMs offer more flexibility in fitting the underlying dynamics of EEG.

EMG-based facial gesture recognition through versatile elliptic basis function neural network.

BACKGROUND: Recently, the recognition of different facial gestures using facial neuromuscular activities has been proposed for human machine interfacing applications. Facial electromyograms (EMGs) analysis is a complicated field in biomedical signal processing where accuracy and low computational cost are significant concerns. In this paper, a very fast versatile elliptic basis function neural network (VEBFNN) was proposed to classify different facial gestures. The effectiveness of different facial EMG time-domain features was also explored to introduce the most discriminating. METHODS: In this study, EMGs of ten facial gestures were recorded from ten subjects using three pairs of surface electrodes in a bi-polar configuration. The signals were filtered and segmented into distinct portions prior to feature extraction. Ten different time-domain features, namely, Integrated EMG, Mean Absolute Value, Mean Absolute Value Slope, Maximum Peak Value, Root Mean Square, Simple Square Integral, Variance, Mean Value, Wave Length, and Sign Slope Changes were extracted from the EMGs. The statistical relationships between these features were investigated by Mutual Information measure. Then, the feature combinations including two to ten single features were formed based on the feature rankings appointed by Minimum-Redundancy-Maximum-Relevance (MRMR) and Recognition Accuracy (RA) criteria. In the last step, VEBFNN was employed to classify the facial gestures. The effectiveness of single features as well as the feature sets on the system performance was examined by considering the two major metrics, recognition accuracy and training time. Finally, the proposed classifier was assessed and compared with conventional methods support vector machines and multilayer perceptron neural network. RESULTS: The average classification results showed that the best performance for recognizing facial gestures among all single/multi-features was achieved by Maximum Peak Value with 87.1% accuracy. Moreover, the results proved a very fast procedure since the training time during classification via VEBFNN was 0.105 seconds. It was also indicated that MRMR was not a proper criterion to be used for making more effective feature sets in comparison with RA. CONCLUSIONS: This work was accomplished by introducing the most discriminating facial EMG time-domain feature for the recognition of different facial gestures; and suggesting VEBFNN as a promising method in EMG-based facial gesture classification to be used for designing interfaces in human machine interaction systems.

An expectation-maximization algorithm based Kalman smoother approach for single-trial estimation of event-related potentials.

This paper applies an expectation-maximization (EM) based Kalman smoother (KS) approach for single-trial event-related potential (ERP) estimation. Existing studies assume a Markov diffusion process for the dynamics of ERP parameters which is recursively estimated by optimal filtering approaches such as Kalman filter (KF). However, these studies only consider estimation of ERP state parameters while the model parameters are pre-specified using manual tuning, which is time-consuming for practical usage besides giving suboptimal estimates. We extend the KF approach by adding EM based maximum likelihood estimation of the model parameters to obtain more accurate ERP estimates automatically. We also introduce different model variants by allowing flexibility in the covariance structure of model noises. Optimal model selection is performed based on Akaike Information Criterion (AIC). The method is applied to estimation of chirp-evoked auditory brainstem responses (ABRs) for detection of wave V critical for assessment of hearing loss. Results shows that use of more complex covariances are better estimating inter-trial variability.

An artifacts removal post-processing for epiphyseal region-of-interest (EROI) localization in automated bone age assessment (BAA).

BACKGROUND: Segmentation is the most crucial part in the computer-aided bone age assessment. A well-known type of segmentation performed in the system is adaptive segmentation. While providing better result than global thresholding method, the adaptive segmentation produces a lot of unwanted noise that could affect the latter process of epiphysis extraction. METHODS: A proposed method with anisotropic diffusion as pre-processing and a novel Bounded Area Elimination (BAE) post-processing algorithm to improve the algorithm of ossification site localization technique are designed with the intent of improving the adaptive segmentation result and the region-of interest (ROI) localization accuracy. RESULTS: The results are then evaluated by quantitative analysis and qualitative analysis using texture feature evaluation. The result indicates that the image homogeneity after anisotropic diffusion has improved averagely on each age group for 17.59%. Results of experiments showed that the smoothness has been improved averagely 35% after BAE algorithm and the improvement of ROI localization has improved for averagely 8.19%. The MSSIM has improved averagely 10.49% after performing the BAE algorithm on the adaptive segmented hand radiograph. CONCLUSIONS: The result indicated that hand radiographs which have undergone anisotropic diffusion have greatly reduced the noise in the segmented image and the result as well indicated that the BAE algorithm proposed is capable of removing the artifacts generated in adaptive segmentation.

Spectral estimation of nonstationary EEG using particle filtering with application to event-related desynchronization (ERD).

This paper proposes non-Gaussian models for parametric spectral estimation with application to event-related desynchronization (ERD) estimation of nonstationary EEG. Existing approaches for time-varying spectral estimation use time-varying autoregressive (TVAR) state-space models with Gaussian state noise. The parameter estimation is solved by a conventional Kalman filtering. This study uses non-Gaussian state noise to model autoregressive (AR) parameter variation with estimation by a Monte Carlo particle filter (PF). Use of non-Gaussian noise such as heavy-tailed distribution is motivated by its ability to track abrupt and smooth AR parameter changes, which are inadequately modeled by Gaussian models. Thus, more accurate spectral estimates and better ERD tracking can be obtained. This study further proposes a non-Gaussian state space formulation of time-varying autoregressive moving average (TVARMA) models to improve the spectral estimation. Simulation on TVAR process with abrupt parameter variation shows superior tracking performance of non-Gaussian models. Evaluation on motor-imagery EEG data shows that the non-Gaussian models provide more accurate detection of abrupt changes in alpha rhythm ERD. Among the proposed non-Gaussian models, TVARMA shows better spectral representations while maintaining reasonable good ERD tracking performance.

Human facial neural activities and gesture recognition for machine-interfacing applications.

The authors present a new method of recognizing different human facial gestures through their neural activities and muscle movements, which can be used in machine-interfacing applications. Human-machine interface (HMI) technology utilizes human neural activities as input controllers for the machine. Recently, much work has been done on the specific application of facial electromyography (EMG)-based HMI, which have used limited and fixed numbers of facial gestures. In this work, a multipurpose interface is suggested that can support 2-11 control commands that can be applied to various HMI systems. The significance of this work is finding the most accurate facial gestures for any application with a maximum of eleven control commands. Eleven facial gesture EMGs are recorded from ten volunteers. Detected EMGs are passed through a band-pass filter and root mean square features are extracted. Various combinations of gestures with a different number of gestures in each group are made from the existing facial gestures. Finally, all combinations are trained and classified by a Fuzzy c-means classifier. In conclusion, combinations with the highest recognition accuracy in each group are chosen. An average accuracy >90% of chosen combinations proved their ability to be used as command controllers.