2pm, Tuesday, 18 Oct 2005:
Dimensionality reduction in neural models: An information-theoretic generalization of spike- triggered average and covariance analysis

Jonathan Pillow
NYU

We describe an information-theoretic framework for characterizing spiking neural responses using a Linear-Nonlinear-Poisson (LNP) cascade model. The linear stage is recovered by sequentially searching for a set of linear filters that maximize the information between stimuli and spiking responses. This optimization problem is made tractable by assuming that both the raw and spike-triggered stimulus ensembles are drawn from Gaussian distributions. The resulting solution has several surprising properties: (1) it is equivalent to maximizing the likelihood of the spike train given the stimuli; (2) like many existing techniques, it is based on initial measurements of spike-triggered average and covariance, but it unifies these in a common framework that provides a set of filters sorted according to their informativeness about the neural response; (3) unlike unrestricted information-theoretic estimators, it is both computationally efficient and robust, and the data requirements for recovery of a linear stage of high dimensionality are relatively modest; (4) it implicitly assumes that the nonlinear stage maps the linear filter responses to a Poisson spike rate using a ratio-of-Gaussians function, which is capable of mimicking neural response in the early visual system. We demonstrate the effectiveness and advantages of the method by applying it to simulated responses of a Hodgkin-Huxley neuron, and the recorded extracelluar responses of macaque retinal ganglion cells and V1 cells.