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.     

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.     

Regression analysis of non-stationary discharges in neurons; Adaptation in the electrosensory afferent of dogfish

Interspike interval histograms, as usually regarded for the estimation of statistical variabilities in neuronal spike trains, were applied to non-stationary dynamic responses of a PD receptor. Sliding mean values were introduced describing the average receptor response on defined, recurrent stimuli; mean spike frequencies and interspike intervals were computed a) for fixed sequential analysis periods (of e.g. 500 ms), b) for analysis periods shifted by every consecutive interspike interval (thus the number of spikes being constant), and c) by fitting the dynamic responses for suitable analytic functions (e.g. exponential functions). With these methods variabilities in the non-stationary neuronal impulse patterns were investigated for electrosensory PD afferents in Lorenzinian ampulla of dogfish (Scyliorhinus canicula) with electric stimuli up to 50 nA and defined temperatures between 7° C and 25° C. In this temperature range all investigated ampullae were spontaneously active, the irregularities in neuronal discharges and averaged spike frequencies depended strongly on temperature, the latter showing maxima between 13° C and 19° C. In preparations with small disturbances we generally found static interspike interval histograms following approximatively a Gaussian distribution. The same was true for the momentary spike frequency and its deviation during the dynamic response to given electrical stimuli. A suprathreshold rectangular current (e.g.-0.5 nA) led to a marked but transient synchronisation in spike generation; the higher the stimulus strength, the smaller the standard deviation (s.d.) from mean spike frequency in the beginning of the dynamic response; during adaptation the s.d. increased up to that of the static response frequency. Relating, however, s.d. for different currents, times, and temperatures to the corresponding mean spike frequency led to fairly constant coefficients of variation; s.d. was approximatively a linear function of the sliding mean value even in the dynamic response of the electroreceptor (scaling).Supported by the Deutsche Forschungsgemeinschaft (Br 310/11)     

Interspike interval dependency from arterial chemoreceptors

The carotid body impulse generator has been previously characterized as a Poisson-type random process. We examined the validity of this characterization by analyzing sinus nerve spike trains for interspike interval dependency. Fifteen single chemoreceptive afferents were recorded in vivo under hypoxic-hypercapnic conditions, and approximately 1,000 consecutive interspike intervals for each fiber were timed and analyzed for serial dependence. The same set of intervals placed in shuffled order served as a control series without serial dependence. The original spike interval trains showed significantly negative first-order serial correlation coefficients and less variability in joint interval distributions than did the shuffled interval trains. These results suggest that the chemoreceptor afferent train is not random and may reflect a negative feedback system operating within the carotid body that limits variation about a mean frequency.     

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.     

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.     

体感皮层神经元放电间隔的概率密度函数与分布参数

本文建立了估计ISI概率密度函数的标准化ISI直方图和分布参数拟合方法,对34例猫体感皮层神经元自发和诱发放电活动进行了统计分析.     

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.     

Near Poisson-type firing produced by concurrent excitation and inhibition

The effect of inhibition on the firing variability is examined in this paper using the biologically-inspired temporal noisy-leaky integrator (TNLI) neuron model. The TNLI incorporates hyperpolarising inhibition with negative current pulses of controlled shapes and it also separates dendritic from somatic integration. The firing variability is observed by looking at the coefficient of variation (C(V)) (standard deviation/mean interspike interval) as a function of the mean interspike interval of firing (delta tM) and by comparing the results with the theoretical curve for random spike trains, as well as looking at the interspike interval (ISI) histogram distributions. The results show that with 80% inhibition, firing at high rates (up to 200 Hz) is nearly consistent with a Poisson-type variability, which complies with the analysis of cortical neuron firing recordings by Softky and Koch [1993, J. Neurosci. 13(1) 334-530]. We also demonstrate that the mechanism by which inhibition increases the C(V) values is by introducing more short intervals in the firing pattern as indicated by a small initial hump at the beginning of the ISI histogram distribution. The use of stochastic inputs and the separation of the dendritic and somatic integration which we model in TNLI, also affect the high firing, near Poisson-type (explained in the paper) variability produced. We have also found that partial dendritic reset increases slightly the firing variability especially at short ISIs.     

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     

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.     

Estimation of the Optimal Statistical Quality Control Sampling Time Intervals Using a Residual Risk Measure

Background

An open problem in clinical chemistry is the estimation of the optimal sampling time intervals for the application of statistical quality control (QC) procedures that are based on the measurement of control materials. This is a probabilistic risk assessment problem that requires reliability analysis of the analytical system, and the estimation of the risk caused by the measurement error.

Methodology/Principal Findings

Assuming that the states of the analytical system are the reliability state, the maintenance state, the critical-failure modes and their combinations, we can define risk functions based on the mean time of the states, their measurement error and the medically acceptable measurement error. Consequently, a residual risk measure rr can be defined for each sampling time interval. The rr depends on the state probability vectors of the analytical system, the state transition probability matrices before and after each application of the QC procedure and the state mean time matrices. As optimal sampling time intervals can be defined those minimizing a QC related cost measure while the rr is acceptable. I developed an algorithm that estimates the rr for any QC sampling time interval of a QC procedure applied to analytical systems with an arbitrary number of critical-failure modes, assuming any failure time and measurement error probability density function for each mode. Furthermore, given the acceptable rr, it can estimate the optimal QC sampling time intervals.

Conclusions/Significance

It is possible to rationally estimate the optimal QC sampling time intervals of an analytical system to sustain an acceptable residual risk with the minimum QC related cost. For the optimization the reliability analysis of the analytical system and the risk analysis of the measurement error are needed.     

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.     

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.     

Temporal pattern sensitive and nonsensitive responses in the cat's retinal ganglion cells

In order to characterize temporal pattern sensitivity in the cat ganglion cells, a new analysis technique by semi-Markov models which was developed in the previous papers (Tsukada et al., 1975–1977) was applied to input-output relations of the receptive-field. Three types of statistical spot stimuli positioned in the center region of receptive fields were used. Each type of stimulus has an identical histogram in the inter-stimulus intervals and therefore the same mean and variance, but different correlations between adjacent inter-stimulus intervals (Type 1, positive; Type 2, negative; and Type 3, independent processes). From the output spike trains of cat retinal ganglion cells to each stimulus, mean, variance, and histogram were computed. As the result of investigating these data, we could draw the following conclusion from the resultant output interval histograms. The receptive-field-center responses of cat ganglion cells can be classified into two groups (Types L and N) according to the difference of responsiveness to the three types of statistical spot stimuli. A Type L response has the same histogram in interspike intervals for all three stimuli, and is not sensitive to the temporal pattern, while a Type N response has three different forms depending on each type of stimulus showing high sensitivity to the temporal pattern. These results were also simulated by the Markov chain model and discussed with relation to neural coding and classification of ganglion cell types.     

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.     

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.     

A Superposition Model of the Spontaneous Activity of Cerebellar Purkinje Cells

Based on physiological evidence for multiple firing zones in the dendritic arborizations of cerebellar Purkinje cells, a superposition model is proposed for spike triggering in these cells. Spike trains from 10 Purkinje cells were analyzed in terms of independence of interspike intervals and the properties of their variance-time curves. The results of this analysis were found consistent with the hypothesis that the spike train of a cerebellar Purkinje cell is the pooled output of a relatively large number of independent component processes. Simplifying assumptions as to the statistical nature of these processes lead to a very rough estimate of the number of firing zones.