用r编写非线性回归方程的正确方法是什么
a0 *(D ^ a1)*(H ^ a2) X ^(b1 (Z ^ 4)+ b2 *(1 / exp(D / H)+ b3 *(X ^ 0.1)+ b4 *(1 / D)+ b5 *(H ^ Q)+ b6 * X)
我尝试使用下面的格式编写它,但是r表示我应该删除公式和数据之间的逗号()。我尝试删除逗号符号,但仍然无法运行数据。
我的数据是带有D,H,Z,X,Q变量的表格格式。而a0,a1,a2,b1,b2,b3,b4,b5,b6都是要估计的参数。
model2 <- nls (Dbh ~ a0*(Dbh^a1)*(Tht^a2)*X^(b1*(Z^4)+b2*((1/exp)^(Dbh/Tht))+b3*(X^0.1)+b4*(1/Dbh)+b5*(Tht^Q)+b6*X)),
data = kodat,
Start = list(a0= 4.6, a1= 6.9, a2= 5.9, b1= 30, b2= 2.1, b3= 7.8, b4= 1.2, b5= 0.2, b6= 8.1)
summary(model2)