你好,游客 登录
背景:
阅读新闻

R语言非线性最小二乘法

[日期:2017-09-12] 来源:yiibai  作者:初生不惑 [字体: ]

当对真实世界数据建模进行回归分析时,我们观察到模型的方程很少是给出线性图的线性方程。 大多数情况下,现实世界数据模型的方程式涉及更高程度的数学函数,如3sin函数的指数。 在这种情况下,模型的曲线给出了曲线而不是线性。线性和非线性回归的目标是调整模型参数的值以找到最接近您的数据的线或曲线。当找到这些值时,我们才能够准确估计响应变量。

在最小二乘回归中,我们建立了一个回归模型,不同点与回归曲线的垂直距离的平方和之和最小化。 我们通常从定义的模型开始,并假设系数的一些值。 然后应用R中的nls()函数来获得更准确的值以及置信区间。

语法

在R中创建非线性最小二乘检验的基本语法是 -

nls(formula, data, start)
R

以下是使用的参数的描述 -

  • formula - 是包含变量和参数的非线性模型公式。
  • data - 是用于评估(计算)公式中的变量的数据帧。
  • start - 是起始估计的命名列表或命名数字向量。

示例

我们将考虑一个假设其系数的初始值的非线性模型。 接下来,我们将看到这些假设值的置信区间是多少,以便可以判断这些值是如何进入模型的。

所以考虑下面这个方程式 -

a = b1*x^2+b2
R

我们假设初始系数为13,并将这些值拟合成nls()函数。

setwd("F:/worksp/R")
xvalues <- c(1.6,2.1,2,2.23,3.71,3.25,3.4,3.86,1.19,2.21)
yvalues <- c(5.19,7.43,6.94,8.11,18.75,14.88,16.06,19.12,3.21,7.58)

# Give the chart file a name.
png(file = "nls.png")


# Plot these values.
plot(xvalues,yvalues)


# Take the assumed values and fit into the model.
model <- nls(yvalues ~ b1*xvalues^2+b2,start = list(b1 = 1,b2 = 3))

# Plot the chart with new data by fitting it to a prediction from 100 data points.
new.data <- data.frame(xvalues = seq(min(xvalues),max(xvalues),len = 100))
lines(new.data$xvalues,predict(model,newdata = new.data))

# Save the file.
dev.off()

# Get the sum of the squared residuals.
print(sum(resid(model)^2))

# Get the confidence intervals on the chosen values of the coefficients.
print(confint(model))
R

当我们执行上述代码时,会产生以下结果 -

[1] 1.081935
Waiting for profiling to be done...
       2.5%    97.5%
b1 1.137708 1.253135
b2 1.497364 2.496484
R

输出的图形如下所示 -

我们可以得出结论,b1的值更接近于1,而b2的值更接近于2而不是3

收藏 推荐 打印 | 阅读: