Parametric inference of neuronal response latency in presence of a background signal

Presynaptic Spontaneous Activity Enhances the Accuracy of Latency Coding

The time to the first spike after stimulus onset typically varies with the stimulation intensity. Experimental evidence suggests that neural systems use such response latency to encode information about the stimulus. We investigate the decoding accuracy of the latency code in relation to the level of noise in the form of presynaptic spontaneous activity. Paradoxically, the optimal performance is achieved at a nonzero level of noise and suprathreshold stimulus intensities. We argue that this phenomenon results from the influence of the spontaneous activity on the stabilization of the membrane potential in the absence of stimulation. The reported decoding accuracy improvement represents a novel manifestation of the noise-aided signal enhancement.

Effect of spontaneous activity on stimulus detection in a simple neuronal model

It is studied what level of a continuous-valued signal is optimally estimable on the basis of first-spike latency neuronal data. When a spontaneous neuronal activity is present, the first spike after the stimulus onset may be caused either by the stimulus itself, or it may be a result of the prevailing spontaneous activity. Under certain regularity conditions, Fisher information is the inverse of the variance of the best estimator. It can be considered as a function of the signal intensity and then indicates accuracy of the estimation for each signal level. The Fisher information is normalized with respect to the time needed to obtain an observation. The accuracy of signal level estimation is investigated in basic discharge patterns modelled by a Poisson and a renewal process and the impact of the complex interaction between spontaneous activity and a delay of the response is shown.

Parameter inference from hitting times for perturbed Brownian motion

A latent internal process describes the state of some system, e.g. the social tension in a political conflict, the strength of an industrial component or the health status of a person. When this process reaches a predefined threshold, the process terminates and an observable event occurs, e.g. the political conflict finishes, the industrial component breaks down or the person dies. Imagine an intervention, e.g., a political decision, maintenance of a component or a medical treatment, is initiated to the process before the event occurs. How can we evaluate whether the intervention had an effect? To answer this question we describe the effect of the intervention through parameter changes of the law governing the internal process. Then, the time interval between the start of the process and the final event is divided into two subintervals: the time from the start to the instant of intervention, denoted by S, and the time between the intervention and the threshold crossing, denoted by R. The first question studied here is: What is the joint distribution of (S,R)? The theoretical expressions are provided and serve as a basis to answer the main question: Can we estimate the parameters of the model from observations of S and R and compare them statistically? Maximum likelihood estimators are calculated and applied on simulated data under the assumption that the process before and after the intervention is described by the same type of model, i.e. a Brownian motion, but with different parameters. Also covariates and handling of censored observations are incorporated into the statistical model, and the method is illustrated on lung cancer data.

A review of the methods for neuronal response latency estimation

Neuronal response latency is usually vaguely defined as the delay between the stimulus onset and the beginning of the response. It contains important information for the understanding of the temporal code. For this reason, the detection of the response latency has been extensively studied in the last twenty years, yielding different estimation methods. They can be divided into two classes, one of them including methods based on detecting an intensity change in the firing rate profile after the stimulus onset and the other containing methods based on detection of spikes evoked by the stimulation using interspike intervals and spike times. The aim of this paper is to present a review of the main techniques proposed in both classes, highlighting their advantages and shortcomings.

On Firing Rate Estimation for Dependent Interspike Intervals

If interspike intervals are dependent, the instantaneous firing rate does not catch important features of spike trains. In this case, the conditional instantaneous rate plays the role of the instantaneous firing rate for the case of samples of independent interspike intervals. If the conditional distribution of the interspikes intervals (ISIs) is unknown, it becomes difficult to evaluate the conditional firing rate. We propose a nonparametric estimator for the conditional instantaneous firing rate for Markov, stationary, and ergodic ISIs. An algorithm to check the reliability of the proposed estimator is introduced, and its consistency properties are proved. The method is applied to data obtained from a stochastic two-compartment model and to in vitro experimental data.

Estimating latency from inhibitory input

Stimulus response latency is the time period between the presentation of a stimulus and the occurrence of a change in the neural firing evoked by the stimulation. The response latency has been explored and estimation methods proposed mostly for excitatory stimuli, which means that the neuron reacts to the stimulus by an increase in the firing rate. We focus on the estimation of the response latency in the case of inhibitory stimuli. Models used in this paper represent two different descriptions of response latency. We consider either the latency to be constant across trials or to be a random variable. In the case of random latency, special attention is given to models with selective interaction. The aim is to propose methods for estimation of the latency or the parameters of its distribution. Parameters are estimated by four different methods: method of moments, maximum-likelihood method, a method comparing an empirical and a theoretical cumulative distribution function and a method based on the Laplace transform of a probability density function. All four methods are applied on simulated data and compared.