问题标签 [loess]

我有几个使用类似数据的黄土模型......

我想要完成的（没有任何运气）是确定平滑线在哪里做两件事：

1）我想确定线与 y=0 相交的 x 值是多少

2）我想确定 2 条黄土线在 x 的什么值处相互交叉。

我一直在寻找类似的问题和这些问题的解决方案太久了，但没有成功。任何帮助将不胜感激。

我想更改示例图中的 span 参数，示例图

在我为得到这个图而编写的代码中，除了默认的 0.2 之外，我无法让跨度做任何事情。我有一种预感，这与将黄土适当地分配给每个组有关，但我尝试过的方法没有奏效。我制作了一些示例代码来演示 span 参数不会影响下面的 loess fit。

尝试编写 LOESS 的 python/r 版本时提出的一些问题。我在网上搜索但没有得到满意的答案。Cleveland 编写的原始 C 版本在他的论文 Cleveland (1979, 1988) 中提出和描述可以在 R 源代码中找到，但是我不习惯 C 代码，而程序有点复杂，这使得它更加困难。

大多数其他语言的实现都没有原始 LOESS 的全部功能，例如 Cleveland (1988) 中描述的使用插值函数的预测函数。并且使用的 kd 树方法（在那篇论文中提到）在这些版本中丢失了。
所以，我的问题是：

1.

在单元顶点处评估黄土拟合，然后混合函数进行插值

（来自 R 中的 LOESS 手册）。

那么如果我写一个python/r版本，模型应该怎么存储呢？每个顶点存储一个拟合值、混合函数、边（窗口、单元格或其他）？

2、如何实现LOESS模型的predict()函数？如何预测输入数据（首先查找它属于哪个窗口的数据，然后使用混合函数计算值）？

我有这些数据：

看起来像：

我想找到一种方法让肘部/膝盖点在附近x=6.5

我认为拟合loess曲线然后取二阶导数可能会起作用，但是：

看起来它不会完成这项工作。

任何想法？

我正在尝试按组设置 proc 模板生成的线条的颜色。

我的数据有以下列：

天数 - 事件发生后的天数（x 轴）

测量 - 用于黄土曲线的数据点（y 轴）

MarkerColor - 我想用于线条的颜色

groupnm - 区分组以生成每一行

我geom_smooth从ggplot2包中使用在时间序列散点图上创建一条平滑线（一年中的每一天一个点，所以我有 365 个点）。其中一个参数称为span，进入帮助文件 ( ?geom_smooth) 会给出以下描述：

span控制默认黄土平滑器的平滑量。较小的数字产生更弯曲的线条，较大的数字产生更平滑的线条。

但是，这实际上并没有告诉我 span 参数控制的是什么。将其设置为 1 是没有用的，而将其设置为 0.1 提供了一些看起来不错的东西。

span = 0.5

span = 0.1

但是，在描述情节时，由于我不完全确定span实际发生了什么变化，因此我不确定如何描述平滑线。任何指针？

我正在尝试过滤掉太接近或低于黄土曲线的点：

结果如下所示：

显然不是想要的结果。

但是，如果我使用 scatter.smooth 函数，我会得到一个正确的曲线：

如何通过我的数据正确拟合黄土曲线？

给定如下所示的数据：

我正在尝试为这些数据拟合几种不同的方法。对于这个问题，我正在研究二阶多项式拟合与黄土拟合。为了获得更平滑的曲线，我想扩展 x 数据以给我更多的点来预测。所以对于我的黄土曲线，我这样做：

当我尝试用二阶多项式拟合做同样的事情时 - 我得到一个错误

显然Poly2和Loess_Fit有一些区别，但我不知道有什么区别。有没有办法像我使用 Loess_Fit 一样平滑 Poly2 贴合度？

在 R 中，使用 ggplot2 包，我有以下两行：

第一个为我的数据提供了一条带有阴影置信区间的未加权平滑曲线，但是当我添加权重时，尽管我声明了，但置信区间区域消失了se = TRUE。

这是我的数据

我正在比较使用 LOESS 回归的两条线。我想清楚地显示两条线的置信区间，我遇到了一些困难。

我尝试过使用各种线型和颜色，但在我看来结果仍然是忙碌和混乱。我认为置信区间之间的阴影可能会使事情变得更清晰，但是考虑到我的编码到目前为止的结构，我在解决这个问题时遇到了一些困难。我已经包含了生成的图、Analysis5k 和 Analysis5kz 两组的数据，以及到目前为止的代码。

我已经看到了一些示例，其中两个多边形重叠以显示置信区间重叠的位置，这似乎是呈现数据的好方法。如果有一种方法可以在两个置信区间共享的区域中绘制多边形，那可能是呈现数据的另一种好方法。

我了解应该如何完成多边形的基本概念，但我发现的示例已应用于更简单的线条和数据。到目前为止，对于一些糟糕的组织来说，部分原因是我自己的错，但由于这一步基本上是我数据呈现的最后润色，我真的不想从头开始重新工作。

非常感谢任何帮助或见解。

更新

我更新了标题。我收到了一些使用 ggplot 的好例子，虽然我想在未来使用 ggplot，但到目前为止我只处理了基础 R。对于这个特定的项目，如果可能的话，想尝试将其保留在基础 R 中。

分析5k

分析5kz

代码