SAS中文论坛

标题: 求助：请高手帮俺写个循环的语句 [打印本页]

作者: shiyiming 时间: 2010-12-7 08:16
标题: 求助：请高手帮俺写个循环的语句
因为要求目标函数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 时间: 2010-12-9 03:29
标题: Re: 求助：请高手帮俺写个循环的语句
you may like try Proc IML.

作者: shiyiming 时间: 2010-12-9 05:28
标题: Re: 求助：请高手帮俺写个循环的语句
[quote:3p1qs3gz]如果做这个循环，能不能在循环中嵌套proc nlp 得到每一步的最优解呢？ [/quote:3p1qs3gz]

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

作者: shiyiming 时间: 2010-12-9 09:31
标题: Re: 求助：请高手帮俺写个循环的语句
他的意思是要人帮忙写个宏循环。

作者: Qiong 时间: 2010-12-10 09:27
标题: Re: 求助：请高手帮俺写个循环的语句
to sxlion
顺便问问宏变量怎么弄出小数啊？
比如
%let a=0.1;
怎么让b=2*a呀?不用data _null_~~

作者: shiyiming 时间: 2010-12-10 12:48
标题: Re: 求助：请高手帮俺写个循环的语句
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 时间: 2010-12-11 05:07
标题: Re: 求助：请高手帮俺写个循环的语句
非常感谢各位高手：） <img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" />

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)