/* Modified Moving Average MMA(k) */
new; cls;
let data[50,1]=
      60
      61
      60
      62
      60
      63
      64
      65
      69
      70
      69
      66
      64
      64
      62
      64
      65
      55
      53
      57
      60
      58
      63
      63
      63
      68
      69
      69
      68
      69
      70
      67
      66
      67
      68
      69
      69
      77
      80
      80
      78
      76
      76
      72
      71
      70
      73
      71
      78
      80
;
data=rev(data);               /* original data in reverse order */
t=seqa(1,1,rows(data));

library pgraph;
graphset;
_pltype=6;
_plegctl=1;
_plegstr="actual\000MMA(5)\000MMA(10)";
title("Modified Moving Average");
xy(t,data~mma(data,5)~mma(data,10));


proc mma(x,k);
   local m1,i,j,w,m2,m;
/* MA(k) part */
   m1=zeros(rows(x),cols(x));
   i=k;
   do while i<=rows(x);
      j=1;
      do while j<=k;
         m1[i,.]=m1[i,.]+x[i-k+j,.];
         j=j+1;
      endo;
      i=i+1;
   endo;
   m1=m1/k;
 /* the rest */
   w=(k-(2*seqa(1,1,k)-1))/2;
   w=rev(w);
   m2=zeros(rows(x),cols(x));
   i=k;
   do while i<=rows(x);
      m2[i,.]=w'x[i-k+1:i,.];
      i=i+1;
   endo;
/* add them up */
   m=m1+6/((k+1)*k)*m2;
   m[1:k-1,.]=miss(zeros(k-1,cols(x)),0);
   retp(m);
endp;