循环部分都什么意思啊？
proc iml ;
A={1363.6 16.86 16.49 10.52 13.39 75364 36.5 371.1 12.8 37031.1 93987 4450.7,
1385.1 19.2 17.4 9.2 12 69807 38.78 338.9 12.9 36349.4 86546 5720,
1423.2 16.1 15.5 7.5 10.8 79322 40.57 370.4 13.2 43396 96435.5 6142,
1456.4 18.4 13.8 8.4 13 71583 40.87 444.9 12.9 54443 108668.6 5958,
1492.7 21.4 13.5 7.7 12.8 82986 41.91 483.5 13.5 55019 102076 5171,
1538 23.2 13.2 6.8 13.4 71600 42 410.7 13.2 70158 112432.7 10435,
1581 24.5 12.8 6.2 13.7 74970 42.5 318 13.4 93198 129138 13844,
1633 23.8 12.4 5.8 13.9 77778 43 483.9 14 98865 129820 16513,
1695 34.2 12 5.2 14.7 80793 43.5 626.3 13.4 104255 132095 19559,
1755 21.8 12 5.2 14.7 86881 44.4 480.6 13.3 109604 136511 11190
};
r=nrow(A);
c=ncol(A);
B=A;
do i=1 to c;
do p=1 to r;
B[p,i]=(A[p,i]-A[><,i])/(A[<>,i]-A[><,i]);
end;
end;
F=B;
do i=1 to ncol(B);
do p=1 to nrow(B);
F[p,i]=B[p,i]/B[+,i];
end;
end;
do i=1 to ncol(B);
do p=1 to nrow(B);
if F[p,i]=0
then F[p,i]=1;
end;
end;
HF=F#log(F);
HF=HF[+,];
H=HF/-log(r);
W=(1-H)/(c-sum(H));
print H W;
quit;
