from NeuroTools import signals
import numpy

spikes = []
nb_cells = 100
frequencies = nb_cells*[5]
for idx in xrange(nb_cells):
    param   = 1./frequencies[idx]
    isi     = numpy.random.exponential(param, 100)
    pspikes = numpy.cumsum(isi)*1000. # To convert the spikes_time in ms
    for spike in pspikes: 
        spikes.append((idx, spike))
    spk = signals.SpikeList(spikes,range(10))
        
print spk.fano_factor(1)        
print spk.fano_factor(5)
print spk.fano_factor(10)
print spk.fano_factor(100)