Covariance Regression for High Dimensional Neural Data via Graph

TL;DR

This paper introduces a nonparametric regression framework for modeling high-dimensional neural data, capturing predictor-dependent mean and covariance while incorporating graph structures and Gaussian processes for smoothness.

Contribution

It proposes a novel predictor-dependent covariance regression model with graph-informed structure and Gaussian process smoothness, tailored for high-dimensional neural data with restricted inputs.

Findings

Validated through simulations demonstrating model effectiveness.

Applied to neural datasets showing meaningful covariance patterns.

Flexible framework adaptable to various high-dimensional data applications.

Abstract

Modern recording techniques enable neuroscientists to simultaneously study neural activity across large populations of neurons, with capturing predictor-dependent correlations being a fundamental challenge in neuroscience. Moreover, the fact that input covariates often lie in restricted subdomains, according to experimental settings, makes inference even more challenging. To address these challenges, we propose a set of nonparametric mean-covariance regression models for high-dimensional neural activity with restricted inputs. These models reduce the dimensionality of neural responses by employing a lower-dimensional latent factor model, where both factor loadings and latent factors are predictor-dependent, to jointly model mean and covariance across covariates. The smoothness of neural activity across experimental conditions is modeled nonparametrically using two Gaussian processes (GPs), applied to both loading basis and latent factors. Additionally, to account for the covariates lying in restricted subspace, we incorporate graph information into the covariance structure. To flexibly infer the model, we use an MCMC algorithm to sample from posterior distributions. After validating and studying the properties of proposed methods by simulations, we apply them to two neural datasets (local field potential and neural spiking data) to demonstrate the usage of models for continuous and counting observations. Overall, the proposed methods provide a framework to jointly model covariate-dependent mean and covariance in high dimensional neural data, especially when the covariates lie in restricted domains. The framework is general and can be easily adapted to various applications beyond neuroscience.