计算5期移动平均的程序？

n	a
1	0.1
2	0.2
3	0.3
4	0.4
5	0.5
6	0.6
7	0.7
8	0.8
9	0.9
10	1
11	1.1
12	1.2
13	1.3
14	1.4
15	1.5

根据上面的表格求下面的表格。
aa是a的5期移动平均值。
aaa是aa的5期移动平均值。

n	a	aa	aaa	aaaa
1	0.1
2	0.2
3	0.3
4	0.4
5	0.5	0.3
6	0.6	0.4
7	0.7	0.5
8	0.8	0.6
9	0.9	0.7	0.5
10	1	0.8	0.6
11	1.1	0.9	0.7
12	1.2	1	0.8
13	1.3	1.1	0.9	0.7
14	1.4	1.2	1	0.8
15	1.5	1.3	1.1	0.9

data have;
input n a ;
datalines;
1 0.1
2 0.2
3 0.3
4 0.4
5 0.5
6 0.6
7 0.7
8 0.8
9 0.9
10 1
11 1.1
12 1.2
13 1.3
14 1.4
15 1.5
;


%macro move_avg(var,num);
   %let lagnum=%eval(&num-1);
   (&var   
     %do i=1 %to &lagnum;
       +lag&i(&var)
   %end;
   )/&num
%mend move_avg;


data wanted;
    set have;
aa=%move_avg(a,5);
aaa=%move_avg(aa,5);
aaaa=%move_avg(aaa,5);
run;

1. data a;
   
2.         do n=1 to 15 by 1;
   
3.                 a+0.1;
   
4.                 output;       
   
5.         end;
   
6. run;
   
7. 
   
8. %macro Cal(n,x,y);
   
9.         array ar_&x(&n)  _temporary_ ;
   
10.         do i=1 to &n-1 by 1;
    
11.                 ar_&x(i)=ar_&x(i+1);
    
12.         end;
    
13.         ar_&x(&n)=&x;
    
14.         &y=(ar_&x(1)%do i=2 %to &n %by 1;+ar_&x(&i)%end;)/&n;
    
15. %mend;
    
16. 
    
17. data b;
    
18.         set a;
    
19.         %CAL(5,a,aa);
    
20.         %CAL(5,aa,aaa);
    
21.         %CAL(5,aaa,aaaa);
    
22.         drop i;
    
23. run;

复制代码

学习了，俺也凑凑热闹，刚刚用R做了下

1. data<-data.frame(n=1:15,v=seq(.1,1.5,by=.1)) #你的原始数据
   
2. leng_1<-length(data$n) #原始数据的长度
   
3. leng<-5 #移动平均步长
   
4. (data2<-unlist(lapply(1:(leng_1-leng+1),function(i) mean(data$v[i+1:leng-1]))))
   
5. leng_2<-length(data2)
   
6. (data3<-unlist(lapply(1:(leng_2-leng+1),function(i) mean(data2[i+1:leng-1]))))
   
7. leng_3<-length(data3)
   
8. (data4<-unlist(lapply(1:(leng_3-leng+1),function(i) mean(data3[i+1:leng-1]))))

复制代码

1. data raw;
   
2.     input n a;
   
3. datalines;
   
4. 1 0.1
   
5. 2 0.2
   
6. 3 0.3
   
7. 4 0.4
   
8. 5 0.5
   
9. 6 0.6
   
10. 7 0.7
    
11. 8 0.8
    
12. 9 0.9
    
13. 10 1
    
14. 11 1.1
    
15. 12 1.2
    
16. 13 1.3
    
17. 14 1.4
    
18. 15 1.5
    
19. ;
    
20. /*1*/
    
21. data avg1;
    
22.     retain sum_aa sum_aaa sum_aaaa;
    
23.     set raw;
    
24.     sum_aa=sum(sum_aa,a,-lag5(a));
    
25.     aa=sum_aa/5;
    
26.     sum_aaa=sum(sum_aaa,aa,-lag5(aa));
    
27.     aaa=sum_aaa/5;
    
28.     sum_aaaa=sum(sum_aaaa,aaa,-lag5(aaa));
    
29.     aaaa=sum_aaaa/5;
    
30.     select;
    
31.         when(_n_ lt 5) call missing(aa,aaa,aaaa);
    
32.         when(_n_ lt 9) call missing(aaa,aaaa);
    
33.         when(_n_ lt 13) call missing(aaaa);
    
34.         otherwise;
    
35.     end;
    
36.     drop sum_:;
    
37. run;
    
38. /*2*/
    
39. proc expand data=raw out=avg2 method=none;
    
40.     id n;
    
41.     convert a=aa / transformin=(movave 5);
    
42. run;
    
43. proc expand data=avg2 out=avg2 method=none;
    
44.     id n;
    
45.     convert aa=aaa / transformin=(movave 5);
    
46. run;
    
47. proc expand data=avg2 out=avg2 method=none;
    
48.     id n;
    
49.     convert aaa=aaaa / transformin=(movave 5);
    
50. run;
    
51. data avg2;
    
52.     length n a aa aaa aaaa 8.;
    
53.     set avg2;
    
54.     select;
    
55.         when(_n_ lt 5) call missing(aa,aaa,aaaa);
    
56.         when(_n_ lt 9) call missing(aaa,aaaa);
    
57.         when(_n_ lt 13) call missing(aaaa);
    
58.         otherwise;
    
59.     end;
    
60. run;

复制代码

peijiamei 发表于 2010-10-9 12:46

n	a
1	0.1
2	0.2
3	0.3
4	0.4
5	0.5
6	0.6
7	0.7
8	0.8
9	0.9
10	1
11	1.1
12	1.2
13	1.3
14	1.4
15	1.5

根据上面的表格求下面的表格。
aa是a的5期移动平均值。
aaa是aa的5期移动平均值。

n	a	aa	aaa	aaaa
1	0.1
2	0.2
3	0.3
4	0.4
5	0.5	0.3
6	0.6	0.4
7	0.7	0.5
8	0.8	0.6
9	0.9	0.7	0.5
10	1	0.8	0.6
11	1.1	0.9	0.7
12	1.2	1	0.8
13	1.3	1.1	0.9	0.7
14	1.4	1.2	1	0.8
15	1.5	1.3	1.1	0.9

I give a simple solution here. You may calculate the aaa and aaaa simularly.

data tmp;
  do a=0.1 to 1.5 by 0.1;
    output;
end;
run;

data tmp2;
  set tmp;
  la0=a; la1=lag1(a); la2=lag2(a); la3=lag3(a); la4=lag4(a);
  if _n_>=5 then aa=sum(a,la1,la2,la3,la4)/n(a,la1,la2,la3,la4);
  keep a aa;
  run;

  proc print; run;

我下下来了~！免费提供大家下载！点击此处下载附件(附件大小：9KB)

5# hopewell
我想具体的算法或许这样，也就是weighted平均数的变化：

1. 
   
2. proc expand data =raw out =avg2 method =none;
   
3.     id n;
   
4.     convert a =a / transformout=(movave 5 trimleft 4);
   
5.   convert a =aa / transformout=(movave (0.04 0.08 0.12 0.16 0.20 0.16 0.12 0.08 0.04) trimleft 8);
   
6.   convert a =aaa / transformout=(movave (0.008 0.024 0.048 0.08 0.12 0.144 0.152 0.144 0.12 0.08 0.048 0.024 0.008) trimleft 12);
   
7. run;

复制代码

请问计算加权移动平均，可否在BASE里写代码实现？而不是调用PROC EXPAND过程。

mark一下