matlab，频率计算和代码审查

Question

我有一个逻辑/统计问题：

我必须找出根据给予神经元的四种不同刺激的发射率是否显着不同。 我在Matlab中通过psth / binning方法计算了频率，但我不确定这是否正确。接下来，我通过jmp进行了方差分析和tukey测试。乍一看看起来不错，但正如我之前提到的，我认为计算不正确。

也许它不是解决我问题的合适论坛，但也许有些g可以找到我的错误或有更好的解决方案。谢谢：D

bins是通过总持续时间（800ms）除以binwidth（10ms）计算出的bin数。 最后，此函数应给我一个直方图，其中随时间（ms）和频率（此处为1x80矢量，每个仓的平均频率）绘制频率。 完成了四个不同的刺激后，我得到了四个向量，放入jmd并完成了tukey。

function [freq] = BinFireRate(data, dur, times_snippet, binwidth)
%function that plots the firing rate of a given dataset via binning method in [hz]
%in: dataset (n x m-matrix), dur as duarion observed from trial
%time_snippet (1,n-vector) for convert data into time values [ms]
%binwidth
%out: histogram of firing rate (freq) over time and frequency [hz]
%[1x80-vector] itself

bins = dur / binwidth;

spiketimes_stim = data .* times_snippet;
spiketimes_stim = spiketimes_stim(spiketimes_stim ~= 0);

[spikes_per_bin, bincenters] = hist(spiketimes_stim, bins);
freq = ((spikes_per_bin / binwidth) / length(data(:, 1))) * 1000;
bar(bincenters, freq);