求助：请高手帮俺写个循环的语句

shiyiming

因为要求目标函数y=P（概率）×log（w×Y+（1-w））的最大值，我利用nlp找了一个解决的方法：
proc nlp;
  max y;
PARMS w sigma;
bounds 0<=w<=1, 0<=sigma;
Y1= exp(sigma*0);
Y2= exp(sigma);
Y3= exp(-sigma);
P1= 1/3;
P2= 1/3;
P3= 1/3;
y=P1*log(w1*Y1+(1-w))+P2*log(w*Y2+(1-w))+P3*log(w*Y3+(1-w));
run;
得到的结果
w=0.5， sigma= 709.08955 ， max y=235.90


我想在想做一个循环，
首先， w 从 0 逐步增加到 1， 步长为 0.1
然后 在 每一个给定的w，  sigma 从 0 逐步增加到710， 步长为1.
最后， 把每一个 w，sigma 相对应的 y值  导出来， 要求包括 （w， sigma， y）。
因为很久没有用sas了，循环语句很生疏了，请高手帮忙，谢谢在先了。

另外，如果做这个循环，能不能在循环中嵌套proc nlp 得到 每一步的最优解呢？

shiyiming

you may like try  Proc IML.

shiyiming

[quote:3p1qs3gz]如果做这个循环，能不能在循环中嵌套proc nlp 得到 每一步的最优解呢？ [/quote:3p1qs3gz]

给定w 和sigma，最优解是什么意思呢？

shiyiming

他的意思是要人帮忙写个宏循环。

Qiong

to sxlion
顺便问问宏变量怎么弄出小数啊？
比如
%let a=0.1;
怎么让b=2*a呀?不用data _null_~~

shiyiming

to edle
for the looping, maybe you want this? [manually compute Y for each given w, sigma]?

data comp_out;
       do w=0 to 1 by 0.1;
             do sigma=0 to 710 by 1;
             Y1= exp(sigma*0);
             Y2= exp(sigma);
             Y3= exp(-sigma);
            /* P1= 1/3; P2= 1/3; P3= 1/3;*/
             y=(log(w1*Y1+(1-w))+log(w*Y2+(1-w))+log(w*Y3+(1-w)))/3;
             output;  keep w sigma Y;
             end;
        end;
run;

don't understand ur second Q. Do u mean initial values?

shiyiming

非常感谢各位高手：） <!-- s:D --><img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" /><!-- s:D -->