问题标签 [loess]

我正在尝试使用内核或黄土平滑方法来平滑我的数据集。但是，他们都不清楚或不是我想要的。几个问题如下。我的 x 数据是“conc”，y 数据是“depth”，即 ex。厘米。

1）内核平滑

k <- kernel("daniell", 150)
plot(k)

K <- kernapply(conc, k)

plot(conc~depth) lines(K, col = "red")

在这里，我的数据按频率 = 150 进行了平滑处理。这意味着每个数据点都由相邻（左右）150 个数据点平均？“丹尼尔”是什么意思？我在网上找不到它的意思。

2）黄土光滑

p<-qplot(depth, conc, data=total)

p1 <- p + geom_smooth(method = "loess", size = 1, level=0.95)

这里，黄土平滑函数的默认值是多少？如果我想像上述情况一样以频率 = 150 平滑我的数据（每 150 个数据点移动平均值），我该如何修改此代码？

3) 为了用对数刻度显示 y 轴，我输入了“log10(conc)”，而不是“conc”，它起作用了。但是，我无法更改 y 轴刻度标签。我尝试在我的代码中使用“scale_y_log10(limits = c(1,1e3))”来显示轴刻度标签，如 10^0、10^1、10^2...，但没有用。请回答我的问题。非常感谢你的帮助。

和

我使用以下代码在 R 中完成了局部线性回归和内核回归：

有没有办法在 Matlab 中获得相同的结果？

使用函数 fit 和 fittype lowess，我接近 R 中局部线性回归 ( degree=1) 的结果。但是，它们并不相等。degree = 0我为具有两个自变量的内核回归 ( ) 找到的最佳命令是ksrmv，但是在输入时x1，x2我无法手动设置带宽/跨度。

谢谢你的帮助！

我ggplot2用来对我的数据进行平滑估计

它工作正常。现在，当尝试loess直接重现这个 using 函数时，我收到一个错误：

我认为该ggplot2版本有效，因为它设置了一些默认参数。

这是我的数据 yy ：

我有一个奇怪的问题，但希望有人能帮助我。我正在尝试创建一个湖底的表面图，然后添加一些显示植物频率的点，以便直观地了解整个湖中水生植物的出现位置。

现在，我正在分别使用 R 中的 scatterplot3d 和 lattice 包在 scatterplot3d 和线框中创建曲面图。为了实现我感兴趣的绘图类型，我已将深度转换为负值（想象湖的水面在z轴上为0），然后通过纬度和经度坐标创建了一个深度的黄土模型。但是，我遇到的一个问题是黄土模型预测的深度为正值（这在湖中当然是不可能的；只能从 0 深度进入水柱）。

例子

正如你所看到的，我的深度从 -14 到几乎 7。有没有办法在黄土模型上设置一个上限，这样我的模型就不会达到这些正值？

感谢您的帮助，
保罗

我试图在具有连续预测变量的二元结果变量的散点图上拟合一条黄土线。

这是我正在使用的代码：

这是我得到的图表：！https://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loessgraph1.png?at=master

我还尝试订购 x 变量：

并得到这张图：！https://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loessgraph2%28ordered_x%29.png?at=master

这是数据：https ://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loesscurvedata.txt?at=master

任何帮助深表感谢。

如何添加将另一列视为权重的黄土平滑？

假设我有以下data.frame：

绘制一个 loess-smoothing 会产生以下结果：

但我想使用列权重，使得 loess 与权重为 10 的每个点都存在 10 倍一样：

我找到了函数loess的参数权重。但我不知道如何从 geom_smooth() 调用它。

我正在尝试重新创建用于从 ggplot2 包中的 stat_smooth 估计平滑函数的过程。举个例子：

平滑线并不真正适合数据，但没关系。现在，让我们从头开始重新创建相同的图表。根据http://www.inside-r.org/r-doc/stats/loess我们需要使用三次加权内核和 2 次多项式（默认情况下）。我发现这篇http://www.maths.manchester.ac.uk/~peterf/MATH38011/NPR%20local%20Linear%20Estimator.pdf文章描述了如何估计光滑的黄土函数。我尝试重新创建此方法并将其用于我的数据：

我不确定参数h应该使用什么，但根据一本书，我有“对于具有公制宽度的三立方内核，h是支撑区域的半径”。因此，我尝试的第一件事是：

这显然不是同一个功能。我一直在使用不同的h值，但无法重新创建相同的曲线，这让我相信我在其他地方犯了错误。

如果之前已经回答过这个问题，我深表歉意，但我找不到适用的例子。

我正在尝试去除一些数据以进行变异函数分析。

我有一个数据框“aa”，其中包含“y”、“long”、“lat”和“z”列。

我正在尝试运行： loess(aa2$y ~ aa$long + aa$lat, aa, degree =2)在因子 z 的每个级别上。

最后，我需要一个“Long”、“Lat”、“Residual”和“Z”的数据框，残差来自多个特定因子的黄土对象。

鉴于我对 R 的了解有限，我无法找出正确的语法来实现这一点。

我假设可以使用 *apply 函数之一，但我不太了解该语言以正确编写它。

感谢您的任何指导或澄清。

我希望能够创建一个仅适合时间序列中本地 PRECEDING 点的黄土模型，因此它不使用来自未来的信息。

我想知道它是否受支持，或者是否有任何其他功能可以执行此操作。我试过使用子集参数，但它似乎是静态的，不适合我想要的。

我有下面的图表，我想通过平滑图表使其更悦目。可能吗 ？

理想图。