sum的问题求助

shiyiming

data temp;
input id y;
cards;
1 104.932
2 107.732
3 103.121
4 103.340
5 104.653
6 101.401
7 95.430
8 98.784
9 102.981
10 98.921
11 94.178
12 90.831
13 98.756
14 93.079
15 99.559
16 99.036
17 98.258
18 90.956
19 89.173
20 89.344
run;
我想做的就是
对于第i个id，算出sum=sum_{yj>yi}（yj ）
我能想到的就是用do loop，但是结果不对

这个要怎么做呢？

shiyiming

下面是不是你要的结果：

data temp;
input id y;
cards;
1 104.932
2 107.732
3 103.121
4 103.340
5 104.653
6 101.401
7 95.430
8 98.784
9 102.981
10 98.921
11 94.178
12 90.831
13 98.756
14 93.079
15 99.559
16 99.036
17 98.258
18 90.956
19 89.173
20 89.344
;
run;

%macro sm;

%let nfile=%sysfunc(open(temp));
%let nrw=%sysfunc(attrn(&nfile,nlobs));
%let rc=%sysfunc(close(&nfile));

proc iml;
use work.temp;
read all into mx;

nrw=nrow(mx);

%do i = 1 %to &nrw;
  s=0;
  m1=mx[&i,2];
  %do j=1 %to &nrw;
    m2=mx[&j,2];
    if (m2>m1) then do;
      s=s+m2;
        end;
    print s;
  %end;
%end;
quit;
%mend;
%sm;