Odelia Schwartz, 6/14/01
Spike-triggered average techniques are effective for
linear characterization of neural responses. But sensory neurons
exhibit important nonlinear behaviors that are not captured by such
analyses. Many of these nonlinear behaviors are consistent with gain
control. We demonstrate an automated methodology for characterizing
a neuron with gain control. We assume a model in which the spike rate of
a neuron is proportional to the halfwave rectified response of a linear
kernel suppressively modulated by a weighted sum of rectified responses
of other linear kernels. First, we recover the linear kernel through spike
triggered averaging. Next, we compute the covariance matrix of the stimuli
eliciting spikes, and perform a principal components decomposition of
this matrix. The principal axes (eigenvectors) associated with small
variance (eigenvalue) correspond to directions in which the response
of the neuron is modulated suppressively. Finally, we recover the
suppressive weights associated with these axes by maximizing the
likelihood of the spike data. We demonstrate this method on a simulated
neuron and on salamander retinal ganglion cell data (Chichilnisky lab).
Analysis of physiological data reveals meaningful suppressive axes
and explains interesting nonlinearities. We believe this method will
be applicable to other sensory areas and modalities.