寻找随机曲线拟合方法

机器算法验证 置信区间 随机过程 异方差 曲线拟合
2022-03-30 06:05:01

我正在寻找一种曲线拟合方法,它为拟合曲线提供上限和下限置信区间,如下图所示。假设红色实心曲线是可以通过普通曲线拟合工具得到的拟合曲线。然而,这样的曲线以某种方式预测了预期值。我正在寻找拟合曲线的上限和下限的虚线曲线。请注意,我的数据具有异方差性。任何帮助深表感谢。 在此处输入图像描述

1个回答

您可以在 R 中使用 ggplot2 库的 geom_smooth() 函数。

这是一个例子:

x = seq(0, 100, by=0.5)
y = sqrt(x)
y = y + rnorm(n = length(y),mean = 0,sd = 3)
df = data.frame(cbind(x,y))
require(ggplot2)
ggplot(data = df, aes(x = x,y = y)) + geom_smooth()

这是输出:

在此处输入图像描述

geom_smooth 调用了一种叫做 loess 的曲线拟合方法,它进行局部回归。使用默认选项,对于每个 x,它会考虑一个包含 75% 数据点的邻域,并使用加权最小二乘法拟合二次方。它假设错误是正态分布的,并计算置信区间,如http://www.netlib.org/a/cloess.pdf第 44-46 页中所述

这是黄土的文档: https ://stat.ethz.ch/R-manual/R-devel/library/stats/html/loess.html