SPT–Data Analysis and Theoretical Routines

These M-files implement the analysis procedures discussed in chapter 9 of “Methods in Neuronal Modeling”. In addition, many theoretical results described there are also implemented in the form of M-files, allowing a direct comparison between simulations and theory. References to the literature are provided in the M-files themselves or through the Matlab help utility. This includes the following routines:

Data Analysis

countprob: returns the probability density of the spike count in a fixed time interval.
dispersion: returns the variance in the spike count as a function of the mean spike count for a range of time intervals.
isidist: computes the interspike interval distribution and returns its mean and coefficient of variation.
psautostim: estimates the power spectrum and autocorrelation function of a random input stimulus.
psautospk: estimates the power spectrum and autocorrelation function of a spike train.
roc: computes from two spike count probability distributions the ROC curve characterizing the discrimination performance of an ideal observer based on these distributions.
stimest: estimates a random stimulus from a spike train by Wiener-Kolmogorov filtering.
wiener1: estimates the transfer function characterizing the transformation between a random stimulus and output spike train.

Theoretical Routines

buttfilt*: plots the transfer function of Butterworth low-pass filter.
fithifref: returns the F-I curve of a perfect integrate-and-fire neuron with refractory period.
fithlifref: returns the F-I curve of a leaky integrate-and-fire neuron with refractory period.
gaussroc: returns the ROC curve characterizing the discrimination performance of an ideal observer between two Gaussian probability distributions.
gfitheor: returns the ISI distribution of an integrate-and-fire neuron which sums n_th Poisson inputs to reach threshold (this also corresponds to the ISI distribution of an integrate-and-fire neuron with gamma distributed random thereshold of order n_th).
gpstheor: returns the power spectrum and autocorrelation function of an integrate-and-fire neuron which sums n_th Poisson inputs to reach threshold.
lgnfilt*: returns the transfer function of an LGN relay cell. neymtacnt: returns the probability density of a Neyman type A spike count distribution.
poisscnt: returns the probability density of a Poisson spike count distribution.
prefdisp: returns the variance in spike count as a function of mean spike count for a Poisson neuron with refractory period.
sndpdisp: returns the approximate value of the variance in spike count as function of the mean spike count for a doubly stochastic Poisson process where the driving rate is a shot noise with a square impulse response.
stesttheor: returns the fraction of a white Gaussian signal encoded in two half-wave rectifying Poisson neurons.

The M-files marked with an asterisk (*) are mainly convienent to set up simulation models.