Identification of neural feature space from spike triggered covariance expressed as a function of PRC
BMC Neuroscience volume 11, Article number: P14 (2010)
For the purpose of elucidating the neural coding process based on the neural excitability mechanism, some researchers have investigated the relationship between the neural dynamics and the spike triggered stimulus ensemble (STE), which indicates what stimuli are more likely or less likely to induce neural spikes. Ermentrout et al. have analytically derived the relational equation between the phase response curve (PRC) and the spike triggered average (STA), which is the average of the STE, when regular spikes with a period are disturbed by sufficiently small white noise, as (1). Here, is the time relative to a spike, is the noise intensity, and is PRC . Furthermore, they showed that Eq. (1) holds true for real neurons. Their study has made meaningful progress in relating the neural dynamics to the neural coding for real neurons. However, the STA is the first cumulant of the STE. In order to approximately identify the distribution of STE as a Gaussian, we should determine its second cumulant, called spike triggered covariance (STC).
where represents the Heaviside function which takes 1/2 at . Moreover, we analyze the eigenfunctions of in order to extract the neural feature space, which is a low dimensional subspace of the full stimulus space characterizing the stimulus encoded by neurons. The eigenfunctions associated with the positive and negative eigenvalues of are called the excitatory and suppressive eigenfunction, respectively. In this case, the stimuli in the subspace spanned by excitatory eigenfunctions cause shorter interspike intervals (ISIs) than , while the stimuli in the subspace spanned by suppressive eigenfunctions cause longer ISIs.
Figure 1 shows the STC of a rat hippocampal CA1 pyramidal neuron as calculated by Eq. (2), where the PRC could be estimated by our algorithm . Note that it is difficult to measure the STC for real neurons directly, because the number of neural spikes required for a stable calculation of STC is nearly square of the number required for the STA. Figure 1 suggests that the neural feature space of this rat hippocampal CA1 pyramidal neuron can be described by the four eigenfunctions in Fig. 1b.
Ermentrout GB, Galan RF, Urban NN: Relating neural dynamics to neural coding. Phys Rev Lett. 2007, 99: 248103-10.1103/PhysRevLett.99.248103.
Hong S, Arcas BA, Fairhall AL: Single Neuron Computation: From Dynamical System to Feature Detector. Neural Comput. 2007, 19 (12): 3133-3172. 10.1162/neco.2007.19.12.3133.
Teramae J, Fukai T: Temporal Precision of Spike Response to Fluctuating Input in Pulse-Coupled Networks of Oscillating Neurons. Phys Rev Lett. 2008, 101: 248105-10.1103/PhysRevLett.101.248105.
Ota K, Tsunoda T, Omori T, Watanabe S, Miyakawa H, Okada M, Aonishi T: Is the Langevin phase equation an efficient model for oscillating neurons?. J Phys Conf. 2009, 197: 012016-10.1088/1742-6596/197/1/012016.
About this article
Cite this article
Ota, K., Omori, T., Watanabe, S. et al. Identification of neural feature space from spike triggered covariance expressed as a function of PRC. BMC Neurosci 11 (Suppl 1), P14 (2010). https://doi.org/10.1186/1471-2202-11-S1-P14