Spike-triggered covariance

Spike-triggered covariance (STC) analysis is a tool for characterizing a neuron's response properties using the covariance of stimuli that elicit spikes from a neuron. STC is related to the spike-triggered average (STA), and provides a complementary tool for estimating linear filters in a linear-nonlinear-Poisson (LNP) cascade model. Unlike STA, the STC can be used to identify a multi-dimensional feature space in which a neuron computes its response.

STC analysis identifies the stimulus features affecting a neuron's response via an eigenvector decomposition of the spike-triggered covariance matrix.[1][2][3][4] Eigenvectors with eigenvalues significantly larger or smaller than the eigenvalues of the raw stimulus covariance correspond to stimulus axes along which the neural response is enhanced or suppressed.

STC analysis is similar to principal components analysis (PCA), although it differs in that the eigenvectors corresponding to largest and smallest eigenvalues are used for identifying the feature space. The STC matrix is also known as the 2nd-order Volterra or Wiener kernel.

Mathematical definition

edit

Standard STC

edit

Let   denote the spatio-temporal stimulus vector preceding the  'th time bin, and   the spike count in that bin. The stimuli can be assumed to have zero mean (i.e.,  ). If not, it can be transformed to have zero-mean by subtracting the mean stimulus from each vector. The spike-triggered covariance (STC) is given by

 

where   is the total number of spikes, and STA is the spike-triggered average. The covariance of the stimulus is given by

 

where   is the number of stimuli   used during the experiment. The eigenvectors of   associated to significantly positive eigenvalues correspond to excitatory vectors, whereas eigenvectors associated to significantly negative eigenvalues are inhibitory eigenvectors.[4]

References

edit
  1. ^ Brenner, N., Bialek, W., & de Ruyter van Steveninck, R.R. (2000).
  2. ^ Schwartz, O., Chichilnisky, E. J., & Simoncelli, E. P. (2002).
  3. ^ Bialek, W. & de Ruyter van Steveninck, R. (2005). Arxiv preprint q-bio/0505003.
  4. ^ a b Schwartz O., Pillow J. W., Rust N. C., & Simoncelli E. P. (2006). Spike-triggered neural characterization. Journal of Vision 6:484-507
edit

Matlab code for STA/STC analysis of neural data