import quantities as pq
import matplotlib.pyplot as plt
from elephant.spike_train_generation import homogeneous_poisson_process
from viziphant.statistics import plot_isi_histogram
np.random.seed(12)
rates = [5, 10, 15] * pq.Hz
spiketrains = [homogeneous_poisson_process(rate=r,
               t_stop=100 * pq.s) for r in rates]
plot_isi_histogram(spiketrains, cutoff=250*pq.ms,
                   legend=rates)
plt.show()