sas 分组连乘

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

sas如何做到分组连乘啊？希望各位大师赐教。。。谢谢哈！

data a; imput x1 x2; cards;
1 2
2 2
3 2
4 3
5 3
6 3
;
run;

我想得到
x2   x3
2     6
3   120
也就是按照x2分组，将x1连乘。。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：分组连乘 cards card Data put 如何

1. data a; input x1 x2; cards;
   
2. 1 2
   
3. 2 2
   
4. 3 2
   
5. 4 3
   
6. 5 3
   
7. 6 3
   
8. ;
   
9. run;
   
10. 
    
11. proc sort data=a out=b;by x2;run;quit;
    
12. 
    
13. data c;
    
14.         retain x2 x3 1;
    
15.         set b;
    
16.         by x2;
    
17.         if first.x2 then x3=1;
    
18.         x3=x3*x1;
    
19.         if last.x2 then output;
    
20.         keep x2 x3;
    
21. run;
    
22. 
    
23. proc print;run;quit;

复制代码

2# soporaeternus2

谢谢哈，如果按照两个变量分组呢！

data a; imput x1 x2 x3; cards;
1 2 1
2 2 1
3 2 2
4 3 1
5 3 2
6 3 2
;
run;

我想得到
x2   x3    x4
2     1     2
2     2     3
3     1     4
3     2     30
也就是按照x2和x3分组，将x1连乘。。

1. data a; input x1 x2 x3; cards;
   
2. 1 2 1
   
3. 2 2 1
   
4. 3 2 2
   
5. 4 3 1
   
6. 5 3 2
   
7. 6 3 2
   
8. ;
   
9. run;
   
10. 
    
11. 
    
12. 
    
13. proc sort data=a out=b;by x2 x3;run;quit;
    
14. 
    
15. data c;
    
16.         retain x2 x3 x4 1;
    
17.         set b;
    
18.         by x2 x3;
    
19.         if first.x3 then x4=1;
    
20.         x4=x4*x1;
    
21.         if last.x3 then output;
    
22.         keep x2 x3 x4;
    
23. run;
    
24. 
    
25. proc print;run;quit;

复制代码

......以此类推了，SQL应该也能写，一下子想不起来了......

4# soporaeternus2

这已经很好了，谢谢哈。。。对sas的编程还是不怎么熟练啊！

5# wdxmahone

SQL应该也能......

For the 2nd example:

1. 
   
2. proc sql;
   
3.     select x2, x3,
   
4.          case
   
5.       when min(abs(x1)) > 0 then (-1)**sum((x1<0))*exp(sum(log(abs(x1))))
   
6.            else 0
   
7.           end as product
   
8.       from A
   
9.       group by x2, x3;
   
10. quit;

复制代码

The basic idea is that the log product is equal to the sum of log values when the values are all positive, thus you need to consider 1. the sign 2. if 0.
One of the pros of using SQL is its much more understandable syntax here.

1. data a; input x1 x2; cards;
   
2. 1 2
   
3. 2 2
   
4. 3 2
   
5. 4 3
   
6. 5 3
   
7. 6 3
   
8. ;
   
9. run;
   
10. 
    
11. proc sort data=a out=b;
    
12. by x2;
    
13. run;
    
14. 
    
15. proc transpose data=b out=c;
    
16. var x1;
    
17. by x2;
    
18. run;
    
19. 
    
20. data d (keep=x2 x3);
    
21. set c;
    
22. x3=col1*col2*col3;
    
23. run;

复制代码

6# jingju11
恩，就是这个意思

1. 
   
2. log(abs(x1)

复制代码

改成

1. 
   
2. log(case when abs(x1)>0 then abs(x1) else 1 end)

复制代码

就没有那个log函数参数无效的NOTE了......

学习了，好强大，谢谢

6# jingju11

学习了 太牛了

6# jingju11

这种转换太强大了。。。这是一个牛人辈出的年代。。。