import pynapple as nap
import numpy as np
dt = 0.001
times = np.arange(0, 4, dt)
rate = 40 * (1 + np.sin(2*np.pi*times)) / 2
spikes = nap.Ts(times[np.random.rand(len(times)) < rate * dt])
counts = spikes.count(0.1)
countsTime (s)
------------------- --
0.082 2
0.182 4
0.28200000000000003 0
0.382 0
0.482 3
0.5820000000000001 1
0.682 0
...
3.3819999999999997 2
3.4819999999999998 3
3.582 1
3.682 1
3.782 0
3.8819999999999997 0
3.9819999999999998 4
dtype: int64, shape: (40,)
