How to make a weighted histogram with specific bins?

Question

Jaclyn 2013-7-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/81805-how-to-make-a-weighted-histogram-with-specific-bins

I'm new and I use MATLAB 2011.

I have a column of data and another column of weights:

d = [33 11 18 ... ] %data
w = [1 0.5 1 ...] %weights

There are only 11 possible values of my data so there would be 11 bins, I need to know the frequency count of each value.

I got a plot to work nicely WITHOUT including the weights by changing the data to ordinal:

bins = [0 11 12 13 14 15 16 17 18 32 33];
counts = droplevels(ordinal(d,[],bins));
hist(counts);
set(gca,'XTick',1:11);

I changed it to ordinal because otherwise, the bars had large spaces between them (because the x-axis ranged from 0-33) or merged together (clumped 11-18 together as one). I tried so many things, I couldn't list them all here.

The point is I need a plot to INCLUDE the weights in the frequency counts. I assume it can't be ordinal so the above code is irrelevant. I've done a lot of googling, everything I tried hasn't worked.

Any help is appreciated it, sorry if I'm confusing.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Jaclyn 2013-7-11

Each value in d has a 1 or a 0.5 associated with it as the weight. Values with a 0.5 are considered half an occurrence. For example, if 11 comes up five times in the dataset but one is weighted 0.5, the frequency count in the histogram for 11 will be 4.5.

dpb 2013-7-11

在 MATLAB Online 中打开

So

length(w) == length(unique(data))

???

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2013-7-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/81805-how-to-make-a-weighted-histogram-with-specific-bins#answer_91517

编辑：dpb 2013-7-11

在 MATLAB Online 中打开

edges=unique(data);          % vector of unique values
n=histc(data,unique(data));  % to get bin counts w/ bins specified
bar(edges,n,'histc')         % plot

Use weights how desired...[edited to add]

If the answer to the question just posed is "yes" then

n=n.*w;

I don't know what bar() will do w/ a half for a count; never tried it. If it doesn't like non-integer values you can either floor() or round() the results depending on which you think is better/more appropriate for you problem.

bar() is ok w/ the fractional count; if you still have trouble w/ the bar location because the x-axis is positioned by value of the bin rather than labeling the bin you can just use x=1:11 and then

set(gca,'xticklab', num2str(bins'))

to label them by the bin numbers. It's a key point w/ handle graphics axes object that the xtick values and the labels do not have to be the same.

doc histc % and friends for further details

5 个评论
显示 3更早的评论隐藏 3更早的评论

dpb 2013-7-12

编辑：dpb 2013-7-14

在 MATLAB Online 中打开

c is the index into data of each unique value in data...it's what accumarray is grouping the weights over to get the total weight per bin You need to arrange it in the accumarray call that c is a column vector for reasons having to do w/ how it functions; not to worry over the details at the moment, it just needs must be a column. In my sample where I tested I had a row, hence the transpose. If you have a column already owing to the orientation of data, then remove that. The weights then also must be a column vector to coincide.

In fact going back I see you say you do have column vectors; my bad -- remove the transpose [ ' ] from both and joy should ensue.

nwt=n.*accumarray(c,w);

Here's a toy example small enough you can see what's going on explicitly...

>> d=[1 1 3]'; w=[1 .5 1]';
>> [n,ib]=histc(d,unique(d));
>> [u,~,c] = unique(d);
>> [n,ib]=histc(d,u);
>> nwt=n.*accumarray(c,w)
nwt =
   3
   1
>> n
n =
   2
   1
>>