参考直线化可以;但随便试也可以,只要每个参数有意义就行,只要有意义,不管初始值是什么,估计值都相同。(建议没意义的初值一个量级一个量级的试)
这是《医学统计学》上的例子,具体解释可以参看我上传 proc nlin的附件(没有论坛币的同学私信我,我邮箱发)。
- /*例子12.1*/
- data eg12_1 ;
- input x y @@ ;
- cards ;
- 60 0.8 80 6.5 100 20.5 120 45.9
- ;
- run;
- symbol value = dot cv = red
- interpol = spline ci =green;
- proc gplot data = eg12_1;
- plot y * x /
- vaxis =0 to 50 by 10
- haxis =60 to 120 by 20;
- run;
- quit ;
- /*直线化*/
- data eg12_1;
- set eg12_1;
- xx = x/20-3;
- yy = log(y+5);
- run;
- symbol value = dot cv = red
- interpol = spline ci =green;
- proc gplot data = eg12_1;
- plot yy * xx ;
- run;
- quit ;
- /*12.1 nlin*/
- proc nlin data = eg12_1;
- parms a=-0.19 b=0.01;
- model y = exp(a+b*x);
- output out= neg12_1 predicted =yp;
- run;
- symbol value = dot cv = red
- interpol = spline ci =green;
- proc gplot data = neg12_1;
- plot yp * x ;
- run;
- quit ;
复制代码