elephant.spike_train_generation.single_interaction_process

elephant.spike_train_generation.single_interaction_process(rate, coincidence_rate, t_stop, n_spiketrains=2, jitter=array(0.) * ms, coincidences='deterministic', t_start=array(0.) * ms, min_delay=array(0.) * ms, return_coincidences=False)[source]

Generates a multidimensional Poisson SIP (single interaction process) plus independent Poisson processes (Kuhn et al., 2003).

A Poisson SIP consists of Poisson time series which are independent except for simultaneous events in all of them. This routine generates a SIP plus additional parallel independent Poisson processes.

Parameters:
t_stoppq.Quantity

Total time of the simulated processes. The events are drawn between 0 and t_stop.

ratepq.Quantity

Overall mean rate of the time series to be generated (coincidence rate coincidence_rate is subtracted to determine the background rate). Can be:

  • a float, representing the overall mean rate of each process. If so, it must be higher than coincidence_rate.

  • an iterable of floats (one float per process), each float representing the overall mean rate of a process. If so, all the entries must be larger than coincidence_rate.

coincidence_ratepq.Quantity

Coincidence rate (rate of coincidences for the n-dimensional SIP). The SIP spike trains will have coincident events with rate coincidence_rate plus independent ‘background’ events with rate rate-rate_coincidence.

n_spiketrainsint, optional

If rate is a single pq.Quantity value, n_spiketrains specifies the number of SpikeTrains to be generated. If rate is an array, n_spiketrains is ignored and the number of SpikeTrains is equal to len(rate). Default: 2

jitterpq.Quantity, optional

Jitter for the coincident events. If jitter == 0, the events of all n correlated processes are exactly coincident. Otherwise, they are jittered around a common time randomly, up to +/- jitter. Default: 0 * pq.ms

coincidences{‘deterministic’, ‘stochastic’}, optional

Whether the total number of injected coincidences must be determin- istic (i.e. rate_coincidence is the actual rate with which coincidences are generated) or stochastic (i.e. rate_coincidence is the mean rate of coincidences):

  • ‘deterministic’: deterministic rate

  • ‘stochastic’: stochastic rate

Default: ‘deterministic’

t_startpq.Quantity, optional

Starting time of the series. If specified, it must be lower than t_stop. Default: 0 * pq.ms

min_delaypq.Quantity, optional

Minimum delay between consecutive coincidence times. Default: 0 * pq.ms

return_coincidencesbool, optional

Whether to return the coincidence times for the SIP process Default: False

Returns:
outputlist

Realization of a SIP consisting of n_spiketrains Poisson processes characterized by synchronous events (with the given jitter). If return_coinc is True, the coincidence times are returned as a second output argument. They also have an associated time unit (same as t_stop).

Examples

>>> import quantities as pq
>>> import elephant.spike_train_generation as stg
>>> sip, coinc = stg.single_interaction_process(
... rate=20*pq.Hz, coincidence_rate=4.*pq.Hz,
... t_stop=1*pq.s, n_spiketrains=10, return_coincidences = True)