%stat351 stat351-rank.m
%Rank order statistics
>> X=unifrnd(0,1,10,1)
X =[
0.9501
0.2311
0.6068
0.4860
0.8913
0.7621
0.4565
0.0185
0.8214
0.4447]

>> Y=sort(X)
Y =
0.0185
0.2311
0.4447
0.4565
0.4860
0.6068
0.7621
0.8214
0.8913
0.9501

for i=1:10
   r(i)=find(Y==X(i))
end;
r =
    10     2     6     5     9     7     4     1     8     3

X=[-1:0.1:1];
Y=sqrt(1-X.^2);
plot(X,Y,'+');
>>corrcoef(X,Y)
ans =
    1.0000    0.0000
    0.0000    1.0000

%wine experts
wine=[1 2
      2 1
      3 5
      4 3
      5 4
      6 7
      7 8
      8 6];
plot(wine(:,1),wine(:,2),'o')
>>corrcoef(wine(:,1),wine(:,2))
ans =
    1.0000    0.8333
    0.8333    1.0000

>>corrcoef(X,r)
ans =
    1.0000    0.9780
    0.9780    1.0000

plot(X,r,'+')