你想变动哪个数, parm 吗? 我对MV模型不熟悉,不太清楚你的问题的背景,否则我可以帮你用宏实现的

不好意思  是我没有说清楚。

我想变动的参数是 ge 的约束条件，以及最后非线性后处理得到最优化的值。

如不是用宏处理的话，我要这样计算每一个值，比较烦琐：

data quad1;

set quad;

if _type_ = "ge" then _rhs_ =0.02 ；(上面的是0.019)

proc nlp inquad=quad1 outset=nlp_out7;

min;

parms x1-x8;

run;

然后又要从nlp-out7提取目标函数。

我的想法其实就是直接建立宏，或者其他可行的方法。一次计算所有的值，得到相应的目标值。

这个不难，首先你在宏定义外面建好你要参数表，ge行的_rhs_可以留空．

然后这样．

%macro main(start,end,increment);

   %let i=&start.;

   %do  %while(&i.<=&end.);

data temp;

 set quad;

 if _type_=’ge’ then _rhs_=&i.;

 run;

proc nlp inquad=quad outest=nlp_out;

 min;

 parms x1-x8;

run;

data function;  ge=&i.;

set nlp_out(where=(_type_=’OBJ’) keep=value);

 keep ge value; run;

value 就是函数值*/;

%if &i.=&start. %then %do;  data base; set function;run;                            

                  %end;

%else do;  proc append base=base data= function force;run;

%let I = %eval(&i.+&increment.);

%end;

%mend;

   %main(0.006,0.02,0.002);

谢谢~~主要

明天我把论坛币汇到你帐户谢谢啊~~~

data base;

set function;run; 与proc append base=base data= function force;run;

有什么作用，我感觉是不是重复了？？我对sas是正处于探索阶段，用它做毕业论文模型数据处理。

[此贴子已经被作者于2008-3-18 23:45:59编辑过]

不重复的，第一步是初始化，如果在base不存在的时候append，结果没错，但是log会出error message