do 循环语句求移动平均代码疑惑

libname ch7 'f:\data_model\book_data\chapt7';
proc sort data=ch7.smooth;by cid month;run;
data smooth;
set ch7.smooth;
by cid;
array lags(12);
lags(1)=lag(balance);
do i=2 to 12;
    lags(i)=lag(lags(i-1));
    end;
do j=12 to 2 by -1;
if j gt cns then lags(j)=.;
  put j=;
end;
if first.cid then do
       cns=0;
       do i=1 to 12;lags(i)=.;
       end;
    end;
cns+1;
mean6=mean(of lags1-lags6);
mean12=mean(of lags1 - lags12);
/*drop i j;*/
run;

这段代码求移动平均，我没完全看懂，哪位大侠能解释下在PDV中怎样运行的？ if语句为什么不能放到前面，第二个do语句是怎么运行的?谢谢各位？

呵呵，提供另一个实现方法，不知道会不会好理解点。。。

1. data a;
   
2. input id close@@;
   
3. cards;
   
4. 1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15
   
5. 2 2 2 4 2 6 2 8 2 10 2 12 2 14 2 16 2 18 2 20 2 22 2 24 2 26 2 28 2 30
   
6. ;
   
7. run;
   
8. 
   
9. proc sort data=a;
   
10.        by id;
    
11. run;
    
12. 
    
13. %macro move_average(n,m);
    
14. data a;
    
15.         do until(last.id);
    
16.         set a;
    
17.         by id;
    
18.         retain sum;
    
19.                 sum=0;
    
20.                 %do i=0 %to &n;
    
21.                         if lag&i(id)^= id then sum=.;
    
22.                         sum=sum+lag&i(close);
    
23.                 %end;
    
24.                 move_average=sum/&m;
    
25.                 output;
    
26.         end;
    
27. run;
    
28. %mend;
    
29. %move_average(4,5);

复制代码

这个程序有点问题。每个cid的最后一个观测无法进入移动平均数的计算。也就是说，MEANS6给出的是该观测之前六个数的均值（如果之前具备6个观测的话）。京剧

求移动平均最好的办法应该是用proc expand的mvave函数吧，可以试一下

可~乐 发表于 2012-2-21 13:36
呵呵，提供另一个实现方法，不知道会不会好理解点。。。

你好！首先感谢你的这段代码！原代码用的数据见附件： smooth.txt (2.58 KB)

2012-2-24 23:19:13 上传

数据文件

你可以试试！