Использование теста Дипа Хартигана для расчета бимодальности

Question

Я пытаюсь использовать Hartigans Dip Test в Matlab для вычисления бимодальности.У меня есть данные о полуширинах потенциалов действия, и мне кажется, что они имеют 2 пика (см. Прилагаемые примеры).Результаты почти одинаковы для обоих случаев.Спасибо!

бимодальный пример

вот код, который я использую: nicprice.net/diptest

halfwidth = zeros(1,size(spikes,1));

  for j = 1:size(spikes,1)                 
    new = spikes(j,:);                     
    half = calc_halfwidth(new);

if isempty(half)==1
    halfwidth(j) = 0;
else
    halfwidth(j) = half;
end     
 end

[f,xi] = ksdensity(halfwidth);

[dip,p] = HartigansDipTest(f)

унимодальный