%stat351 stat351-signrank.m
%lengths of steel rods
X=[9.83 10.09 9.72 9.87 10.04 9.95 9.82 9.73 9.79 9.90]
Y=abs(X'-10);
Z=sort(Y);

for i=1:10
   r(i)=find(Z==Y(i))
end;

>>[X' (X'-10) Y  r' sign(X'-10)]
ans =

    9.8300   -0.1700    0.1700    6.0000   -1.0000
   10.0900    0.0900    0.0900    3.0000    1.0000
    9.7200   -0.2800    0.2800   10.0000   -1.0000
    9.8700   -0.1300    0.1300    5.0000   -1.0000
   10.0400    0.0400    0.0400    1.0000    1.0000
    9.9500   -0.0500    0.0500    2.0000   -1.0000
    9.8200   -0.1800    0.1800    7.0000   -1.0000
    9.7300   -0.2700    0.2700    9.0000   -1.0000
    9.7900   -0.2100    0.2100    8.0000   -1.0000
    9.9000   -0.1000    0.1000    4.0000   -1.0000

>>find(sign(X'-10)==1)
ans =
     2
     5
>>sum(r(find(sign(X'-10)==1)))
ans =
     4                                                                      
>>sum(r(find(sign(X'-10)==-1)))
ans =
    51