"Wyatt " <email@example.com> wrote in message news:firstname.lastname@example.org... > Hello All, > > I have encountered some weird behavior with matlab when attempting to rank > a list of numbers. > > lets assume that I generate a vector of random numbers.
> If I then sort these numbers
> and wish to obtain their ranking
> I should get back a permutation of the vector temp in terms of the ordered > indices of sorted values from A, instead this code simply returns the I > vector > > the correct execution of the command rank = temp(I) should return > > rank = [7 10 3 1 8 9 4 6 5 2] > > Am I missing something here, or mistaken about how matlab should behave in > this case? Or is there simply a bug?
You have it backwards. You want "rank(I) = temp" instead.
S = 100*rand(1, 10) [value, order] = sort(S) restoreS(order) = value isequal(restoreS, S)
"[B,IX] = sort(A,...) also returns an array of indices IX, where size(IX) == size(A). If A is a vector, B = A(IX)." [or exchanging the two sides of the equation, A(IX) = B.]
So the 1st element of B, the sorted array, corresponds to element IX(1) in A -- and that's the assignment you need to make to "undo" the sorting. To generalize, B(k) corresponds to A(IX(k)), or if you vectorize newA(IX) = B undoes the sort and makes A and newA agree.