求助SAS有条件下的求和

%macro test(n=);
data ex1;
        merge ex %do i=2 %to %eval(&n+1); ex(keep=b rename=(b=b&i) firstobs=&i) %end;;
        sumb=sum(of b:);
run;
proc sql noprint;
        create table want as
        select a,  sumb, sum(sumb) as sumall
        from ex1
        having a=1;
quit;
%mend;

%test(n=5);

1. data ex ex1(rename=(b=c) keep=b);
   
2. do i =1 to 100;
   
3. a=round(ranuni(0),1);
   
4. b=ranuni(0);
   
5. output;
   
6. end;
   
7. drop i;
   
8. run;
   
9. 
   
10. data want(drop=i c);
    
11. set ex;
    
12. sum=0;
    
13. if a then do;
    
14.         do i=_N_ to min(_N_+4, aa);
    
15.                 set ex1 nobs=aa point=i;
    
16.                 sum+c;
    
17.         end;
    
18. end;
    
19. sum1+sum;
    
20. run;

复制代码

jackyloda 发表于 2014-12-23 11:23
这个倒不必，我主要是想知道算法的实现，像2L那样随机生成的做Sample就行。主要是觉得sas里每个变量的第几 ...

1. /* It is not easy to understant what do you want */
   
2. /* the following code produce one number of sum */
   
3. 
   
4. data want;
   
5.   array aa{200} _temporary_;
   
6.   array bb{200} _temporary_;
   
7.   if _n_ = 1 then do;
   
8.      i = 1;
   
9.    do until (last);
   
10.       set ex end=last;
    
11.         aa[i] = a;
    
12.         bb[i] = b;
    
13.      i+1;
    
14.      end;
    
15.      end;
    
16.    sum_a=0;
    
17.     do i = 1 to 200;
    
18.       a = aa[i]; b=bb[i];
    
19.     if a = 1 then do;
    
20.         sum_6b= sum(bb[i], bb[i+1], bb[i+2], bb[i+3],bb[i+4], bb[i+5]) ;
    
21.         sum_a + sum_6b;
    
22.         put i= aa[i]= sum_6b=  sum_a=;
    
23.     end;
    
24.     end;
    
25. run;

复制代码

1. /*举个例子*/
   
2. data a;
   
3. input A B @@;
   
4. cards;
   
5. 1 1 1 2 1 3 1 4 1 5 1 8 1 3
   
6. 0 2 0 3
   
7. ;
   
8. run;
   
9. data b;
   
10. set a;
    
11. id=_n_;
    
12. run;
    
13. proc sql;
    
14. create table c as select b.id,b.A,b.B,b1.B as B1 from b,b as b1 where 0<=b1.id-b.id<=5 and b.A=b1.A and b.A=1 order by b.id;
    
15. create table d as select distinct id,A,B,sum(B1) as sum_5,count(*) as count
    
16.                   from c group by A,id; /*计算A=1时相应B在这个位置的值与后5个的B的和，并统计个数，注意后面有些未满6个，特用count标注*/
    
17. create table e as select distinct A,sum(sum_5) as sum_all
    
18.                   from d;/*最终结果，注意其中包括未满6个的B对应的和，如果需要删除不满足条件的count，可以考虑加where count<=6*/
    
19. quit;

复制代码

jackyloda 发表于 2014-12-23 09:42
你好，第二段的代码我还不是太理解，你是一个什么思路呢，g是一个数组吗？

g是一个矩阵，也可以当数组来理解，loc函数找出所有A=1的位置，比如有一个位置在N，求N：N+5数据的和，
当然也要考虑N+5与总长度的大小，如果N+5超过总长度会报错，所以取两个之间较小的那个

李会超 发表于 2014-12-23 12:03
这个用二楼那哥们的pro iml肯定更简洁，只是很多人不太懂iml模块，我这个思路容易理解，就是代码过长，如 ...

可以简单介绍下Iml模块的思路吗，或者哪里有比较简单介绍这个的，我身边的SAS书都没讲到这个

ziyenano 发表于 2014-12-22 23:50
data ex;
do i =1 to 100;
a=round(ranuni(0),1);

非常感谢你的答案，对我帮助很大！还有几个问题麻烦解答下
res=j(ncol(location),1,0);
这个j是什么意思呢，另外ncol是什么函数？

jackyloda 发表于 2014-12-24 10:43
非常感谢你的答案，对我帮助很大！还有几个问题麻烦解答下
res=j(ncol(location),1,0);
这个j是什么意思 ...

j(m,n,a)生成m行 n列 元素全为a的矩阵，ncol(M)返回矩阵M的行数

ziyenano 发表于 2014-12-24 12:23
j(m,n,a)生成m行 n列 元素全为a的矩阵，ncol(M)返回矩阵M的行数

了解了，多谢！以后有SAS问题再多多请教你！

ziyenano 发表于 2014-12-24 12:23
j(m,n,a)生成m行 n列 元素全为a的矩阵，ncol(M)返回矩阵M的行数

你好，如果现在问题变为还有一个C列也想求和，代码能写成类似这样么，
g[location,2]+g[location:location+5,3][+]这种形式吗？