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.

Variable Selection Using Nonlocal Priors in High-Dimensional Generalized Linear Models With Application to fMRI Data Analysis

Bayesian variable selection using spike‐and‐slab priors with application to high dimensional electroencephalography data by local modelling

Variable Selection for Generalized Varying Coefficient Partially Linear Models with Diverging Number of Parameters

Variable Selection and Minimax Prediction in High-dimensional Functional Linear Model

Bayesian variable selection in linear regression models with instrumental variables

Variable Selection for High-dimensional Generalized Linear Models using an Iterated Conditional Modes/Medians Algorithm

A Selective Overview of Variable Selection in High Dimensional Feature Space (Invited Review Article)

Variable selection for partially linear models via Bayesian subset modeling with diffusing prior

A High-dimensional M-estimator Framework for Bi-level Variable Selection

Variable selection in high-dimensional varying-coefficient models with global optimality

Feature screening and variable selection for partially linear models with ultrahigh-dimensional longitudinal data

Regularization Methods for High-Dimensional Instrumental Variables Regression With an Application to Genetical Genomics

Ultrahigh dimensional variable selection: beyond the linear model

A robust and efficient estimation and variable selection method for partially linear models with large-dimensional covariates

Bayesian Variable Selection and Estimation Based on Global-Local Shrinkage Priors

Thresholded Multiscale Gaussian Processes with Application to Bayesian Feature Selection for Massive Neuroimaging Data

Advances in Bayesian model selection consistency for high-dimensional generalized linear models

Variable Selection for Latent Dirichlet Allocation

Nonparametric Covariance Regression for Massive Neural Data on Restricted Covariates via Graph

On the Generalizability of Linear and Non-Linear Region of Interest-Based Multivariate Regression Models for fMRI Data

Nonstationary Gaussian Process Discriminant Analysis With Variable Selection for High-Dimensional Functional Data