У меня логическая / статистическая проблема:
Я должен выяснить, существенно ли различаются скорострельность в соответствии с четырьмя различными стимулами, передаваемыми нейрону.Я рассчитал частоты с помощью метода psth / binning в Matlab, и я не уверен, что это был правильный путь.После я сделал анову и тест Тьюки через JMP.На первый взгляд это выглядит хорошо, но, как я уже говорил, я не думаю, что расчет был правильным.
Возможно, это не тот форум, на котором можно решить мою проблему, но, возможно, некоторые люди найдут мою ошибку или найдут лучшее решение.Спасибо: D
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);