call symput 多次使用问题

我在用proc nlin做nelson siegel 模型，在估计rate时遇到了以下问题在用如下程序生成了Parmds表格后

1. %macro NS;
   
2. %do i=1 %to 3;
   
3. data para;
   
4. set rate(keep=rate&i);
   
5. set time;
   
6. run;
   
7. 
   
8. proc nlin data=para method=newton;
   
9. parms b0 = 0.05
   
10.       b1 = 0.05
    
11.       b2 = 0.05
    
12.       b3 = 0.05
    
13.       tao1 = 1
    
14.       tao2 = 1
    
15.    ;
    
16. 
    
17. model rate&i=b0+b1*(1-exp(-ytm/tao1))/(ytm/tao1)+b2*(1-exp(-ytm/tao1))/((ytm/tao1)-exp(-ytm/tao1))+b3*(1-exp(-ytm/tao2))/((ytm/tao2)-exp(-ytm/tao2));
    
18. ods output parameterestimates = parmd;
    
19. run;
    
20. 
    
21. 
    
22. data parmds;
    
23.   %if &i^=1
    
24.   %then %do;
    
25. set parmds;
    
26.   %end;
    
27. set parmd(rename=(estimate=estimate&i));
    
28. drop StdErr Alpha LowerCL UpperCL tValue Probt;
    
29. run;
    
30. 
    
31. %end;
    
32. %mend NS;

复制代码

我的想法是每次将估计值esitmate&i利用call symput赋予parameter值，带入function计算

1. data _null_;
   
2. set parmds;
   
3. call symput(parameter, estimate1);
   
4. run;
   
5. 
   
6. 
   
7. proc fcmp outlib = work.func.nelson;
   
8. function nelson_siegel(years);
   
9. return(  
   
10. &b0+&b1*(1-exp(-years/&tao1))/(years/&tao1)+&b2*(1-exp(-years/&tao1))/((years/&tao1)-exp(-years/&tao1))+&b3*(1-exp(-years/&tao2))/((years/&tao2)-exp(-years/&tao2))
    
11. );
    
12. endsub;
    
13. run;
    
14. options cmplib=work.func;

复制代码

可是call symput只能使用一次，在第二次赋值时

1. data _null_;
   
2. set parmds;
   
3. call symput(parameter, esitmate2);
   
4. run;
   
5. 
   
6. 
   
7. data  rate2;
   
8. do years =1 to 30;
   
9. predyield= nelson_siegel(years);
   
10. output;
    
11. end;
    
12. run;

复制代码

data rate2 的结果仍然是第一次的结果，纠结了一下午无解。如能帮助，万分感谢！

许多的细节都衔接不上。如果只是要预测值，我想，程序拼凑起来很多。或许可以简化成一两行。京剧

jingju11 发表于 2013-11-6 10:39
许多的细节都衔接不上。如果只是要预测值，我想，程序拼凑起来很多。或许可以简化成一两行。京剧

谢谢，问题已经解决了