histc and bin-width determination
显示 更早的评论
Hi
Given I have the following data:
A = [100 150 190 200 250 300 350 370 390 400 400]
And I want 3 bins: From 100 to 200, 200 to 300, and 300 to 400.
How do I do that with histc?
What I tried so far was:
edges = linspace(min(stim_durations),max(stim_durations),4);
counts = histc(stim_durations, edges);
But that results in 4 bins with "400" having its own bin..how can I solve that?
采纳的回答
The simplest solution on a recent enough version of matlab (2014b or newer) is to use histcounts which behaves exactly as you want:
>>A = [100 150 190 200 250 300 350 370 390 400 400];
>>hiscounts(A, [100 200 300 400])
ans =
3 2 6
As per histc documentation, the last edge is its own bin and you can't do anything about that. If you want 400 to be included in the third bin and still use histc, you have to shift the edge of the 4th value slightly above 400 (and discard that last bin):
A = [100 150 190 200 250 300 350 370 390 400 400];
edges = linspace(min(A),max(A),4);
edges(end) = edges(end) + eps(edges(end)); %or any value greater than eps(400);
counts = histc(A, edges);
counts = counts(1:end-1)
3 个评论
Walter Roberson
2015-5-7
Another approach other than adding eps() is to add another bin at the edge with value of infinity, and discard the 0 count that occurs for it.
edges(end+1) = inf;
MiauMiau
2015-5-7
Image Analyst
2015-5-7
histc IS available in older versions. It's now been deprecated in favor of histcounts, meaning that in new versions they recommend you use histcounts and not histc. I didn't mention it in my answer because it seems like hardly anyone except me has the latest version.
更多回答（1 个）
Walter Roberson
2015-5-7
edges = linspace(min(stim_durations),max(stim_durations),4);
edges(end) = [];
That is, construct 4 bins and discard the last of them. The final bin will be everything from the previous value upwards.
Please pay attention to whether you want your bins to be 100 <= x <= 200, 200 <= x <= 300, or if you want them to be 100 <= x < 200, 200 <= x < 300 or 100 < x <= 200, 200 < x <= 300 . Those are slightly different boundary conditions that can be pretty meaningful.
类别
在 帮助中心 和 File Exchange 中查找有关 Data Distribution Plots 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
