如何在R中找到第四和第五回归系数?

How to find fourth and fifth regression coefficients in R?

我想计算 5 回归 coefficients.I 通过互联网搜索,但我没有找到任何相关信息。 我的数据:

    y=c(2,13,0.4,5,8,10,13)
    x=c(2,13,0.004,5,8,1,13)
    z=c(2,3,0.004,15,8,10,1)

正规方程:

   y=a1x+a2z+a3

x,z,自变量,y为因变量,a1 a3,a2为模型参数。

这个等式的正常拟合::

       fit=lm(y~x+z)

但我的等式看起来像:

      y=((a1x+a2z+a3)^a4)+a5

我的问题是我们如何找到这个模型的其他参数 a4a5(常量)在显示的等式中 above.what 我应该在 R 中使用它们来找到它们吗?

下面使用 optim 怎么样?

f <- function(p){
          sum((y - (p[1]*x+p[2]*z+p[3])^p[4]+p[5])^2)
     }

p <- optim(rep(.5, 5), f)$par
p
[1]  3.5539397  0.8423521  0.1872422  0.6287906 -0.1863242

所以,a5-0.1863242。拟合值如下所示:

plot(seq_along(y), y)
lines((p[1]*x+p[2]*z+p[3])^p[4]+p[5])