SAS中文论坛

 找回密码
 立即注册

扫一扫,访问微社区

查看: 651|回复: 2
打印 上一主题 下一主题

请教多元非线性回归

[复制链接]

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
楼主
 楼主| 发表于 2010-6-9 21:00:01 | 只看该作者

请教多元非线性回归

请教高手 我在进行多元非线性回归时只迭代一两步就终止,显示 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;
求高手解答!!
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
沙发
 楼主| 发表于 2010-6-11 09:17:08 | 只看该作者

Re: 请教多元非线性回归

nlin 过程默认的采用 Gauss迭代求解方法。
PROC NLIN method=NEWTON; 改成newton方法,就可以收敛了;
如果试遍所有地方法,都不能收敛,就只能增加观测值了
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
板凳
 楼主| 发表于 2010-6-11 14:23:25 | 只看该作者

Re: 请教多元非线性回归

有时候,对A1-A6赋初值,很重要。先要有个预估值,不要与实际值差距太大,太大也经常导致不收敛。
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|SAS中文论坛  

GMT+8, 2026-2-4 05:01 , Processed in 0.090812 second(s), 22 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表