imported>Jonathan Tweet: the lead should summarize the topic, adding information from the body. see wp:lead

2020-06-14T22:33:16Z

the lead should summarize the topic, adding information from the body. see wp:lead

New page

The '''linear-nonlinear-Poisson (LNP) cascade model''' is a simplified functional model of neural spike responses.<ref name="chichilnisky01">Chichilnisky, E. J., [http://www.stat.columbia.edu/~liam/teaching/neurostat-fall13/papers/Chichilnisky-2001.pdf A simple white noise analysis of neuronal light responses.] {{webarchive|url=https://web.archive.org/web/20081007232914/http://www.stat.columbia.edu/~liam/teaching/neurostat-spr07/papers/Chichilnisky-2001.pdf |date=2008-10-07 }} Network: Computation in Neural Systems 12:199–213. (2001)</ref><ref name="simoncelli">Simoncelli, E. P., Paninski, L., Pillow, J. & Swartz, O. (2004). [http://www.cns.nyu.edu/~eero/ABSTRACTS/simoncelli03c-abstract.html Characterization of Neural Responses with Stochastic Stimuli] in (Ed. M. Gazzaniga) ''The Cognitive Neurosciences 3rd edn'' (pp 327–338) MIT press.</ref><ref name="schwartz06">Schwartz O., Pillow J. W., Rust N. C., & Simoncelli E. P. (2006). Spike-triggered neural characterization. ''Journal of Vision'' 6:484–507</ref> It has been successfully used to describe the response characteristics of neurons in early sensory pathways, especially the visual system. The LNP model is generally implicit when using reverse correlation or the [[spike-triggered average]] to characterize neural responses with white-noise stimuli.

[[Image:LNPModelDiagram.png|right|thumb|300px|The Linear-Nonlinear-Poisson Cascade Model]]
There are three stages of the LNP cascade model. The first stage consists of a linear filter, or linear [[receptive field]], which describes how the neuron integrates stimulus intensity over space and time. The output of this filter then passes through a nonlinear function, which gives the neuron's instantaneous spike rate as its output. Finally, the spike rate is used to generate spikes according to an inhomogeneous [[Poisson process]].

The linear filtering stage performs [[dimensionality reduction]], reducing the high-dimensional spatio-temporal stimulus space to a low-dimensional [[feature space]], within which the neuron computes its response. The nonlinearity converts the filter output to a (non-negative) spike rate, and accounts for nonlinear phenomena such as spike threshold (or rectification) and response saturation. The Poisson spike generator converts the continuous spike rate to a series of spike times, under the assumption that the probability of a spike depends only on the instantaneous spike rate.

The model offers a useful approximation of neural activity, allowing scientists to derive reliable estimates from a mathematically simple formula.

== Mathematical formulation ==

===Single-filter LNP===
Let <math>\mathbf{x}</math> denote the spatio-temporal stimulus vector at a particular instant, and
<math>\mathbf{k}</math> denote a linear filter (the neuron's linear receptive field), which is a vector with the same number of elements as <math>\mathbf{x}</math>. Let <math>f</math> denote the nonlinearity, a scalar function with non-negative output. Then the LNP model specifies that, in the limit of small time bins,
: <math>P(\textrm{spike}) \propto f(\mathbf{k} \cdot \mathbf{x})</math>.
For finite-sized time bins, this can be stated precisely as the probability of observing ''y'' spikes in a single bin:
: <math> P(y \textrm{~spikes}) = \frac{\left(\Delta \lambda\right)^y}{y!} e^{- \Delta \lambda}</math>
: where <math> \lambda = f(\mathbf{k}\cdot\mathbf{x})</math>, and <math>\Delta</math> is the bin size.

===Multi-filter LNP===
For neurons sensitive to multiple dimensions of the stimulus space, the linear stage of the LNP model can be generalized to a bank of linear filters, and the nonlinearity becomes a function of multiple inputs. Let <math>\mathbf{k_1}, \mathbf{k_2}, \ldots, \mathbf{k_n}</math> denote the set of linear filters that capture a neuron's stimulus dependence. Then the multi-filter LNP model is described by
: <math>P(\textrm{spike}) \propto f(\mathbf{k_1}\!\cdot\!\mathbf{x},\; \mathbf{k_2}\!\cdot\!\mathbf{x},\; \ldots,\; \mathbf{k_n}\!\cdot\!\mathbf{x})</math>
or
: <math>P(\textrm{spike}) \propto f(K\mathbf{x}),</math>
where <math>K</math> is a matrix whose columns are the filters <math>\mathbf{k_i}</math>.

== Estimation ==
The parameters of the LNP model consist of the linear filters <math>\{{k_i}\}</math> and the nonlinearity <math>f</math>. The estimation problem (also known as the problem of ''neural characterization'') is the problem of determining these parameters from data consisting of a time-varying stimulus and the set of observed spike times. Techniques for estimating the LNP model parameters include:
* moment-based techniques, such as the [[spike-triggered average]] or [[spike-triggered covariance]]<ref name="chichilnisky01"/><ref name="simoncelli"/><ref name="schwartz06"/><ref name="brenner00">Brenner, N., Bialek, W., & de Ruyter van Steveninck, R. R. (2000).</ref>
* with information-maximization or [[maximum likelihood]] techniques.<ref name="Paninski:2004">Paninski, L. (2004) Maximum likelihood estimation of cascade point-process neural encoding models. In ''Network: Computation in Neural Systems''.</ref><ref name="Mirbagheri:2012">Mirbagheri M. (2012) [https://ieeexplore.ieee.org/document/6288342/?arnumber=6288342&tag=1 Dimension reduction in regression using Gaussian Mixture Models.] In ''Proceedings of International Conference on Acoustics, Speech and Signal Processing (ICASSP)''.</ref>

== Related models ==
* The LNP model provides a simplified, mathematically tractable approximation to more biophysically detailed [[Biological neuron model|single-neuron models]] such as the [[integrate-and-fire]] or [[Hodgkin–Huxley model]].
* If the nonlinearity <math>f</math> is a fixed invertible function, then the LNP model is a [[generalized linear model]]. In this case, <math>f</math> is the inverse link function.
* An alternative to the LNP model for neural characterization is the [[Volterra series|Volterra kernel]] or [[Wiener kernel]] series expansion, which arises in classical nonlinear systems-identification theory.<ref>Marmarelis & Marmerelis, 1978. ''Analysis of Physiological Systems: The White Noise Approach.'' London: Plenum Press.</ref> These models approximate a neuron's input-output characteristics using a polynomial expansion analogous to the [[Taylor series]], but do not explicitly specify the spike-generation process.

==See also==
*[[Random neural network]]
*[[Spike-triggered average]]
*[[Spike-triggered covariance]]

==References==
<references/>

[[Category:Computational neuroscience]]
[[Category:Stochastic models]]

Linear-nonlinear-Poisson cascade model - Revision history

imported>Jonathan Tweet: the lead should summarize the topic, adding information from the body. see wp:lead