The human nervous system is an ensemble of connected neuronal networks. Modeling and system identification of the human nervous system helps us understand how the brain processes sensory input and controls responses at the systems level. This study aims to propose an advanced approach based on a hierarchical neural network and non-linear system identification method to model neural activity in the nervous system in response to an external somatosensory input. The proposed approach incorporates basic concepts of Non-linear AutoRegressive Moving Average Model with eXogenous input (NARMAX) and neural network to acknowledge non-linear closed-loop neural interactions. Different from the commonly used polynomial NARMAX method, the proposed approach replaced the polynomial non-linear terms with a hierarchical neural network. The hierarchical neural network is built based on known neuroanatomical connections and corresponding transmission delays in neural pathways. The proposed method is applied to an experimental dataset, where cortical activities from ten young able-bodied individuals are extracted from electroencephalographic signals while applying mechanical perturbations to their wrist joint. The results yielded by the proposed method were compared with those obtained by the polynomial NARMAX and Volterra methods, evaluated by the variance accounted for (VAF). Both the proposed and polynomial NARMAX methods yielded much better modeling results than the Volterra model. Furthermore, the proposed method modeled cortical responded with a mean VAF of 69.35% for a three-step ahead prediction, which is significantly better than the VAF from a polynomial NARMAX model (mean VAF 47.09%). This study provides a novel approach for precise modeling of cortical responses to sensory input. The results indicate that the incorporation of knowledge of neuroanatomical connections in building a realistic model greatly improves the performance of system identification of the human nervous system.

In hemiparetic stroke, functional recovery of paretic limb may occur with the reorganization of neural networks in the brain. Neuroimaging techniques, such as magnetic resonance imaging (MRI), have a high spatial resolution which can be used to reveal anatomical changes in the brain following a stroke. However, low temporal resolution of MRI provides less insight of dynamic changes of brain activity. In contrast, electro-neurophysiological techniques, such as electroencephalography (EEG), have an excellent temporal resolution to measure such transient events, however are hindered by its low spatial resolution. This proof-of-principle study assessed a novel multimodal brain imaging technique namely Variational Bayesian Multimodal Encephalography (VBMEG), which aims to improve the spatial resolution of EEG for tracking the information flow inside the brain and its changes following a stroke. The limitations of EEG are complemented by constraints derived from anatomical MRI and diffusion weighted imaging (DWI). EEG data were acquired from individuals suffering from a stroke as well as able-bodied participants while electrical stimuli were delivered sequentially at their index finger in the left and right hand, respectively. The locations of active sources related to this stimulus were precisely identified, resulting in high Variance Accounted For (VAF above 80%). An accurate estimation of dynamic information flow between sources was achieved in this study, showing a high VAF (above 90%) in the cross-validation test. The estimated dynamic information flow was compared between chronic hemiparetic stroke and able-bodied individuals. The results demonstrate the feasibility of VBMEG method in revealing the changes of information flow in the brain after stroke. This study verified the VBMEG method as an advanced computational approach to track the dynamic information flow in the brain following a stroke. This may lead to the development of a quantitative tool for monitoring functional changes of the cortical neural networks after a unilateral brain injury and therefore facilitate the research into, and the practice of stroke rehabilitation.

The essential task of a motor imagery brain–computer interface (BCI) is to extract the motor imagery-related features from electroencephalogram (EEG) signals for classifying motor intentions. However, the optimal frequency band and time segment for extracting such features differ from subject to subject. In this work, we aim to improve the multi-class classification and to reduce the required EEG channel in motor imagery-based BCI by subject-specific time-frequency selection. Our method is based on a criterion namely Fisher discriminant analysis-type F-score to simultaneously select the optimal frequency band and time segment for multi-class classification. The proposed method uses only few Laplacian EEG channels (C3, Cz and C4) located around the sensorimotor area for classification. Applied to a standard multi-class BCI dataset (BCI competition III dataset IIIa), our method leads to better classification performance and smaller standard deviation across subjects compared to the state-of-art methods. Moreover, adding artifacts contaminated trials to the training dataset does not necessarily deteriorate our classification results, indicating that our method is tolerant to artifacts.

Communication between neuronal populations is facilitated by synchronization of their oscillatory activity. Although nonlinearity has been observed in the sensorimotor system, its nonlinear connectivity has not been widely investigated yet. This study investigates nonlinear connectivity during the human stretch reflex based on neuronal synchronization. Healthy participants generated isotonic wrist flexion while receiving a periodic mechanical perturbation to the wrist. Using a novel cross-frequency phase coupling metric, we estimate directional nonlinear connectivity, including time delay, from the perturbation to brain and to muscle, as well as from brain to muscle. Nonlinear phase coupling is significantly stronger from the perturbation to the muscle than to the brain, with a shorter time delay. The time delay from the perturbation to the muscle is 33 ms, similar to the reported latency of the spinal stretch reflex at the wrist. Source localization of nonlinear phase coupling from the brain to the muscle suggests activity originating from the motor cortex, although its effect on the stretch reflex is weak. As such nonlinear phase coupling between the perturbation and muscle activity is dominated by the spinal reflex loop. This study provides new evidence of nonlinear neuronal synchronization in the stretch reflex at the wrist joint with respect to spinal and transcortical loops.

Neural systems can present various types of nonlinear input-output relationships, such as harmonic, subharmonic, and/or intermodulation coupling. This paper aims to introduce a general framework in frequency domain for detecting and characterizing nonlinear coupling in neural systems, called the cross-frequency coherence framework (CFCF). CFCF is an extension of classic coherence based on higher-order statistics. We demonstrate an application of CFCF for identifying nonlinear interactions in human motion control. Our results indicate that CFCF can effectively characterize nonlinear properties of the afferent sensory pathway. We conclude that CFCF contributes to identifying nonlinear transfer in neural systems.