Ticket #22: fano_factor.py
| File fano_factor.py, 0.5 KB (added by pierre, 4 years ago) |
|---|
| Line | |
|---|---|
| 1 | from NeuroTools import signals |
| 2 | import numpy |
| 3 | |
| 4 | spikes = [] |
| 5 | nb_cells = 100 |
| 6 | frequencies = nb_cells*[5] |
| 7 | for idx in xrange(nb_cells): |
| 8 | param = 1./frequencies[idx] |
| 9 | isi = numpy.random.exponential(param, 100) |
| 10 | pspikes = numpy.cumsum(isi)*1000. # To convert the spikes_time in ms |
| 11 | for spike in pspikes: |
| 12 | spikes.append((idx, spike)) |
| 13 | spk = signals.SpikeList(spikes,range(10)) |
| 14 | |
| 15 | print spk.fano_factor(1) |
| 16 | print spk.fano_factor(5) |
| 17 | print spk.fano_factor(10) |
| 18 | print spk.fano_factor(100) |
