Analyzing multiple spike trains with nonparametric granger causality

Simultaneous recordings of spike trains from multiple single neurons are becoming commonplace. Understanding the interaction patterns among these spike trains remains a key research area. A question of interest is the evaluation of information flow between neurons through the analysis of whether one spike train exerts causal influence on another. For continuous-valued time series data, Granger causality has proven an effective method for this purpose. However, the basis for Granger causality estimation is autoregressive data modeling, which is not directly applicable to spike trains. Various filtering options distort the properties of spike trains as point processes. Here we propose a new nonparametric approach to estimate Granger causality directly from the Fourier transforms of spike train data. We validate the method on synthetic spike trains generated by model networks of neurons with known connectivity patterns and then apply it to neurons simultaneously recorded from the thalamus and the primary somatosensory cortex of a squirrel monkey undergoing tactile stimulation.     

Estimating summary statistics in the spike-train space

Estimating sample averages and sample variability is important in analyzing neural spike trains data in computational neuroscience. Current approaches have focused on advancing the use of parametric or semiparametric probability models of the underlying stochastic process, where the probabilistic distribution is characterized at each time point with basic statistics such as mean and variance. To directly capture and analyze the average and variability in the observation space of the spike trains, we focus on a data-driven approach where statistics are defined and computed in a function space in which the spike trains are viewed as individual points. Based on the definition of a “Euclidean” metric, a recent paper introduced the notion of the mean of a set of spike trains and developed an efficient algorithm to compute it under some restrictive conditions. Here we extend this study by: (1) developing a novel algorithm for mean computation that is quite general, and (2) introducing a notion of covariance of a set of spike trains. Specifically, we estimate the covariance matrix using the geometry of the warping functions that map the mean spike train to each of the spike trains in the dataset. Results from simulations as well as a neural recording in primate motor cortex indicate that the proposed mean and covariance successfully capture the observed variability in spike trains. In addition, a “Gaussian-type” probability model (defined using the estimated mean and covariance) reasonably characterizes the distribution of the spike trains and achieves a desirable performance in the classification of the spike trains.     

Quantification of Clustering in Joint Interspike Interval Scattergrams of Spike Trains

Joint interval scattergrams are usually employed in determining serial correlations between events of spike trains. However, any inherent structures in such scattergrams that are often seen in experimental records are not quantifiable by serial correlation coefficients. Here, we develop a method to quantify clustered structures in any two-dimensional scattergram of pairs of interspike intervals. The method gives a cluster coefficient as well as clustering density function that could be used to quantify clustering in scattergrams obtained from first- or higher-order interval return maps of single spike trains, or interspike interval pairs drawn from simultaneously recorded spike trains. The method is illustrated using numerical spike trains as well as in vitro pairwise recordings of rat striatal tonically active neurons.     

Factors Influencing the Estimates of Correlation between Motor Unit Activities in Humans

Background

Alpha motoneurons receive common synaptic inputs from spinal and supraspinal pathways. As a result, a certain degree of correlation can be observed between motoneuron spike trains during voluntary contractions. This has been studied by using correlation measures in the time and frequency domains. These measures are interpreted as reflecting different types of connectivity in the spinal networks, although the relation between the degree of correlation of the output motoneuron spike trains and of their synaptic inputs is unclear.

Methodology/Principal Findings

In this study, we analyze theoretically this relation and we complete this analysis by simulations and experimental data on the abductor digiti minimi muscle. The results demonstrate that correlation measures between motoneuron output spike trains are inherently influenced by the discharge rate and that this influence cannot be compensated by normalization. Because of the influence of discharge rate, frequency domain measures of correlation (coherence) do not identify the full frequency content of the common input signal when computed from pairs of motoneurons. Rather, an increase in sampling rate is needed by using cumulative spike trains of several motoneurons. Moreover, the application of averaging filters to the spike trains influences the magnitude of the estimated correlation levels calculated in the time, but not in the frequency domain (coherence).

Conclusions

It is concluded that the analysis of coherence in different frequency bands between cumulative spike trains of a sufficient number of motoneurons provides information on the spectrum of the common synaptic input. Nonetheless, the absolute values of coherent peaks cannot be compared across conditions with different cumulative discharge rates.     

Stimulus dependent neural correlations in the auditory midbrain of the grassfrog (Rana temporaria L.)

Few-unit recordings were obtained using metal microelectrodes. Separation into single-unit spike trains was based on differences in spike amplitude and spike waveform. For that purpose a hardware microprocessor based spike waveform analyser was designed and built. Spikes are filtered by four matched filters and filter outputs at the moments of spike occurrence are read by a computer and used for off-line separation and spike waveform reconstruction. Thirty-one double unit recordings were obtained and correlation between the separated spike trains was determined. After stimulus correction correlation remained in only 8 of the double unit records. It appeared that in most cases this neural correlation was stimulus dependent. Continuous noise stimulation resulted in the strongest neural correlation remaining after correction for stimulus coupling, stimulation with 48 ms duration tonepips presented once per second generally did not result in a significant neural correlation after the correction procedure for stimulus lock. The usefulness of the additive model for neural correlation and the correction procedure based thereupon is discussed.     

Neural coding properties based on spike timing and pattern correlation of retinal ganglion cells

Correlation between spike trains or neurons sometimes indicates certain neural coding rules in the visual system. In this paper, the relationship between spike timing correlation and pattern correlation is discussed, and their ability to represent stimulus features is compared to examine their coding strategies not only in individual neurons but also in population. Two kinds of stimuli, natural movies and checkerboard, are used to arouse firing activities in chicken retinal ganglion cells. The spike timing correlation and pattern correlation are calculated by cross-correlation function and Lempel–Ziv distance respectively. According to the correlation values, it is demonstrated that spike trains with similar spike patterns are not necessarily concerted in firing time. Moreover, spike pattern correlation values between individual neurons’ responses reflect the difference of natural movies and checkerboard; neurons cooperate with each other with higher pattern correlation values which represent spatiotemporal correlations during response to natural movies. Spike timing does not reflect stimulus features as obvious as spike patterns, caused by their particular coding properties or physiological foundation. As a result, separating the pattern correlation out of traditional timing correlation concept uncover additional insight in neural coding.     

Multi-scale detection of rate changes in spike trains with weak dependencies

The statistical analysis of neuronal spike trains by models of point processes often relies on the assumption of constant process parameters. However, it is a well-known problem that the parameters of empirical spike trains can be highly variable, such as for example the firing rate. In order to test the null hypothesis of a constant rate and to estimate the change points, a Multiple Filter Test (MFT) and a corresponding algorithm (MFA) have been proposed that can be applied under the assumption of independent inter spike intervals (ISIs). As empirical spike trains often show weak dependencies in the correlation structure of ISIs, we extend the MFT here to point processes associated with short range dependencies. By specifically estimating serial dependencies in the test statistic, we show that the new MFT can be applied to a variety of empirical firing patterns, including positive and negative serial correlations as well as tonic and bursty firing. The new MFT is applied to a data set of empirical spike trains with serial correlations, and simulations show improved performance against methods that assume independence. In case of positive correlations, our new MFT is necessary to reduce the number of false positives, which can be highly enhanced when falsely assuming independence. For the frequent case of negative correlations, the new MFT shows an improved detection probability of change points and thus, also a higher potential of signal extraction from noisy spike trains.     

Sequential interval histogram analysis of non-stationary neuronal spike trains

The spike interval histogram, a commonly used tool for the analysis of neuronal spike trains, is evaluated as a statistical estimator of the probability density function (pdf) ofinterspike intervals. Using a mean square error criterion, it is concluded that a Parzen convolution estimate of the pdf is superior to the conventional histogram procedure. The Parzen estimate using a Gaussian weighting function reduces the number of intervals required to achieve a given error by a factor of 5–10. The Parzen estimation procedure has been implemented in the sequential interval histogram (SQIH) procedure for analysis of non-stationary spike trains. Segments of the spike train are defined using a moving window and the pdf for each segment is estimated sequentially. The procedure which we have found most practical is interactive with the user and utlizes the theoretical results of the error analysis as guidelines for the evolution of an estimation strategy. The SQIH procedure appears useful both as a criterion for stationarity and as a means to characterize non-stationary activity.Portions of this work were presented at the Symposium on Computer Technology in Neuroscience Research, West Virginia University Medical Center, Morgantown, West Virginia, USA, April, 1975.     

Identification and estimation algorithm for stochastic neural system

An algorithm for the estimation of stochastic processes in a neural system is presented. This process is defined here as the continuous stochastic process reflecting the dynamics of the neural system which has some inputs and generates output spike trains. The algorithm proposed here is to identify the system parameters and then estimate the stochastic process called neural system process here. These procedures carried out on the basis of the output spike trains which are supposed to be the data observed in the randomly missing way by the threshold time function in the neural system. The algorithm is constructed with the well-known Kalman filters and realizes the estimation of the neural system process by cooperating with the algorithm for the parameter estimation of the threshold time function presented previously (Nakao et al., 1983). The performance of the algorithm is examined by applying it to the various spike trains simulated by some artificial models and also to the neural spike trains recorded in cat's optic tract fibers. The results in these applications are thought to prove the effectiveness of the algorithm proposed here to some extent. Such attempts, we think, will serve to improve the characterizing and modelling techniques of the stochastic neural systems.     

Changes in interspike interval during propagation: Quantitative description

The intervals between nerve impulses can change substantially during propagation because of conduction velocity aftereffects of previous impulse activity. Effects of such changes on interval histograms and on statistical parameters of spike trains were evaluated for Poisson spike trains and for trains generated by a clock with added random delays. The distribution of short intervals was significantly changed during propagation for these spike trains. Substantial changes in serial correlation coefficients were found in trains with certain initial interval distributions. The relevance of these effects to neural coding is discussed.     

Influence of temporal correlation of synaptic input on the rate and variability of firing in neurons

The spike trains that transmit information between neurons are stochastic. We used the theory of random point processes and simulation methods to investigate the influence of temporal correlation of synaptic input current on firing statistics. The theory accounts for two sources for temporal correlation: synchrony between spikes in presynaptic input trains and the unitary synaptic current time course. Simulations show that slow temporal correlation of synaptic input leads to high variability in firing. In a leaky integrate-and-fire neuron model with spike afterhyperpolarization the theory accurately predicts the firing rate when the spike threshold is higher than two standard deviations of the membrane potential fluctuations. For lower thresholds the spike afterhyperpolarization reduces the firing rate below the theory's predicted level when the synaptic correlation decays rapidly. If the synaptic correlation decays slower than the spike afterhyperpolarization, spike bursts can occur during single broad peaks of input fluctuations, increasing the firing rate over the prediction. Spike bursts lead to a coefficient of variation for the interspike intervals that can exceed one, suggesting an explanation of high coefficient of variation for interspike intervals observed in vivo.     

The relationship of post-stimulus time and interval histograms to the timming characteristics of spike trains.

PST (post-stimulus time) and interval histograms computed from recorded spike trains are related to an average timing characteristics of the spike train. The exact nature of this relationship varies with recording parameters, interfering signals, the histogram bin width, and the duration of the measurement interval. This work describes the conditions under which a PST histogram can serve as an unbiased estimate of the ensemble average of a spike train's intensity and an interval histogram can serve as an unbiased estimate of the probability density function of the interspike intervals. Simulation studies are used to confirm the validity of the theoretical results. As an example of an application, these results are used to analyze recordings of singleunit activity in the eight cranial nerve.     

Alias-free sampling of neuronal spike trains

Summary Spectral analysis provides powerful techniques for describing the lower order moments of a stochastic process and interactions between two or more stochastic processes. A major problem in the application of spectral analysis to neuronal spike trains is how to obtain equispaced samples of the spike trains which will give unbiased and alias-free spectral estimates. Various sampling methods, which treat the spike train as a continuous signal, a point process and as a series of Dirac delta-functions, are reviewed and their limitations discussed. A new sampling technique, which gives unbiased and alias-free estimates, is described. This technique treats the spike train as a series of delta functions and generates samples by digital filtering. Implementation of this technique on a small computer is simple and virtually on-line.     

A model for the neuronal implementation of selective visual attention based on temporal correlation among neurons

We propose a model for the neuronal implementation of selective visual attention based on temporal correlation among groups of neurons. Neurons in primary visual cortex respond to visual stimuli with a Poisson distributed spike train with an appropriate, stimulus-dependent mean firing rate. The spike trains of neurons whose receptive fields donot overlap with the focus of attention are distributed according to homogeneous (time-independent) Poisson process with no correlation between action potentials of different neurons. In contrast, spike trains of neurons with receptive fields within the focus of attention are distributed according to non-homogeneous (time-dependent) Poisson processes. Since the short-term average spike rates of all neurons with receptive fields in the focus of attention covary, correlations between these spike trains are introduced which are detected by inhibitory interneurons in V4. These cells, modeled as modified integrate-and-fire neurons, function as coincidence detectors and suppress the response of V4 cells associated with non-attended visual stimuli. The model reproduces quantitatively experimental data obtained in cortical area V4 of monkey by Moran and Desimone (1985).     

Variability and coding efficiency of noisy neural spike encoders

Encoding synaptic inputs as a train of action potentials is a fundamental function of nerve cells. Although spike trains recorded in vivo have been shown to be highly variable, it is unclear whether variability in spike timing represents faithful encoding of temporally varying synaptic inputs or noise inherent in the spike encoding mechanism. It has been reported that spike timing variability is more pronounced for constant, unvarying inputs than for inputs with rich temporal structure. This could have significant implications for the nature of neural coding, particularly if precise timing of spikes and temporal synchrony between neurons is used to represent information in the nervous system. To study the potential functional role of spike timing variability, we estimate the fraction of spike timing variability which conveys information about the input for two types of noisy spike encoders--an integrate and fire model with randomly chosen thresholds and a model of a patch of neuronal membrane containing stochastic Na(+) and K(+) channels obeying Hodgkin-Huxley kinetics. The quality of signal encoding is assessed by reconstructing the input stimuli from the output spike trains using optimal linear mean square estimation. A comparison of the estimation performance of noisy neuronal models of spike generation enables us to assess the impact of neuronal noise on the efficacy of neural coding. The results for both models suggest that spike timing variability reduces the ability of spike trains to encode rapid time-varying stimuli. Moreover, contrary to expectations based on earlier studies, we find that the noisy spike encoding models encode slowly varying stimuli more effectively than rapidly varying ones.     

The correlation grid: analysis of synchronous spiking in multi-dimensional spike train data and identification of feasible connection architectures

This paper presents a visualization technique specifically designed to support the analysis of synchronous firings in multiple, simultaneously recorded, spike trains. This technique, called the correlation grid, enables investigators to identify groups of spike trains, where each pair of spike trains has a high probability of generating spikes approximately simultaneously or within a constant time shift. Moreover, the correlation grid was developed to help solve the following reverse problem: identification of the connection architecture between spike train generating units, which may produce a spike train dataset similar to the one under analysis. To demonstrate the efficacy of this approach, results are presented from a study of three simulated, noisy, spike train datasets. The parameters of the simulated neurons were chosen to reflect the typical characteristics of cortical pyramidal neurons. The schemes of neuronal connections were not known to the analysts. Nevertheless, the correlation grid enabled the analysts to find the correct connection architecture for each of these three data sets.     

Transmission of spike trains by nerve fibers after near-threshold stimulation

Spike trains in group A nerve fibers were studied in anesthetized cats in response to stimulation of the hind-limb nerves by random (Poisson) and regular sequences of stimuli. In response to above-threshold stimulation spike trains in the nerve fibers were shown to differ from the stimulating trains purely in the absence of intervals less than 1–1.5 msec in duration, as a result of the presence of a refractory period. With near-threshold stimulation with an average frequency of over 10 per second, spike trains differed significantly from the stimulating trains, as reflected in histograms of interspike intervals, the shape of the intensity function, and the magnitude of the coefficient of correlation for successive intervals. It is postulated that changes in the structure of the spike trains conveyed by a nerve fiber are attributable to the presence of after-activity.A. A. Bogomolets Institute of Physiology, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Neirofiziologiya, Vol. 8, No. 1, pp. 91–98, January–February, 1976.     

A cross-interval spike train analysis: the correlation between spike generation and temporal integration of doublets

A stochastic spike train analysis technique is introduced to reveal the correlation between the firing of the next spike and the temporal integration period of two consecutive spikes (i.e., a doublet). Statistics of spike firing times between neurons are established to obtain the conditional probability of spike firing in relation to the integration period. The existence of a temporal integration period is deduced from the time interval between two consecutive spikes fired in a reference neuron as a precondition to the generation of the next spike in a compared neuron. This analysis can show whether the coupled spike firing in the compared neuron is correlated with the last or the second-to-last spike in the reference neuron. Analysis of simulated and experimentally recorded biological spike trains shows that the effects of excitatory and inhibitory temporal integration are extracted by this method without relying on any subthreshold potential recordings. The analysis also shows that, with temporal integration, a neuron driven by random firing patterns can produce fairly regular firing patterns under appropriate conditions. This regularity in firing can be enhanced by temporal integration of spikes in a chain of polysynaptically connected neurons. The bandpass filtering of spike firings by temporal integration is discussed. The results also reveal that signal transmission delays may be attributed not just to conduction and synaptic delays, but also to the delay time needed for temporal integration. Received: 3 March 1997 / Accepted in revised form: 6 November 1997     

Use of `relative-phase' analysis to assess correlation between neuronal spike trains

 We propose a new method of studying the correlation between neuronal spike trains. This technique is based on the analysis of relative phase between two point processes. Relative phase here is defined as the relative timing difference between two spike trains normalized by the associated interspike interval of one cell. This phase measurement is intended to reveal the relative timing relationship between spike trains atdifferent firing rates. We apply this method to a numerical example and an example from two cerebellar neuronal spike trains of a behaving rat. The results are compared with classical cross-correlation analysis. We show that the technique can avoid some of the limitations of cross-correlation methods, reveal certain statistical dependencies that cannot be shown by cross correlation, and provide information as to the direction of influence between two spike trains. Received: 8 November 2001 / Accepted: 30 September 2002 / Published online: 24 January 2003 Correspondence to: Y. Chen (e-mail: chen@nsi.edu, Fax: + 1-858-626-2099) Acknowledgements. Research for this paper was supported by the Alafi Family Foundation and the Neurosciences Research Foundation.     

An information-geometric framework for statistical inferences in the neural spike train space

Statistical inferences are essentially important in analyzing neural spike trains in computational neuroscience. Current approaches have followed a general inference paradigm where a parametric probability model is often used to characterize the temporal evolution of the underlying stochastic processes. To directly capture the overall variability and distribution in the space of the spike trains, we focus on a data-driven approach where statistics are defined and computed in the function space in which spike trains are viewed as individual points. To this end, we at first develop a parametrized family of metrics that takes into account different warpings in the time domain and generalizes several currently used spike train distances. These new metrics are essentially penalized L ^p norms, involving appropriate functions of spike trains, with penalties associated with time-warping. The notions of means and variances of spike trains are then defined based on the new metrics when p = 2 (corresponding to the “Euclidean distance”). Using some restrictive conditions, we present an efficient recursive algorithm, termed Matching-Minimization algorithm, to compute the sample mean of a set of spike trains with arbitrary numbers of spikes. The proposed metrics as well as the mean spike trains are demonstrated using simulations as well as an experimental recording from the motor cortex. It is found that all these methods achieve desirable performance and the results support the success of this novel framework.