如何从R中的平滑内核估计函数中获取任何（x，y）点？

Question

我有一个散点图，我使用ksmooth函数来获得平滑线。 现在我想通过从平滑线中减去实数来获得残差。 也许有人知道如何为任何x找到y的值？或者可能还有另一种方法来获取残差？

plot(x,y)
kernel <- ksmooth(x,y, kernel="normal", bandwidth=0.01)
lines(kernel, col=2)

Here is the result

Answer 1

您可以直接从ksmooth返回的模型对象访问拟合值。这是一个MWE：

x <- 1:100
y <- rnorm(100, mean=(x/2000)^2)
plot(x,y)
kernel <- ksmooth(x,y, kernel="normal", bandwidth=10, x.points=x)
lines(kernel, col=2)
resid <- kernel$y - y
print(resid)

然后运行

all.equal(kernel$x, x)