Abstract:The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μ ECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets. Extracting insight from data into the dynamic processes that produce observed phenomena is a ubiquitous challenge in biology. For example, brain functions are generated by the collective activity of many neurons. Extracting insight from biological time-series data can be challenging because the number of observations can be large, activities can be noisy, and collective dynamics can be complex. Here, we developed a new latent processes Bayesian regression model, the orthogonal stochastic linear mixing model (OSLMM) to address these challenges. Compared to competing methods, we demonstrate that OSLMM has superior recovery performance of latent trajectories in synthetic experiments and is computationally efficient. In two diverse neural data sets, we find that the latent spaces extracted from OSLMM are directly structured by task parameters. Our results demonstrate that latent dynamics can be structured by representations and suggest that OSLMM will be useful for data-driven discovery in large-scale biological time-series data.

Brain Waves Analysis Via a Non-Parametric Bayesian Mixture of Autoregressive Kernels

Modeling Spectral Properties in Stationary Processes of Varying Dimensions with Applications to Brain Local Field Potential Signals

Full Bayesian Modeling for fMRI Group Analysis

Bayesian Kernel Methods for Analysis of Functional Neuroimages

Bayesian Model Selection Maps for Group Studies Using M/EEG Data

Bayesian Nonparametric Models for Synchronous Brain-Computer Interfaces

BSD: a Bayesian framework for parametric models of neural spectra

A hierarchical Bayesian approach for learning sparse spatio-temporal decompositions of multichannel EEG.

A machine learning approach for real-time cortical state estimation

Bayesian analysis of phase data in EEG and MEG

Evolutionary State-Space Model and Its Application to Time-Frequency Analysis of Local Field Potentials

Spatio-Spectral Analysis of Brain Signals

Multiscale Bayesian State Space Model for Granger Causality Analysis of Brain Signal

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

Hierarchical Bayesian Modeling of Inter-Trial Variability and Variational Bayesian Learning of Common Spatial Patterns from Multichannel EEG

Nonlinear brain dynamics as macroscopic manifestation of underlying many-body field dynamics

Recovering Arrhythmic EEG Transients from Their Stochastic Interference

Bayesian Analysis of fMRI data with Spatially-Varying Autoregressive Orders

Functional-Coefficient Models for Multivariate Time Series in Designed Experiments: with Applications to Brain Signals

Mapping single-trial EEG records on the cortical surface through a spatiotemporal modality

Symposium Title: Computational Psychophysiology