Open r-harsha opened 11 years ago
Hi Harsha, thanks for raising the issue. Can you please be more specific? What's the input (spiketrain) and the expected output (mean rate value)? The mean_rate function extracts the time of the first and last spikes and divides this time interval by the number of spikes in the spiketrains. You might want to calculate instead the mean interspike interval (ISI) and take 1./ISI as your mean_rate. The two definitions can be much different.
Fabio
Hi Fabio, Thanks for your reply. My input is a discontinuous spike train of pulses and pauses from regular spike generator for 300 ms. I tried to obtain the mean rates over the whole population using two functions, mean_rate and firing_rate. I chose the whole stimulus duration as my time_window for firing_rate function, which makes it equivalent to mean_rate function with my t_stop is again the whole duration of 300 ms. These functions gave same results (as I expected) for the whole population of 100 neurons but they returned different results when i used only one neuron in my population. So I counted the spikes manually and found that mean_rate() function didnt give me the exact spike count for the population of one neuron.
The function monitor.sl.mean_rate() didnt give expected results for a single neuron population.