Statistical Inference for Assessing Functional Connectivity of Neuronal Ensembles With Sparse Spiking Data

Abstract

The ability to accurately infer functional connectivity between ensemble neurons using experimentally acquired spike train data is currently an important research objective in computational neuroscience. Point process generalized linear models and maximum likelihood estimation have been proposed as effective methods for the identification of spiking dependency between neurons. However, unfavorable experimental conditions occasionally results in insufficient data collection due to factors such as low neuronal firing rates or brief recording periods, and in these cases, the standard maximum likelihood estimate becomes unreliable. The present studies compares the performance of different statistical inference procedures when applied to the estimation of functional connectivity in neuronal assemblies with sparse spiking data. Four inference methods were compared: maximum likelihood estimation, penalized maximum likelihood estimation, using either ℓ 2 or ℓ 1 regularization, and hierarchical Bayesian estimation based on a variational Bayes algorithm. Algorithmic performances were compared using well-established goodness-of-fit measures in benchmark simulation studies, and the hierarchical Bayesian approach performed favorably when compared with the other algorithms, and this approach was then successfully applied to real spiking data recorded from the cat motor Correspondence to: Zhe Chen, zhechen@neurostat.mit.edu. cortex. The identification of spiking dependencies in physiologically acquired data was encouraging, since their sparse nature would have previously precluded them from successful analysis using traditional methods. NIH Public Access

A Bayesian nonparametric approach for uncovering rat hippocampal population codes during spatial navigation

“…The goal of Bayesian inference is to incorporate prior knowledge and constraints of the problem and to infer the posterior distribution of unobserved variables of interest (Gelman et al, 2013). In recent years, cutting-edge Bayesian methods have become increasingly popular for data analyses in neuroscience, medicine and biology (Mishchenko et al, 2011; Chen et al, 2011; Chen, 2013; Davidson et al, 2009; Kloosterman et al, 2014; Yau et al, 2011). Specifically, thanks to ever-growing computing power, Markov chain Monte Carlo (MCMC) methods have been widely used in Bayesian inference.…”

Predicting Single-Neuron Activity in Locally Connected Networks

“…One imagines there is a relationship between the temporal filters one attains and the connectivity of the network. Whether one can reverse the process and infer connectivity from predictive power would be of significant interest (Okatan et al, 2005; Nykamp, 2007; Kim, Putrino, Ghosh, & Brown, 2011; Chen, Putrino, Ghosh, Barbieri, & Brown, 2011). Indeed the setting of a reduced computational model, where connectivity and network architecture are under direct control, would be a promising place to begin.…”

Identification of functional synaptic plasticity from spiking activities using nonlinear dynamical modeling

“…Therefore, this MISO model can be termed as a generalized Volterra model (GVM) (Song et al, 2009a, b). As a special case, the first-order GVM is equivalent to the commonly used generalized linear models (Paninski et al, 2004; Okatan et al, 2005; Pillow et al, 2005; Truccolo et al, 2005; Eldawlatly et al, 2009; Chen et al, 2011; Zhao et al, 2012). …”