Functional data analysis for BCI and biomedical signals

Vadim, 2023

Brain-computer interfaces require a sophisticated forecasting model. This model fits heterogeneous data. The signals come from ECoG, ECG, fMRI, hand and eye movements, and audio-video sources. The model must reconstruct hidden dependencies in these signals and establish relations between brain signals and limb motions. My research focuses on the construction of BCI forecasting models. I use deep learning models. The main challenges of the study are phase space construction, dimensionality reduction, manifold learning, heterogeneous modeling, and knowledge transfer. Since the measured data are stochastic and contain errors, I actively use and develop Bayesian model selection methods. These methods infer criteria to optimize model structure and parameters. They aim to select an accurate and robust BCI model.

Brain signals and dimensionality reduction

Intracranial electroencephalography [1] signals are tensors or time-related tensor fields. They have several indexes for physical space, time, and frequency. The multi-index structure of time series causes redundancy of space features and multi-correlation [2]. It turns out to increase the complexity of the model and obtain unstable forecasts [3]. I address the dimensionality reduction problem for high-dimensional data. The essential methods are tensor and high-order singular value decomposition [4]. We proposed a feature selection method to reveal hidden dependencies in data [5]. It minimizes multi-correlation in the source space and maximizes the relation between source and target spaces. This solution shrinks the number of model parameters tenfold and stabilizes the forecast. To boost the accuracy, I plan to investigate dimensionality reduction. I consider deep learning models in discrete time: the stacks of autoencoders and recurrent neural networks, and in continuous time: the neural ordinal differential equations [6].

Biomedical signal decoding and multi-modeling

The BCI models are the signal decoding models. It is a special class of models that includes canonical correlation analysis for multivariate and tensor variables. I plan to study the problem of model selection to restore hidden dependencies in the source and target spaces. For example, periodical limb movements cause multiple-correlations in the target space. We proposed to reduce the dimension by projecting the source and target in the latent space. Linear and non-linear methods match predictive models in spaces of high dimensions. Recently we proposed a feature selection algorithm for linear models and tested it on ECoG signals. I plan to develop this algorithm for tensor dimensionality reduction. The base method to compare is High-order Partial Least Squares. An exemplary problem for development is manifold learning. The manifold is solution to the neural partial differential equations. The challange is the find an optimal dimensionality of the manifold.

Heterogeneous data and multi-modeling

The new studies of brain activity fruitfully deliver a variety of measurements. For a group of patients, they contain audio, video, iEEEG-ECoG, ECG, fMRI, and hand or eye movements. These data sets require multi models. Each patient has its peculiarities. It requires a method to transfer knowledge from one patient's model to another. Knowledge transfer for heterogenous models is an important part of my investigation. I use Bayesian inference for multimodel selection to construct an ensemble of models and teacher-student pairs. The information, gained by the properly trained models serves as a prior distribution for a student model.

Continous-time physical activity recognition

A forecast of limb motions stands on the precedents. These precedents, quasi-periodic time series, form a phase trajectory. It is a basic cycle of motion. This trajectory is a loop whose parameters define a class of movement. A sequence of these classes forms the physical human behavior pattern. Recently we proposed human activity recognition algorithm based on the data from wearable sensors. The solution is based on the hierarchical representation of activities as sets of low-level actions. The hierarchical representation provides an interpretable description of studied activities in terms of actions.

Functional data analysis

The brain functional mapping methods verify the signal diffusion hypothesis. It tells that changes in cortical activity zones over the intracranial space control limb movements. The model must consider the spatial structure of the signals. Due to the lack of a common definition of the neighborhood on the spherical surface of the brain, convolutional neural networks cannot be effectively applied to account for spatial information. We proposed a graph representation of the signal. It reveals interrelationships of different areas of intracranial activity and provides a neurobiological interpretation of the functional connections. I plan to develop various methods for constructing a connectivity matrix that defines a graph structure. Estimating connectivity relies on correlation, spectral analysis, and canonic correlation analysis. The matrix is a metric tensor that defines a Riemannian space. The forecasting model is a composition of a graph convolution for aggregating spatial information and a recurrent or neural ODE model.