Ranks calculation

[robmio]

This version works fine.

Code:


#cas
ranking(m):=
BEGIN
LOCAL hi, m2, m3, lo, j, ii, r, k;
hi:= length(m);
m2:=zip(m, range(hi).+1);
m2:=sort(m2);
lo:=1;
m3:=m;
FOR ii FROM 1 TO hi DO
j:=ii+1;
IF ii≠hi AND m2[ii][1]==m2[j][1] THEN
CONTINUE;
END;
r:=(lo+ii)/2;
FOR k FROM lo TO ii DO
m3[m2[k][2]]:=r;
END;
lo:=j;
END;
RETURN m3;
END;
#end