Quantifying Spike Train Oscillations: Biases,Distortions and Solutions

Estimation of the power spectrum is a common method for identifying oscillatory changes in neuronal activity. However, the stochastic nature of neuronal activity leads to severe biases in the estimation of these oscillations in single unit spike trains. Different biological and experimental factors cause the spike train to differentially reflect its underlying oscillatory rate function. We analyzed the effect of factors, such as the mean firing rate and the recording duration, on the detectability of oscillations and their significance, and tested these theoretical results on experimental data recorded in Parkinsonian non-human primates. The effect of these factors is dramatic, such that in some conditions, the detection of existing oscillations is impossible. Moreover, these biases impede the comparison of oscillations across brain regions, neuronal types, behavioral states and separate recordings with different underlying parameters, and lead inevitably to a gross misinterpretation of experimental results. We introduce a novel objective measure, the "modulation index", which overcomes these biases, and enables reliable detection of oscillations from spike trains and a direct estimation of the oscillation magnitude. The modulation index detects a high percentage of oscillations over a wide range of parameters, compared to classical spectral analysis methods, and enables an unbiased comparison between spike trains recorded from different neurons and using different experimental protocols.     

Note on the estimation of the correlation function of neural spike trains

The use of time-bins in the estimation of the correlation function of neural spike trains has a filtering effect on the estimate and results in distortion and aliasing. Prior low-pass filtering of the spike trains, on the other hand, and computation of the correlation function of the emerging waveforms in the standard way result in an estimate that is also a filtered version of the original function but distortion- and alias-free. In addition, the correlation function so computed can be normalized. An analogous definition of the correlation coefficient for the first technique enables the comparison of these various correlation estimates and clarifies their properties.     

Neuronal spike trains and stochastic point processes. I. The single spike train

In a growing class of neurophysiological experiments, the train of impulses (“spikes”) produced by a nerve cell is subjected to statistical treatment involving the time intervals between spikes. The statistical techniques available for the analysis of single spike trains are described and related to the underlying mathematical theory, that of stochastic point processes, i.e., of stochastic processes whose realizations may be described as series of point events occurring in time, separated by random intervals. For single stationary spike trains, several orders of complexity of statistical treatment are described; the major distinction is that between statistical measures that depend in an essential way on the serial order of interspike intervals and those that are order-independent. The interrelations among the several types of calculations are shown, and an attempt is made to ameliorate the current nomenclatural confusion in this field. Applications, interpretations, and potential difficulties of the statistical techniques are discussed, with special reference to types of spike trains encountered experimentally. Next, the related types of analysis are described for experiments which involve repeated presentations of a brief, isolated stimulus. Finally, the effects of nonstationarity, e.g. long-term changes in firing rate, on the various statistical measures are discussed. Several commonly observed patterns of spike activity are shown to be differentially sensitive to such changes. A companion paper covers the analysis of simultaneously observed spike trains.     

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.     

Reconstruction of Underlying Nonlinear Deterministic Dynamics Embedded in Noisy Spike Trains

An experimentally recorded time series formed by the exact times of occurrence of the neuronal spikes (spike train) is likely to be affected by observational noise that provokes events mistakenly confused with neuronal discharges, as well as missed detection of genuine neuronal discharges. The points of the spike train may also suffer a slight jitter in time due to stochastic processes in synaptic transmission and to delays in the detecting devices. This study presents a procedure aimed at filtering the embedded noise (denoising the spike trains) the spike trains based on the hypothesis that recurrent temporal patterns of spikes are likely to represent the robust expression of a dynamic process associated with the information carried by the spike train. The rationale of this approach is tested on simulated spike trains generated by several nonlinear deterministic dynamical systems with embedded observational noise. The application of the pattern grouping algorithm (PGA) to the noisy time series allows us to extract a set of points that form the reconstructed time series. Three new indices are defined for assessment of the performance of the denoising procedure. The results show that this procedure may indeed retrieve the most relevant temporal features of the original dynamics. Moreover, we observe that additional spurious events affect the performance to a larger extent than the missing of original points. Thus, a strict criterion for the detection of spikes under experimental conditions, thus reducing the number of spurious spikes, may raise the possibility to apply PGA to detect endogenous deterministic dynamics in the spike train otherwise masked by the observational noise.     

自发和诱发神经放电序列的内在时序特性

运用Fano因子分析法,考察豚鼠听神经单纤维的自发放电序列、小鼠海马CAl区神经元的自发放电序列以及蟾蜍缝匠肌肌梭传入神经的诱发放电序列的时序特性,结果显示自发和诱发放电时间序列均存在Fano因子随计算窗口时间的增大而持续增长的特点,而原始数据的随机重排替代数据则没有这一特性,说明这些神经放电时间序列与一般的随机点过程不同,存在长时程相关性,在时序上具有某种结构特征。进一步的研究表明,这一时序结构特征可以通过将随机产生的一维正态分布序列数据,与神经放电时间序列数据进行跟随排序后而体现,提示这一特征与放电间隔的分布特点无关。     

STOCHASTIC RESONANCE MEASURED BY THEMUTUAL INFORMATION IN THE NEURONSTHAT TRANSMIT CHAOTIC SPIKE TRAINS

The spike trains generated by a neuron model are studied by the methods of nonlinear time series analysis. The results show that the spike trains are chaotic. To investigate effect of noise on transmission of chaotic spike trains, this chaotic spike trains are used as a discrete subthreshold input signal to the integrate-and-fire neuronal model and the FitzHugh-Nagumo(FHN) neuronal model working in noisy environment. The mutual information between the input spike trains and the output spike trains is calculated, the result shows that the transformation of information encoded by the chaotic spike trains is optimized by some level of noise, and stochastic resonance(SR) measured by mutual information is a property available for neurons to transmit chaotic spike trains.     

STOCHASTIC RESONANCE MEASURED BY THEMUTUAL INFORMATION IN THE NEURONSTHAT TRANSMIT CHAOTIC SPIKE TRAINS

１ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌｋｎｏｗｎｔｈａｔｎｅｒｖｅｃｅｌｌｓｗｏｒｋｉｎｎｏｉｓｙｅｎｖｉｒｏｎｍｅｎｔ,ａｎｄｎｏｉｓｅｓｏｕｒｃｅｓｒａｎｇｉｎｇｆｒｏｍｉｎｔｅｒｎａｌｔｈｅｒｍａｌｎｏｉｓｅｔｏｅｘｔｅｒｎａｌｐｅｒｔｕｒｂａｔｉｏｎ．Ｏｎｅｐｕｚｚｌｉｎｇｐｒｏｂｌｅｍｉｓｈｏｗｄｏｎｅｒｖｅｃｅｌｌｓａｃｃｏｍｍｏｄａｔｅｎｏｉｓｅｉｎｃｏｄｉｎｇａｎｄｔｒａｎｓｆｏｒｍｉｎｇｉｎｆｏｒｍａｔｉｏｎ,ｒｅｃｅｎｔｒｅｓｅａｒｃｈｓｈｏｗｓｔｈａｔｎｏｉｓｅｍａｙｐ…     

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.     

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.     

Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains

The statistical analysis of two simultaneously observed trains of neuronal spikes is described, using as a conceptual framework the theory of stochastic point processes.The first statistical question that arises is whether the observed trains are independent; statistical techniques for testing independence are developed around the notion that, under the null hypothesis, the times of spike occurrence in one train represent random instants in time with respect to the other. If the null hypothesis is rejected—if dependence is attributed to the trains—the problem then becomes that of characterizing the nature and source of the observed dependencies. Statistical signs of various classes of dependencies, including direct interaction and shared input, are discussed and illustrated through computer simulations of interacting neurons. The effects of nonstationarities on the statistical measures for simultaneous spike trains are also discussed. For two-train comparisons of irregularly discharging nerve cells, moderate nonstationarities are shown to have little effect on the detection of interactions.Combining repetitive stimulation and simultaneous recording of spike trains from two (or more) neurons yields additional clues as to possible modes of interaction among the monitored neurons; the theory presented is illustrated by an application to experimentally obtained data from auditory neurons.A companion paper covers the analysis of single spike trains.     

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.     

A joint interspike interval difference stochastic spike train analysis: detecting local trends in the temporal firing patterns of single neurons

We introduce a stochastic spike train analysis method called joint interspike interval difference (JISID) analysis. By design, this method detects changes in firing interspike intervals (ISIs), called local trends, within a 4-spike pattern in a spike train. This analysis classifies 4-spike patterns that have similar incremental changes. It characterizes the higher-order serial dependence in spike firing relative to changes in the firing history. Mathematically, this spike train analysis describes the statistical joint distribution of consecutive changes in ISIs, from which the serial dependence of the changes in higher-order intervals can be determined. It is similar to the joint interspike interval (JISI) analysis, except that the joint distribution of consecutive ISI differences (ISIDs) is quantified. The graphical location of points in the JISID scatter plot reveals the local trends in firing (i.e., monotonically increasing, monotonically decreasing, or transitional firing). The trajectory of these points in the serial-JISID plot traces the time evolution of these trends represented by a 5-spike pattern, while points in the JISID scatter plot represent trends of a 4-spike pattern. We provide complete theoretical interpretations of the JISID analysis. We also demonstrate that this method indeed identifies firing trends in both simulated spike trains and spike trains recorded from cultured neurons. Received: 13 May 1997 / Accepted in revised form: 9 December 1998     

Statistical Analysis of Motor Unit Firing Patterns in a Human Skeletal Muscle

A statistical analysis of the firing pattern of single motor units in the human brachial biceps muscle is presented. Single motor unit spike trains are recorded and analyzed. The statistical treatment of these spike trains is as stochastic point processes, the theory of which is briefly discussed. Evidence is presented that motor unit spike trains may be modelled by a renewal process with an underlying gaussian probability density. Statistical independence of successive interspike intervals is shown using scatter diagrams; the hypothesis of a gaussian distribution is accepted at the 99th percentile confidence limit, chi-square test, in 90% of the units tested. A functional relationship between the mean and standard deviation is shown and discussed; its implications in obtaining sample size are presented in an appendix.The results of higher order analysis in the form of autocorrelograms and grouped interval histograms are presented. Grouped interval histograms are discussed in the context of motor unit data, and used to confirm the hypothesis that a stable probability density function does not represent a good model of the data at this level of analysis.     

Identification of connectivity in neural networks.

Analytical and experimental methods are provided for estimating synaptic connectivities from simultaneous recordings of multiple neurons. The results are based on detailed, yet flexible neuron models in which spike trains are modeled as general doubly stochastic point processes. The expressions derived can be used with nonstationary or stationary records, and can be readily extended from pairwise to multineuron estimates. Furthermore, we show analytically how the estimates are improved as more neurons are sampled, and derive the appropriate normalizations to eliminate stimulus-related correlations. Finally, we illustrate the use and interpretation of the analytical expressions on simulated spike trains and neural networks, and give explicit confidence measures on the estimates.     

神经放电序列模式的一种识别方法及应用

提出了神经放电序列模式识别的一种新方法。首先,把放电序列用阶梯状的响应函数来表示,然后定义了其一阶、二阶形式导数以及形式积分。这三个特征量均有着不同的几何和物理意义,因此采用这三个特征量来刻画神经放电序列的模式,就可以较全面地表示其特征。对神经放电序列的重构也表明通过这几个特征量可以很好地反映序列中所包含的信息。作为应用例子,这种量化方法用来研究冷热感受器模型所产生的放电模式,结果表明它能够识别在不同温度条件下的放电模式。     

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     

Model based decoding of spike trains

Reliably decoding neuronal responses requires knowing what aspects of neuronal responses are stimulus related, and which aspects act as noise. Recent work shows that spike trains can be viewed as stochastic samples from the rate variation function, as estimated by the time dependent spike density function (or normalized peristimulus time histogram). Such spike trains are exactly described by order statistics, and can be decoded millisecond-by-millisecond by iterative application of order statistics.     

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.     

Fourier analysis of spike train data

A formal representation of nerve spike trains in the form of a sum of rectangular functions is proposed. This formal instantaneous frequency function can be Fourier analyzed. The resulting algorithm has the useful properties of spike by spike calculations and an insensitivity to the mean (carrier) spike rate. The technique is also useful for producing a smooth (filtered) reconstruction of a spike train.