SAS中文论坛

标题: 请教多元非线性回归 [打印本页]

作者: shiyiming 时间: 2010-6-9 21:00
标题: 请教多元非线性回归
请教高手我在进行多元非线性回归时只迭代一两步就终止，显示 WARNING: Step size shows no improvement.这是怎么回事？
DATA FN;
INPUT X1 X2 X3 Y;
CARDS;
1.155 81.691 0.207 199.882
1.155 83.378 0.159 225.491
1.155 84.970 0.198 185.626
1.155 92.234 0.170 233.459
1.154 97.989 0.200 206.253
1.155 85.039 0.199 192.722
1.155 90.748 0.197 185.617
1.154 94.948 0.222 168.092
1.155 96.539 0.234 153.488
1.154 99.311 0.195 193.969
1.144 87.947 0.179 207.582
1.155 90.961 0.199 176.422
1.145 92.252 0.207 183.556
1.155 92.419 0.182 188.362
1.154 92.864 0.199 188.878
1.155 85.447 0.209 195.508
1.155 85.847 0.222 182.129
1.154 89.537 0.218 178.722
1.149 90.722 0.236 171.609
1.154 91.093 0.216 177.620
1.148 91.775 0.205 174.295
;
PROC NLIN;
PARMS A0=130 TO 150 BY 0.01
   A1=-3.00
   A2=3.50
   A3=0.20
   A4=-0.18
   A5=0.50
   A6=1.50;
MODEL Y=2/SQRT(3)*A0*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4);
DER.A0=2/SQRT(3)*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4);
DER.A1=2/SQRT(3)*A0*X1*EXP(A1*X1+A2)*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4);
DER.A2=2/SQRT(3)*A0*EXP(A1*X1+A2)*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4);
DER.A3=2/SQRT(3)*A0*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4)*(X2/10)**(A3*X1+A4)*LOG(X2/10)*X1;
DER.A4=2/SQRT(3)*A0*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4)*(X2/10)**(A3*X1+A4)*LOG(X2/10);
DER.A5=2/SQRT(3)*A0*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*(A6*(X3/0.4)**A5-(A6-1)*X3/0.4)*(X3/0.4)**A5*LOG(X3/0.4)*A6;
DER.A6=2/SQRT(3)*A0*EXP(A1*X1+A2)*(X2/10)**(A3*X1+A4)*((X3/0.4)**A5-X3/0.4);

RUN;
求高手解答！！

作者: shiyiming 时间: 2010-6-11 09:17
标题: Re: 请教多元非线性回归
nlin 过程默认的采用 Gauss迭代求解方法。
PROC NLIN method=NEWTON; 改成newton方法，就可以收敛了；
如果试遍所有地方法，都不能收敛，就只能增加观测值了

作者: shiyiming 时间: 2010-6-11 14:23
标题: Re: 请教多元非线性回归
有时候，对A1-A6赋初值，很重要。先要有个预估值，不要与实际值差距太大，太大也经常导致不收敛。

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