我正在使用指数非线性回归来拟合我的数据。
这是我的模特:
O-acetyl content = A*exp(-B*Time)+C
with:
- A: amplitude of decrease (difference between the intercept T0 and the plateau);
- B: speed of decrease (curvature);
- C: low plateau;
- O-acetyl content: expressed in µmol/dose;
- Time: expressed in hours.
我必须将我的因变量(O-乙酰基含量)传递给ln。我已经获得了系数,但是如何解释结果呢?我正在使用系数估计的Newthon方法。
例如,这些结果用于参数估计:
A = 0.062
B=0.0573
C=0.0309
我的SAS脚本(一部分):
proc nlin data=donn method=newton;
parms A=0.8815 B=0.0124 C=4.4067;
if (CCR="before_change") then do;
beta_1 = A1;
beta_2 = B1;
beta_0 = C1;
end;
else do;
beta_1 = A2;
beta_2 = B2;
beta_0 = C2;
end;
model l_AgHA=beta_0 + beta_1*exp(-beta_2*hours); *variable dependante en ln;
output out=nlinexp_full predicted=pred l95m=l95mean u95m=u95mean
l95=l95ind u95=u95ind;
ods output ParameterEstimates=pest_full;
ods output Anova=aovred_full(rename=(ss=ssred ms=msred df=dfred));
run;
非常感谢。