如何在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
我的问题是我们如何找到这个模型的其他参数 a4
和 a5
(常量)在显示的等式中 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])
我想计算 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
我的问题是我们如何找到这个模型的其他参数 a4
和 a5
(常量)在显示的等式中 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])