作者Agingu (小土)
站内Statistics
标题[问题] 请教利用SAS作Non-linear regression的语法？
时间Wed May 2 14:55:46 2007
大家好！我是化学所研究生，最近想利用实验得到的已知数据（X,Y,Z），经由Gauss-Newton regression去fit一个反应model：
Z=(a*b*c*X*Y)/((1+b*X)*(1+c*Y))
想得到a,b,c三个参数。但在语法上碰到困难，希望先进们不啬指教。
经由实验可以得到三个column的数据，其中X(PCO),Y(PO2),Z(RATE)
利用下列指令：
data LH;
infile 'c:\data80.csv' delimiter=',' dsd;
input PCO PO2 RATE ;
run;
proc print data = LH;
run;
可以得到
The SAS System
2007年04月25日 星期三 下午09时47分25秒 6
Obs PCO PO2 RATE
1 1.000 0.370 .0000210
2 1.000 0.563 .0000393
3 1.000 0.709 .0000504
4 1.000 0.822 .0000568
5 1.000 2.498 .0001390
6 1.000 4.986 .0001630
7 1.000 7.900 .0001720
8 1.000 1.000 .0000721
9 0.194 1.000 .0000148
10 0.434 1.000 .0000252
11 0.663 1.000 .0000493
12 0.915 1.000 .0000610
13 0.256 1.000 .0000218
14 0.470 1.000 .0000318
15 0.842 1.000 .0000650
16 1.773 1.000 .0001390
17 2.941 1.000 .0002210
18 4.440 1.000 .0002460
接下来我利用”proc nlin”的指令，却失败，我相信我的语法有问题，恳请告知正确的语法！！非常感谢～～
proc nlin data=LH;
parameters a=1 b=1 c=1 ;
model Rate=(a*b*c*PCO*PO2/((1+b*PCO)*(1+c*PO2));
run;
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.154.122