How to properly use accumarray in MATLAB with sorted indices and weights?

  • Context: MATLAB 
  • Thread starter Thread starter atrus_ovis
  • Start date Start date
  • Tags Tags
    Matlab
Click For Summary
SUMMARY

The discussion focuses on the proper use of the accumarray function in MATLAB, particularly when handling sorted indices and weights. The user encountered discrepancies between the output of accumarray and the unique values of the input indices. Key insights include understanding how accumarray processes input indices and weights, creating bins based on unique indices, and accumulating values accordingly. The output size is determined by the maximum index value, leading to potential zero entries for indices not represented in the input.

PREREQUISITES
  • Familiarity with MATLAB programming language
  • Understanding of the accumarray function and its parameters
  • Knowledge of sorting algorithms and their implementation in MATLAB
  • Basic concepts of unique values and binning in data processing
NEXT STEPS
  • Explore advanced features of accumarray in MATLAB documentation
  • Learn about MATLAB's data types and how they affect function performance
  • Investigate the use of sort and unique functions in MATLAB for data manipulation
  • Study examples of histogram creation using accumarray for better visualization of data distribution
USEFUL FOR

Data analysts, MATLAB users, and researchers working with numerical data who need to efficiently aggregate and analyze weighted indices.

atrus_ovis
Messages
99
Reaction score
0
accumarray documentation
Is anyone proficient in the use of accumarray?
I supply two vectors of same length, indices L,weights W, L positive integers as required.
I sort the indices to ascending.Then, the result of accumarray is not equal to unique(L).How could this be?
Code for illustration:
Code:
i = [46    47    47    46    48    49    48    48    48]';
w = 2*rand(size(i));
Y=[i';w'];
[x y] = sort(Y(1,:));
Y = Y(:,y);
Y,pause
uw = accumarray((Y(1,:))',Y(2,:)');

ui = unique(i);
numel(uw),numel(ui)

--- OUTPUT :

Y =

   46.0000   46.0000   47.0000   47.0000   48.0000   48.0000   48.0000   48.0000   49.0000
    1.5844    0.0714    1.9190    1.3115    1.6983    1.3575    1.5155    1.4863    1.8680ans =

    49ans =

     4
I get 49 weight vectors for the 9 indices, the unique of which are 4!
 
Physics news on Phys.org
It took me a few passes to read through the documentation and understand it, but I think the key is where it explains what the function actually does (and then apply that to the first example).

MATLAB accumarray documentation said:
The function processes the input as follows:

  1. Find out how many unique indices there are in subs. Each unique index defines a bin in the output array. The maximum index value in subs determines the size of the output array.
  2. Find out how many times each index is repeated.
    This determines how many elements of vals are going to be accumulated at each bin in the output array.
  3. Create an output array. The output array is of size max(subs) or of size sz.
  4. Accumulate the entries in vals into bins using the values of the indices in subs and apply fun to the entries in each bin.
  5. Fill the values in the output for positions not referred to by subs. Default fill value is zero; use fillval to set a different value.

In the case of the example, the output size is the maximum of the subs array--4. The subs array acts as an array of subscripts applied to the values contained in val. 101 is given a subscript of 1, 102 and 104 given subscripts of 2, 103 and 105 given subscripts of 4, and none given subscripts of 3. Bin labels, if you want to think of it in histogram terms.

When you apply accumarray to them, the first entry will be the sum of those elements with subscript 1, the second with subscript 2, and so on. Which is why the output array has a sum of zero for the third element (none of the elements of val were subscripted with three).

In your case:
  • uw has 49 elements because that's the largest of the values contained in i.
  • uw probably only has six non-zero elements--uw(46) through uw(49)
  • uw(1:45) are zero

I hope this helps in your understanding of what's going on in this portion of your work, because I don't have any experience with the larger problem that you're working on.
 
Yes i got it eventually.
In my example,to add you just use the weights
uw(ui)
 

Similar threads

  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 1 ·
Replies
1
Views
9K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 2 ·
Replies
2
Views
15K
  • · Replies 12 ·
Replies
12
Views
9K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 9 ·
Replies
9
Views
5K