某事件发生时，前三个，后三个的累乘祈求@jingju11等大牛现身

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

附件中有这样一组数据，小弟想算出每当facpr不为0或者空值是他前三个return和后三个ret的累乘值，请问各位大牛不吝指教：）jingju大牛现身啊：）   

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：ING RETURN 不吝指教 turn RET return

@jingju11, 大人现身呀：）

顶起：）

肯定有更有效率的写法,也没有仔细验证.权且抛砖引玉.
京剧

1. data _null_;
   
2.    call symputx('nobs', nobs); set have nobs =nobs; stop;
   
3.    run;
   
4. data have2;
   
5.    array t[&nobs, 3] _temporary_; array r[&nobs,2] _temporary_;
   
6.    if _n_ =1 then do;
   
7.       do p =1 to &nobs;
   
8.          set have point =p;
   
9.          t[p,1] =p; t[p,2] =ret; t[p,3] =facpr;
   
10.          end;
    
11.       do p =1 to &nobs;
    
12.          if t[p,3] then do;
    
13.             if p >3       then r[p,1] =t[p-1,2]*t[p-2,2]*t[p-3,2];*3 days before;
    
14.             if p <&nobs-3 then r[p,2] =t[p+1,2]*t[p+2,2]*t[p+3,2];*3 days after-;
    
15.             end;
    
16.          end;      
    
17.       end;
    
18.    set have;
    
19.    r_b3 =r[_n_,1]; r_a3=r[_n_, 2];
    
20. run;

复制代码

是照着某大牛的葫芦画瓢，自己也不知道其中的逻辑，最后错位了也不知道怎么改，还请多多指点。
还数@jingju11大人的正统。

1. data tt;
   
2.     set t;
   
3.         array r{7} _temporary_;
   
4.         do i=1 to 6;
   
5.             r{i}=r{i+1};
   
6.         end;
   
7.         r{7}=ret;
   
8.         if lag3(FACPR) gt 0 then do;
   
9.             lead_3=r{1}*r{2}*r{3};
   
10.             lag_3=r{5}*r{6}*r{7};
    
11.         end;
    
12.         drop i;
    
13. run;

复制代码

boe 发表于 2013-7-29 01:57
是照着某大牛的葫芦画瓢，自己也不知道其中的逻辑，最后错位了也不知道怎么改，还请多多指点。
还数@jingj ...

很喜欢这样巧妙的程序.

一点意见:如果错位不算什么,有一点或许值得注意:比如,如果最后一条记录的的FACPR=1的话,应该可以算出之前的3个乘积,但是因为程序滞后三位的原因,也将无法显示.
个人意见.也没有充分考虑.
最最重要的是,我不是大牛...

京剧

jingju11 发表于 2013-7-29 02:28
很喜欢这样巧妙的程序.

一点意见:如果错位不算什么,有一点或许值得注意:比如,如果最后一条记录的的F ...

您说的很有道理啊，我确实没有考虑最后一条记录大于零的问题。
改成下面这样可以解决，但lag_3还是错位，可能真得用merge把r5-r7先上移3行。
不知道有没有更好的方法，还请您多指点。
您的方法已经很完美了。
真是越牛越谦虚，佩服！

1. data tt;
   
2.     set t;
   
3. 
   
4.         array r{7} _temporary_;
   
5.         do i=1 to 4; r{i}=r{i+1}; end; r{4}=ret;
   
6.         do i=5 to 6; r{i}=r{i+1}; end; r{7}=ret;
   
7. 
   
8.         if FACPR gt 0 then lead_3=r{1}*r{2}*r{3};
   
9.         if lag3(FACPR) gt 0 then lag_3=r{5}*r{6}*r{7};
   
10. 
    
11.         drop i;
    
12. run;

复制代码

谢谢两位大牛的精彩答案，小弟受益匪浅：）

boe 发表于 2013-7-29 05:12
您说的很有道理啊，我确实没有考虑最后一条记录大于零的问题。
改成下面这样可以解决，但lag_3还是错位 ...

感谢指点.修改了一下程序:京剧

1. data _null_;
   
2.    call symputx('nobs', nobs); set have nobs =nobs; stop;
   
3.    run;
   
4. data have1;
   
5.    array t[&nobs, 2] _temporary_;
   
6.    if _n_ =1 then do p =1 to &nobs;
   
7.       set have point =p;
   
8.       t[p,1] =p; t[p,2] =ret;
   
9.       end;
   
10.    set have;
    
11.    if facpr then do;
    
12.       if _n_ >3       then r_b3 =t[_n_-1,2]*t[_n_-2,2]*t[_n_-3,2];*3 days before;
    
13.       if _n_ <&nobs-3 then r_a3 =t[_n_+1,2]*t[_n_+2,2]*t[_n_+3,2];*3 days after-;      
    
14.       end;
    
15.    run;

复制代码

强大！
好好学习一下。