Matlab, расчет частоты и проверка кода

у меня есть логическая/статистическая проблема:

я должен выяснить, существенно ли различается скорость стрельбы в соответствии с четырьмя различными стимулами, подаваемыми на нейрон. Я рассчитал частоты с помощью метода psth/binning в Matlab, и я не уверен, что это был правильный путь. После этого я сделал anova и tukey-test через jmp. На первый взгляд все выглядит хорошо, но, как упоминалось ранее, я не думаю, что расчет был правильным.

Возможно, это не тот форум для моей проблемы, но, возможно, кто-нибудь найдет мою ошибку или предложит лучшее решение. Спасибо

bins — это количество бинов, рассчитанное по общей продолжительности (800 мс), деленное на ширину бина (10 мс). В конце эта функция должна дать мне гистограмму, построенную с частотой во времени (мс) и частотами (здесь вектор 1x80 со средней частотой на бин). Сделано для четырех разных стимулов, я получил 4 вектора, вставил jmd и сделал туки.

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);