Analysis¶
Statistics of spike trains¶
-
brian.
firing_rate
(spikes)¶ Rate of the spike train.
-
brian.
CV
(spikes)¶ Coefficient of variation.
-
brian.
correlogram
(T1, T2, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns a cross-correlogram with lag in [-width,width] and given bin size. T is the total duration (optional) and should be greater than the duration of T1 and T2. The result is in Hz (rate of coincidences in each bin).
N.B.: units are discarded. TODO: optimise?
-
brian.
autocorrelogram
(T0, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns an autocorrelogram with lag in [-width,width] and given bin size. T is the total duration (optional) and should be greater than the duration of T1 and T2. The result is in Hz (rate of coincidences in each bin).
N.B.: units are discarded.
-
brian.
CCF
(T1, T2, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns the cross-correlation function with lag in [-width,width] and given bin size. T is the total duration (optional). The result is in Hz**2: CCF(T1,T2)=<T1(t)T2(t+s)>
N.B.: units are discarded.
-
brian.
ACF
(T0, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns the autocorrelation function with lag in [-width,width] and given bin size. T is the total duration (optional). The result is in Hz**2: ACF(T0)=<T0(t)T0(t+s)>
N.B.: units are discarded.
-
brian.
CCVF
(T1, T2, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns the cross-covariance function with lag in [-width,width] and given bin size. T is the total duration (optional). The result is in Hz**2: CCVF(T1,T2)=<T1(t)T2(t+s)>-<T1><T2>
N.B.: units are discarded.
-
brian.
ACVF
(T0, width=20.0 * msecond, bin=1.0 * msecond, T=None)¶ Returns the autocovariance function with lag in [-width,width] and given bin size. T is the total duration (optional). The result is in Hz**2: ACVF(T0)=<T0(t)T0(t+s)>-<T0>**2
N.B.: units are discarded.
-
brian.
total_correlation
(T1, T2, width=20.0 * msecond, T=None)¶ Returns the total correlation coefficient with lag in [-width,width]. T is the total duration (optional). The result is a real (typically in [0,1]): total_correlation(T1,T2)=int(CCVF(T1,T2))/rate(T1)
-
brian.
spike_triggered_average
(spikes, stimulus, max_interval, dt, onset=None, display=False)¶ Spike triggered average reverse correlation. spikes is an array containing spike times stimulus is an array containing the stimulus max_interval (second) is the duration of the averaging window dt (second) is the sampling period onset (second) before which the spikes are discarded. Note: it will be at least as long as max_interval display (default=False) display the number of spikes processed out of the total number output the spike triggered average and the corresponding time axis