Characterization of macaque retinal ganglion cell responses using spike-triggered covarianceJ W Pillow, E P Simoncelli and E J ChichilniskyPublished in Computational and Systems Neuroscience (CoSyNe), (109), Mar 2004. |
Multi-electrode recordings from macaque RGCs were obtained from isolated retinas stimulated with one-dimensional spatiotemporal white noise (i.e. flickering bars). Spike-triggered average (STA) analysis revealed center-sur-round spatial organization and biphasic temporal integration expected from RGCs. Eigenvector analysis of the STC revealed components with spatial and temporal structure distinct from the STA. Excitatory components (i.e. those associated with large eigenvalues) exhibited temporal structure similar to the STA, but finer spatial structure, consistent with receptive field subunits. Suppressive components exhibited spatial structure similar to the STA or to excitatory components, but temporal structure that was time-delayed relative to the STA, consistent with spike generation or contrast gain control nonlinearities.
The detailed spatial and temporal structure of these components was probed with a subspace-conditional analysis. STC components were examined in a subspace corresponding to either temporal or spatial stimulus features, obtained by filtering the stimulus with the spatial or temporal profile of the STA. The reduced stimulus dimensionality increased statistical power and revealed a larger number of significant STC components: multiple excitatory spatial components and one or two suppressive temporal components in most cells.
The STC components were then used to constrain a model of RGC light response. The model consisted of identical spatially shifted subunits whose outputs were combined nonlinearly, and a divisive temporal feedback signal. Subunit profiles were constrained to lie within the subspace spanned by spatial and temporal STC components, and the parameters of the nonlinear combination rule and feedback signal were fit using maximum likelihood. This model accurately reproduced the observed STA, contrast-response function, and detailed spatial and temporal structure of STC components. It also provided a substantially more accurate prediction of novel spike trains than a single-filter model characterized using only the STA.