R 我是否可以将SSlogis用作nls方程的一部分,以及如何找到SSlogis的系数
我正在尝试一个R 我是否可以将SSlogis用作nls方程的一部分,以及如何找到SSlogis的系数,r,nls,coefficients,R,Nls,Coefficients,我正在尝试一个nls方程来预测土壤呼吸速率(SR)。我使用土壤温度、植被指数和土壤含水量作为独立因子。我假设温度(T10)、植被指数(VI)和土壤含水量(VWC)分别与R呈指数、线性和逻辑相关。所以我的R代码是 R1<-nls(SR ~ I(Rref*exp(log(Q)*T10) + a*VI+b+SSlogis(VWC,5,12,1)), data=SR, start=list(Rref=0.82, Q=1.5, a=10, b=-1),
nls
方程来预测土壤呼吸速率(SR)。我使用土壤温度、植被指数和土壤含水量作为独立因子。我假设温度(T10)、植被指数(VI)和土壤含水量(VWC)分别与R呈指数、线性和逻辑相关。所以我的R代码是
R1<-nls(SR ~ I(Rref*exp(log(Q)*T10) + a*VI+b+SSlogis(VWC,5,12,1)),
data=SR, start=list(Rref=0.82, Q=1.5, a=10, b=-1),
trace=T)
如何获取
SSlogis
中的系数值,即Asym、xmid和scal?所示代码不使用5,12,1作为起始值。它将它们固定在这些值上。感谢您的推荐。那么如何正确地重写代码呢?我尝试代码``r2,然后运行`r2尝试写出Asym/(1+exp((xmid输入)/scal)),而不是使用SSlogis。显示的代码不使用5,12,1作为起始值。它将它们固定在这些值上。感谢您的推荐。那么如何正确地重写代码呢?我尝试代码``r2,然后运行`r2尝试写出Asym/(1+exp((xmid输入)/scal)),而不是使用SSlogis。
Parameters:
Estimate Std. Error t value Pr(>|t|)
Rref 1.71770 3.11436 0.552 0.5836
Q 1.06064 0.05253 20.191 <2e-16 ***
a 8.47050 3.45392 2.452 0.0176 *
b -8.48737 3.50457 -2.422 0.0190 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 1.5 on 52 degrees of freedom
Number of iterations to convergence: 46
Achieved convergence tolerance: 8.868e-06