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 上传

数据文件

你可以试试！

楼主能否把原始数据给我看看，否则我不知道数据中有什么变量~

%let p=12;
data e (drop=i j);
        set a;
        if balance ne '.';
        n+1;
        by cid;
/*        (n指代每组的序号)*/
        if first.cid then n=1;
        array a[&p.] b1-b&p.;
        a[1]=lag(balance);
        do i=2 to &p.;
                a[i]=lag(a[i-1]);
        end;
       
        j=&p.;
        do until (j=0);
        if(n<=j) then a(j)=.;
        j=j-1;
        end;

        mean6=mean(of b1-b6);
        mean12=mean(of b1-b12);
run;

您可以尝试上面这段代码

do j=12 to 2 by -1;
if j gt cns then lags(j)=.;
  put j=;
end;

我猜楼主是在说这段代码？
这段代码很精彩，比如说刚开始，j=12，在这个时候，1<12，则lags[12]第一行数据=.；同理2、3、4、5、6、7、8、9、10、11都小于12，因此lags[12]的前11个数据都为.；
但是我认为，楼主要求moving average by 12，因此这里if j gt cns应改写为if j ge cns，这样就可以实现lags[12]前12个数据都为.了，您觉得呢？
然后当j=11时，lags[11]的前10个数据都为.；
一次类推，希望能够帮到您

现在我们来看楼主的第一个问题：
if j gt cns then lags(j)=.;
  put j=;
end;
这个if是不可以放在do i=12 to 2 by -1前面的，否则j gt cn无从谈起；
if first.cid then do
       cns=0;
       do i=1 to 12;lags(i)=.;
       end;
    end;
cns+1;

这个if的位置是为了找到每一组的第一条数据，然后为每组数据进行赋值；但是如果放在后面也不是不可以，我想如果放在后面的话程序跑不出cns的数值，那加一个retain应该可解决问题。


希望可以帮到您