使用Hartigan的Dip测试计算双峰

Question

我正在尝试在matlab中使用Hartigans Dip Test来计算双峰态。我有动作电位半角宽度的数据，在我看来它有2个峰值（请参阅附带的示例）。两种情况下的结果几乎相同。 谢谢！

bimodal example

这是我正在使用的代码：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)

unimodal