对于我的数据的特定 x,如何提取黄土函数适合它使用的多项式函数的参数?
例如,在:
cars.lo <- loess(dist ~ speed, cars)
cars.lo
当 cars.lo$x == 5 时它适合什么?
更新:我想要多项式函数的参数,而不是黄土的预测(predict)。
我要求它估算该点的斜率。
谢谢, 塔尔
我将计算 x、x+eps、x-eps 处的预测值,然后对结果进行二次拟合。它在计算机时间上效率不高,但如果你不必经常这样做,它在程序员时间上是非常有效的。
我会说
predict(cars.lo, data.frame(speed=5))
[1] 7.797353
我认为您将独自一人。我会先阅读 中的参考资料lowess
,然后阅读 Clowess
函数的源代码。我建议从那里开始而不是从那里开始,loess
因为它是一种较旧且稍微简单的算法,您可能会发现它更容易适应您的需求。