Spike Triggered Covariance in Strongly Correlated Gaussian Stimuli

Many biological systems perform computations on inputs that have very large dimensionality. Determining the relevant input combinations for a particular computation is often key to understanding its function. A common way to find the relevant input dimensions is to examine the difference in variance between the input distribution and the distribution of inputs associated with certain outputs. In systems neuroscience, the corresponding method is known as spike-triggered covariance (STC). This method has been highly successful in characterizing relevant input dimensions for neurons in a variety of sensory systems. So far, most studies used the STC method with weakly correlated Gaussian inputs. However, it is also important to use this method with inputs that have long range correlations typical of the natural sensory environment. In such cases, the stimulus covariance matrix has one (or more) outstanding eigenvalues that cannot be easily equalized because of sampling variability. Such outstanding modes interfere with analyses of statistical significance of candidate input dimensions that modulate neuronal outputs. In many cases, these modes obscure the significant dimensions. We show that the sensitivity of the STC method in the regime of strongly correlated inputs can be improved by an order of magnitude or more. This can be done by evaluating the significance of dimensions in the subspace orthogonal to the outstanding mode(s). Analyzing the responses of retinal ganglion cells probed with Gaussian noise, we find that taking into account outstanding modes is crucial for recovering relevant input dimensions for these neurons.

[1]  D. Thouless,et al.  The eigenvalue spectrum of a large symmetric random matrix with a finite mean , 1978 .

[2]  D. Ringach.  Mapping receptive fields in primary visual cortex , 2004, The Journal of physiology.

[3]  S. Edwards,et al.  The eigenvalue spectrum of a large symmetric random matrix , 1976 .

[4]  V. Plerou,et al.  Random matrix approach to cross correlations in financial data. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  Ryan J. Rowekamp,et al.  Analyzing multicomponent receptive fields from neural responses to natural stimuli , 2011, Network.

[6]  J. Bouchaud,et al.  Noise Dressing of Financial Correlation Matrices , 1998, cond-mat/9810255.

[7]  Inés Samengo,et al.  Spike-triggered covariance: geometric proof, symmetry properties, and extension beyond Gaussian stimuli , 2012, Journal of Computational Neuroscience.

[8]  E. Wigner.  On the Distribution of the Roots of Certain Symmetric Matrices , 1958 .

[9]  N. C. Singh,et al.  Modulation spectra of natural sounds and ethological theories of auditory processing. , 2003, The Journal of the Acoustical Society of America.

[10]  A. Fairhall,et al.  Encoding properties of haltere neurons enable motion feature detection in a biological gyroscope , 2010, Proceedings of the National Academy of Sciences.

[11]  M. Rattray,et al.  Principal-component-analysis eigenvalue spectra from data with symmetry-breaking structure. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  J. Touryan,et al.  Spatial Structure of Complex Cell Receptive Fields Measured with Natural Images , 2005, Neuron.

[13]  Weihua Zhuang,et al.  Joint Power-Frequency-Time Resource Allocation in Clustered Wireless Mesh Networks , 2008, IEEE Network.

[14]  D Marr,et al.  Theory of edge detection , 1979, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[15]  D. Ruderman,et al.  Statistics of cone responses to natural images: implications for visual coding , 1998 .

[16]  Eero P. Simoncelli,et al.  Natural image statistics and neural representation. , 2001, Annual review of neuroscience.

[17]  Kenneth D. Miller,et al.  Adaptive filtering enhances information transmission in visual cortex , 2006, Nature.

[18]  A. Fairhall,et al.  Shifts in Coding Properties and Maintenance of Information Transmission during Adaptation in Barrel Cortex , 2007, PLoS biology.

[19]  C. Atencio,et al.  Hierarchical computation in the canonical auditory cortical circuit , 2009, Proceedings of the National Academy of Sciences.

[20]  William Bialek,et al.  Real-time performance of a movement-sensitive neuron in the blowfly visual system: coding and information transfer in short spike sequences , 1988, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[21]  T. Sharpee,et al.  Estimating linear–nonlinear models using Rényi divergences , 2009, Network.

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

[23]  Feng Qi Han,et al.  Cortical Sensitivity to Visual Features in Natural Scenes , 2005, PLoS biology.

[24]  Michael J. Berry,et al.  Selectivity for multiple stimulus features in retinal ganglion cells. , 2006, Journal of neurophysiology.

[25]  W. Bialek,et al.  Features and dimensions: Motion estimation in fly vision , 2005, q-bio/0505003.

[26]  M. Meister,et al.  Dynamic predictive coding by the retina , 2005, Nature.

[27]  D. Paul.  ASYMPTOTICS OF SAMPLE EIGENSTRUCTURE FOR A LARGE DIMENSIONAL SPIKED COVARIANCE MODEL , 2007 .

[28]  Eero P. Simoncelli,et al.  Spatiotemporal Elements of Macaque V1 Receptive Fields , 2005, Neuron.

[29]  N. C. Singh,et al.  Estimating spatio-temporal receptive fields of auditory and visual neurons from their responses to natural stimuli , 2001 .

[30]  Noureddine El Karoui.  Spectrum estimation for large dimensional covariance matrices using random matrix theory , 2006, math/0609418.

[31]  Katherine I. Nagel,et al.  Two-dimensional adaptation in the auditory forebrain. , 2011, Journal of neurophysiology.

[32]  Eero P. Simoncelli,et al.  Dimensionality reduction in neural models: an information-theoretic generalization of spike-triggered average and covariance analysis. , 2006, Journal of vision.

[33]  Tatyana O. Sharpee,et al.  Second Order Dimensionality Reduction Using Minimum and Maximum Mutual Information Models , 2011, PLoS Comput. Biol..

[34]  Najeeb M. Halabi,et al.  Protein Sectors: Evolutionary Units of Three-Dimensional Structure , 2009, Cell.

[35]  Aurel A. Lazar,et al.  System identification of Drosophila olfactory sensory neurons , 2011, Journal of Computational Neuroscience.

[36]  D. Ruderman,et al.  Independent component analysis of natural image sequences yields spatio-temporal filters similar to simple cells in primary visual cortex , 1998, Proceedings of the Royal Society of London. Series B: Biological Sciences.

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

[38]  William Bialek,et al.  Maximally Informative “Stimulus Energies” in the Analysis of Neural Responses to Natural Signals , 2012, PloS one.

[39]  Il Memming Park,et al.  Bayesian Spike-Triggered Covariance Analysis , 2011, NIPS.

[40]  Eero P. Simoncelli,et al.  Spike-triggered neural characterization. , 2006, Journal of vision.

[41]  William Bialek,et al.  Statistics of Natural Images: Scaling in the Woods , 1993, NIPS.

[42]  Michael J. Berry,et al.  Recording spikes from a large fraction of the ganglion cells in a retinal patch , 2004, Nature Neuroscience.

[43]  C. Atencio,et al.  Cooperative Nonlinearities in Auditory Cortical Neurons , 2008, Neuron.

[44]  Todd M. Allen,et al.  Coordinate linkage of HIV evolution reveals regions of immunological vulnerability , 2011, Proceedings of the National Academy of Sciences.

[45]  T. Sejnowski,et al.  Color opponency is an efficient representation of spectral properties in natural scenes , 2002, Vision Research.

[46]  J. Touryan,et al.  Isolation of Relevant Visual Features from Random Stimuli for Cortical Complex Cells , 2002, The Journal of Neuroscience.

[47]  D J Field,et al.  Relations between the statistics of natural images and the response properties of cortical cells. , 1987, Journal of the Optical Society of America. A, Optics and image science.