Maximum Likelihood Estimation of a Stochastic Integrate-and-Fire Neural Model

Recent work has examined the estimation of models of stimulus-driven neural activity in which some linear filtering process is followed by a nonlinear, probabilistic spiking stage. We analyze the estimation of one such model for which this nonlinear step is implemented by a noisy, leaky, integrate-and-fire mechanism with a spike-dependent aftercurrent. This model is a biophysically plausible alternative to models with Poisson (memory-less) spiking, and has been shown to effectively reproduce various spiking statistics of neurons in vivo. However, the problem of estimating the model from extracellular spike train data has not been examined in depth. We formulate the problem in terms of maximum likelihood estimation, and show that the computational problem of maximizing the likelihood is tractable. Our main contribution is an algorithm and a proof that this algorithm is guaranteed to find the global optimum with reasonable speed. We demonstrate the effectiveness of our estimator with numerical simulations.

[1]  Wulfram Gerstner,et al.  Spiking Neuron Models: An Introduction , 2002 .

[2]  Adrienne L. Fairhall,et al.  What Causes a Neuron to Spike? , 2003, Neural Computation.

[3]  Yuguo Yu,et al.  Dynamical mechanisms underlying contrast gain control in single neurons. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Emery N. Brown,et al.  The Time-Rescaling Theorem and Its Application to Neural Spike Train Data Analysis , 2002, Neural Computation.

[5]  R. Reid,et al.  Predicting Every Spike A Model for the Responses of Visual Neurons , 2001, Neuron.

[6]  John P. Miller,et al.  Broadband neural encoding in the cricket cereal sensory system enhanced by stochastic resonance , 1996, Nature.

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

[8]  Eero P. Simoncelli,et al.  Biases in white noise analysis due to non-Poisson spike generation , 2003, Neurocomputing.

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

[10]  B. Knight,et al.  The Power Ratio and the Interval Map: Spiking Models and Extracellular Recordings , 1998, The Journal of Neuroscience.

[11]  Eero P. Simoncelli,et al.  Maximum Likelihood Estimation of a Stochastic Integrate-and-Fire Neural Encoding Model , 2004, Neural Computation.

[12]  Wulfram Gerstner,et al.  SPIKING NEURON MODELS Single Neurons , Populations , Plasticity , 2002 .

[13]  A. Genz.  Numerical Computation of Multivariate Normal Probabilities , 1992 .

[14]  Liam Paninski,et al.  Noise-driven adaptation: in vitro and mathematical analysis , 2003, Neurocomputing.

[15]  Lawrence Sirovich,et al.  The Approach of a Neuron Population Firing Rate to a New Equilibrium: An Exact Theoretical Result , 2000, Neural Computation.

[16]  Michael J. Berry,et al.  Refractoriness and Neural Precision , 1997, The Journal of Neuroscience.

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

[18]  Jonathan D. Victor,et al.  How the brain uses time to represent and process visual information 1 1 Published on the World Wide Web on 16 August 2000. , 2000, Brain Research.

[19]  Michael E. Rudd,et al.  Noise Adaptation in Integrate-and-Fire Neurons , 1997, Neural Computation.

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .

[21]  S. Karlin,et al.  A second course in stochastic processes , 1981 .