In the next few sections, we introduce some important concepts commonly used for the statistical description of neuronal spike trains. Central notions will be the interspike interval distribution, (Section 7.3), the noise spectrum (Section 7.4), but most importantly the concept of ‘firing rate’ which we discuss first.
A quick glance at the experimental literature reveals that there is no unique and well-defined concept of ‘mean firing rate’. In fact, there are at least three different notions of rate which are often confused and used simultaneously. The three definitions refer to three different averaging procedures: either an average over time, or an average over several repetitions of the experiment, or an average over a population of neurons. The following three subsections will revisit in detail these three concepts.
The first and most commonly used definition of a firing rate refers to a temporal average. An experimentalist observes in trial the spikes of a given neuron (see Fig. 7.6). The firing rate in trial is the spike count in an interval of duration divided by .
The length of the time window is set by the experimenter and depends on the type of neuron and the stimulus. In practice, to get sensible averages, several spikes should occur within the time window. Typical values are ms or ms, but the duration may also be longer or shorter.
This definition of rate has been successfully used in many preparations, particularly in experiments on sensory or motor systems. A classical example is the stretch receptor in a muscle spindle (9). The number of spikes emitted by the receptor neuron increases with the force applied to the muscle.
If the same experiment is repeated several times, the measured spike count varies between one trial and the next. Let us denote the spike count in trial by the variable , its mean by and deviations from the mean as . Variability of the spike count measure is characterized by the Fano Factor, defined as the variance of the spike count divided by its mean
In experiments, the mean and variance are estimated by averaging over trials and .
If we find on average spikes in a long temporal window of duration , the mean interval between two subsequent spikes is . Indeed, using the notion of interspike-interval distribution to be introduced below (Section 7.3), we can make the following statement: the firing rate defined here as spike count divided by the measurement time is identical to the inverse of the mean interspike interval. We will come back to interspike intervals in Section 7.3.
It is tempting, but misleading, to consider the inverse interspike interval as a ’momentary firing rate’: if a first spike occurs at time and the next one at time , we could artificially assign a variable for all times . However, the temporal average of over a much longer time is not the same as a the mean rate defined here as spike count divided by , simply because . A practical definition of ’instantaneous firing rate’ will be given below in Section 7.2.2.
An experimenter records from a neuron while stimulating with some input sequence. The same stimulation sequence is repeated several times and the neuronal response is reported in a Peri-Stimulus-Time Histogram (PSTH) with bin width ; see Fig. 7.7. The time is measured with respect to the start of the stimulation sequence and defines the time bin for generating the histogram, it is typically on the order of milliseconds.
The number of occurrences of spikes summed over all repetitions of the experiment divided by the number of repetitions is a measure of the typical activity of the neuron between time and . A further division by the interval length yields the spike density
Sometimes the result is smoothed to get a continuous (time-dependent) rate variable, usually reported in units of Hz. As an experimental procedure, the PSTH measure is a useful method to evaluate neuronal activity, in particular in the case of time-dependent stimuli; see Fig. 7.4. We call it the time-dependent firing rate.
In order to see the relation of Eq. 7.6 to a time-dependent firing rate, we recall that spikes are formal events characterized by their firing time where counts the spikes. In Chapter 1 we have defined (Eq. (1.14)) the spike train as a sum of -functions:
If each stimulation can be considered as an independent sample from the identical stochastic process, we can define an instantaneous firing rate as an expectation over trials
An expectation value over -functions may look weird to the reader not used to seeing such mathematical objects. Let us therefore consider the experimental procedure to estimate the expectation value. First, in each trial , we count the number of spikes that occur in a short time interval by integrating the spike train over time, where the lower index denotes the trial number. Note that integration removes the -function. Obviously, if the time bin is small enough we will find at most one spike so that is either zero or one. Second, we average over the trials and divide by in order to have the empirical estimate
The PSTH, defined as spike count per time bin averaged over several trials and divided by the bin length (the right-hand side of Eq. (7.9)), provides therefore an empirical estimate of the instantaneous firing rate (the left-hand side).
The number of neurons in the brain is huge. Often many neurons have similar properties and respond to the same stimuli. For example, neurons in the primary visual cortex of cats and monkeys are arranged in columns of cells with similar properties (230). Let us idealize the situation and consider a population of neurons with identical properties. In particular, all neurons in the population should have the same pattern of input and output connections. The spikes of the neurons in a population are sent off to another population . In our idealized picture, each neuron in population receives input from all neurons in population . The relevant quantity, from the point of view of the receiving neuron, is the proportion of active neurons in the presynaptic population ; see Fig. 7.8A. Formally, we define the population activity
where is the size of the population, is the number of spikes (summed over all neurons in the population) that occur between and where is a small time interval; see Fig. 7.8. Eq. (7.11) defines a variable with units inverse time – in other words, a rate.
© Cambridge University Press. This book is in copyright. No reproduction of any part of it may take place without the written permission of Cambridge University Press.