SAS分组求移动平均值

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

数据集（data）概貌：
代码（id）  日期（Tradate）  成交量（Volume）  换手率（Rate）
000001        2011-01-01         1000                      1
000001        2011-01-02          200                      0.5
                                     .
                                     .
600000        2011-01-01           500                     0.1
600000        2011-01-02          3000                    0.3

data test;
set data;
if id='000001' then output;
run;
proc sort data=test;
run;
%macro test3;
data test1;
set test;
by id tradate;
array Tavg(30);
array Hrate(30);
Tavg1=Volume;
Hrate1=Rate;
%do i=2 %to 30;
           Tavg(&i.)=lag(Tavg(&i.-1));
          Hrate(&i.)=lag(Hrate(&i.-1));
%end;
/*求10、20、30天的成交量及换手率移动平均值*/
%do j=10 %to 30 %by 10;
           Avg_T&j.=round(mean(of Tavg1-Tavg&j.),1);
           Avg_H&j.=round(mean(of Hrate1-Hrate&j.),0.01);
%end;
run;
%mend test3;
%test3;
-------------------------------------------------------
问题：对求单ID=‘000001’对应的移动平均值会求，但不会对ID等于多值时进行批量处理（即，根据ID进行分组球Volume及Rate的移动平均值），请各位大虾指导，谢谢！

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：移动平均值 移动平均 平均值 volume output 换手率 平均值 成交量

上传一份样例数据一试

420948492 发表于 2014-8-15 17:32
上传一份样例数据一试

附件为样本数据，拜托了

420948492 发表于 2014-8-15 17:32
上传一份样例数据一试

如果把其中lag好的数据，在置成“.”是可以，但是效率低，同时增加代码了，非常不可取

%macro m();
data test4;
set test1;
by id tradate;
if first.id  then do;
n=1;
sum_var=volumn;
end;
else do;
n+1;
sum_var+volumn;
end;
%do inn=10 %to 30 %by 10;
lag&inn._sum=lag&inn.(sum_var);
if n<= &inn. then do ;
mean_var&inn.=round(sum_var/n,1);end;
else do;
mean_var&inn.=round((sum_var-lag&inn._sum)/&inn.,1);
end;
%end;
drop sum_var: lag:;
run;
%mend;
%m();

王者归来14 发表于 2014-8-15 17:44
附件为样本数据，拜托了

1. data data;
   
2. infile "E:\data.txt"  dlm='09'X firstobs=2;
   
3. input id :$10. Tradate YYMMDD10. Volume rate;
   
4. format tradate YYMMDD10.;
   
5. run;
   
6. 
   
7. options mprint mlogic;
   
8. %macro test3(id);
   
9. data test;
   
10. set data;
    
11. if id=&id. then output;
    
12. run;
    
13. proc sort data=test;
    
14. by id tradate;
    
15. run;
    
16. data test1;
    
17. set test;
    
18. by id tradate;
    
19. array Tavg(30);
    
20. array Hrate(30);
    
21. Tavg1=Volume;
    
22. Hrate1=Rate;
    
23. %do i=2 %to 30;
    
24. Tavg(&i.)=lag(Tavg(&i.-1));
    
25. Hrate(&i.)=lag(Hrate(&i.-1));
    
26. %end;
    
27. /*求10、20、30天的成交量及换手率移动平均值*/
    
28. %do j=10 %to 30 %by 10;
    
29. Avg_T&j.=round(mean(of Tavg1-Tavg&j.),1);
    
30. Avg_H&j.=round(mean(of Hrate1-Hrate&j.),0.01);
    
31. %end;
    
32. run;
    
33. proc append base=enddata data=test1 force;
    
34. run;
    
35. %mend test3;
    
36. 
    
37. proc sql noprint;
    
38. select count(distinct id) into:nid  from data;
    
39. create table did as
    
40. select distinct id  from data;
    
41. quit;
    
42. 
    
43. data _null_;
    
44. set did;
    
45. call symput("id_"||left(_n_),id);
    
46. run;
    
47. %macro dataend;
    
48. %do m=1 %to %eval(&nid.);
    
49. %test3(&&id_&m.);
    
50. %end;
    
51. run;
    
52. %mend;
    
53. 
    
54. %dataend;

复制代码

gongxun1234 发表于 2014-8-16 01:34
%macro m();
data test4;
set test1;

你好，非常感谢你的回复，看了你的代码，获益匪浅，我平时用惯了数组，产生很多不必要的变量，呵呵风格180度转变了要。

420948492 发表于 2014-8-16 12:11

感谢你的回复，从中又学到了不少

建议，最好用proc expand。

学习一下~