作者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