Bayesian Spike-Triggered Covariance Analysis

Neurons typically respond to a restricted number of stimulus features within the high-dimensional space of natural stimuli. Here we describe an explicit model-based interpretation of traditional estimators for a neuron's multi-dimensional feature space, which allows for several important generalizations and extensions. First, we show that traditional estimators based on the spike-triggered average (STA) and spike-triggered covariance (STC) can be formalized in terms of the "expected log-likelihood" of a Linear-Nonlinear-Poisson (LNP) model with Gaussian stimuli. This model-based formulation allows us to define maximum-likelihood and Bayesian estimators that are statistically consistent and efficient in a wider variety of settings, such as with naturalistic (non-Gaussian) stimuli. It also allows us to employ Bayesian methods for regularization, smoothing, sparsification, and model comparison, and provides Bayesian confidence intervals on model parameters. We describe an empirical Bayes method for selecting the number of features, and extend the model to accommodate an arbitrary elliptical nonlinear response function, which results in a more powerful and more flexible model for feature space inference. We validate these methods using neural data recorded extracellularly from macaque primary visual cortex.

[1]  Max Welling,et al.  Extreme Components Analysis , 2003, NIPS.

[2]  Ian H. Stevenson,et al.  Bayesian Inference of Functional Connectivity and Network Structure From Spikes , 2009, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[3]  S. David,et al.  Estimating sparse spectro-temporal receptive fields with natural stimuli , 2007, Network.

[4]  Eero P. Simoncelli,et al.  Spatiotemporal Elements of Macaque V1 Receptive Fields , 2005, Neuron.

[5]  N. C. Singh,et al.  Estimating spatio-temporal receptive fields of auditory and visual neurons from their responses to natural stimuli , 2001 .

[6]  Sarah M. N. Woolley,et al.  A Generalized Linear Model for Estimating Spectrotemporal Receptive Fields from Responses to Natural Sounds , 2011, PloS one.

[7]  William Bialek,et al.  Real-time performance of a movement-sensitive neuron in the blowfly visual system: coding and information transfer in short spike sequences , 1988, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[8]  Eero P. Simoncelli,et al.  Characterizing Neural Gain Control using Spike-triggered Covariance , 2001, NIPS.

[9]  Eero P. Simoncelli,et al.  Spike-triggered neural characterization. , 2006, Journal of vision.

[10]  E J Chichilnisky,et al.  A simple white noise analysis of neuronal light responses , 2001, Network.

[11]  Maneesh Sahani,et al.  Evidence Optimization Techniques for Estimating Stimulus-Response Functions , 2002, NIPS.

[12]  David P. Wipf,et al.  A New View of Automatic Relevance Determination , 2007, NIPS.

[13]  Kamiar Rahnama Rad,et al.  Efficient, adaptive estimation of two-dimensional firing rate surfaces via Gaussian process methods , 2010, Network.

[14]  L. Paninski Convergence Properties of Some Spike-Triggered Analysis Techniques , 2002 .

[15]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Eero P. Simoncelli,et al.  Dimensionality reduction in neural models: an information-theoretic generalization of spike-triggered average and covariance analysis. , 2006, Journal of vision.

[17]  William Bialek,et al.  Analyzing Neural Responses to Natural Signals: Maximally Informative Dimensions , 2002, Neural Computation.

[18]  Matthias Bethge,et al.  Bayesian Inference for Generalized Linear Models for Spiking Neurons , 2010, Front. Comput. Neurosci..

[19]  Michael E. Tipping,et al.  Probabilistic Principal Component Analysis , 1999 .

[20]  Julian J. Bussgang,et al.  Crosscorrelation functions of amplitude-distorted gaussian signals , 1952 .

[21]  Christopher M. Bishop,et al.  Bayesian PCA , 1998, NIPS.

[22]  Matthias Bethge,et al.  Bayesian Inference for Spiking Neuron Models with a Sparsity Prior , 2007, NIPS.

[23]  Liam Paninski,et al.  Convergence properties of three spike-triggered analysis techniques , 2003, NIPS.

[24]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[25]  de Ruyter van Steveninck,et al.  Real-time performance of a movement-sensitive neuron in the blowfly visual system , 1986 .

[26]  Eero P. Simoncelli,et al.  To appear in: The New Cognitive Neurosciences, 3rd edition Editor: M. Gazzaniga. MIT Press, 2004. Characterization of Neural Responses with Stochastic Stimuli , 2022 .

[27]  L. Paninski Maximum likelihood estimation of cascade point-process neural encoding models , 2004, Network.

[28]  George Eastman House,et al.  Sparse Bayesian Learning and the Relevan e Ve tor Ma hine , 2001 .

[29]  Tom Minka,et al.  Automatic Choice of Dimensionality for PCA , 2000, NIPS.

[30]  Alexander J. Smola,et al.  Neural Information Processing Systems , 1997, NIPS 1997.

[31]  M. Welling,et al.  Bayesian extreme components analysis , 2009, IJCAI 2009.

[32]  C. Stein,et al.  Estimation with Quadratic Loss , 1992 .

[33]  William Bialek,et al.  Synergy in a Neural Code , 2000, Neural Computation.