问：如何在data步中快速计算分组内循环之对数和

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我有一个方法，需要编写程序，需要改内容
data a;
input x y;
cards;
1 1
1 2
1 3
1 4
1 5
2 3
2 5
2 6
2 7
;
需要按照x分组，以x=1这组来说，我需要log(2/1)+log(3/1）+log(4/1)+log(5/1)+log(3/2)+log(4/2)+log(5/2)+log(4/3)+log(5/3)+log(5/4)这样一个组合的之和，请问该如何最有效率

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：data步 Data cards Input 编写程序 编写程序 如何

本质是做有条件的笛卡儿积

HASH DATA-STEP SQL都可以做的

是不是这个结果

1. data a;
   
2. input x y;
   
3. cards;
   
4. 1 1
   
5. 1 2
   
6. 1 3
   
7. 1 4
   
8. 1 5
   
9. 2 3
   
10. 2 5
    
11. 2 6
    
12. 2 7
    
13. ; run;
    
14. 
    
15. proc transpose data=a out=aa (drop=_name_) prefix=y;
    
16.    by x;
    
17. run;
    
18. 
    
19. data aaa;
    
20.    array yy{5} y1-y5;
    
21.    set aa;
    
22.    m=1;
    
23.    do I = 1 to dim(yy)-1 while( yy[I] ^=.);
    
24.         do J = I+1 to dim(yy) while (yy[J] ^=.);
    
25.            m = m*yy[J]/yy[I];
    
26.    end;  end;
    
27.    Log_value = log10(m);
    
28. 
    
29.    do I = 1 to dim(yy) while(yy[i]^=.);
    
30.       y = yy[i];
    
31.       output;
    
32.    end;
    
33.    drop i j y1-y5 m;
    
34. run;

复制代码

爱萌挖了个大坑，让大家跳呢。
log(a*b)=loga+logb吧。

1. data a;
   
2. input x y;
   
3. cards;
   
4. 1 1
   
5. 1 2
   
6. 1 3
   
7. 1 4
   
8. 1 5
   
9. 2 3
   
10. 2 5
    
11. 2 6
    
12. 2 7
    
13. ;
    
14. 
    
15. data b;
    
16.         do count=0 by 1 until(last.x);
    
17.                 set a(keep=x);
    
18.                 by x;
    
19.         end;
    
20.                 cumsum=0;
    
21.         do _n_ = -count to count by 2;
    
22.                 set a;
    
23.                 cumsum = cumsum + _n_*log(y);
    
24.         end;
    
25.         keep x cumsum;
    
26. run;

复制代码

邓贵大 发表于 2013-12-27 03:32

处理的很巧妙，非常感谢，最少的代码干最多的事情

用笛卡尔积就行，只不过数据量不能很大了

proc sql noprint;
        create table result as
                select a.x,a.y,b.x as x1,b.y as y1
                from a as a,a as b
                where a.x = b.x and a.y < b.y;
quit;

proc sql noprint;
        create table result1 as
                select distinct x,SUM(LOG(y1)/LOG(y)) as logsum
                from result
                group by x;
quit;

playmore 发表于 2013-12-27 13:06
用笛卡尔积就行，只不过数据量不能很大了

proc sql noprint;

楼主貌似将SUM(LOG(y1/y))敲成SUM(LOG(y1)/LOG(y))了