这样可以么？
%let f=0.9,0.7,0.5,0.2;  /* Input the Value.(0<=F<=.99)*/


data addyr;
set data end=end;
output;
if end then do until(yr=2005);  /*Add the End Year.*/
yr+1;  y=.;
output;
end;
run;

data Q1;
set  addyr;
retain ytop 0;
do f=&f;
                if _n_=1 then ytop=y;
                else do;
                ytop=ytop*(1-f)+y;
                end;
                aaa=0; ybar=0;
                do m=1 to _n_;
                aa= (1-f)**(m-1);
                aaa=aa+aaa;   
                end;
                ybar=ytop/aaa;
                output;
   end;
run;

data Q2;
set Q1;
retain  r1 0;   
r1=lag(ybar);
r=ybar/r1;
run;

data Q3;
set  Q2;
retain rtop 0;
do f=&f;
                if _n_=2 then rtop=r;
                else do;
                rtop=(1-f)*rtop+r;
                end;
                aaa=0; rbar=0;
                do m=2 to _n_;
                aa= (1-f)**(m-2);
                aaa=aa+aaa;   
                end;
                rbar=rtop/aaa;
          output;
        end;
run;

ouput 後很奇怪呢....

Obs     yr      y       py

                                                 1    1988    55       .
                                                 2    1988    55       .
                                                 3    1988    55       .
                                                 4    1988    55       .
                                                 5    1988    55       .
                                                 6    1988    55     55.000
                                                 7    1988    55     55.000
                                                 8    1988    55     55.000
                                                 9    1988    55     55.000
                                                10    1988    55     55.000
                                                11    1988    55     56.269
                                                12    1988    55     61.050
                                                13    1988    55     71.256
                                                14    1988    55     61.105
                                                15    1988    55     54.207
                                                16    1988    55     52.957
                                                17    1989    70     52.926
                                                18    1989    70     91.602
                                                19    1989    70     81.993
                                                20    1989    70     76.664
                                                21    1989    70     67.971
                                                22    1989    70     85.355
                                                23    1989    70     71.893
                                                24    1989    70     64.635
                                                25    1989    70     51.817
                                                26    1989    70     93.004
                                                27    1989    70     80.829
                                                28    1989    70     72.166
                                                29    1989    70     53.681
                                                30    1989    70    116.059
                                                31    1989    70    101.177
                                                32    1989    70     89.738
                                                33    1990    90     62.833
                                                34    1990    90    114.748
                                                35    1990    90     97.656
                                                36    1990    90     85.835
                                                37    1990    90     57.208
                                                38    1990    90     91.882
                                                39    1990    90     77.647
                                                40    1990    90     68.010
                                                41    1990    90     43.636
                                                42    1990    90     93.079


還有很多...

为什么会没有变量f？

8# pandayan
对于数据集中的每一行，根据每一个f的值会输出一条记录。
如原数据集中有n行，如有5个不同的f值，则会输出n*5行记录。

另外，ytop变量在计算中会被不断修改而使结果错误，需要变一下retain的变量名：

retain temp 0;

  if _n_=1 then temp=y;
  else do;
  ytop=temp*(1-f)+y;


Hope this helps!

thanks so much!!!
我研究一番!!!謝謝你的幫助!!!

二楼的解决办法很实用，谢谢！