r loess：全局“参数”项的系数答案

【问题标题】：r loess: coefficients of global "parametric" termsr loess：全局“参数”项的系数
【发布时间】：2016-07-01 07:25:00
【问题描述】：

有没有办法在局部回归建模中提取全局拟合项的系数？

也许我确实误解了全局拟合项在函数 loess 中的作用，但我想要的是以下内容：

# baseline:
x <- sin(seq(0.2,0.6,length.out=100)*pi)

# noise:
x_noise <- rnorm(length(x),0,0.1)

# known structure:
x_1 <- sin(seq(5,20,length.out=100))

# signal:
y <- x + x_1*0.25 + x_noise

# fit loess model:
x_seq <- seq_along(x)
mod <- loess(y ~ x_seq + x_1,parametric="x_1")

拟合完成得很完美，但是，如何提取全局拟合项 x_1 的估计值（即上面示例中接近 0.25 的某个值）？

【问题讨论】：

标签： r coefficients loess

【解决方案1】：

最后，我使用 gam 包中的函数 gam 找到了解决问题的方法：

require(gam)
mod2 <- gam(y ~ lo(x_seq,span=0.75,degree=2) + x_1)

但是，两个模型的拟合并不完全相同（可能是由于控制设置不同？）...

【讨论】：