R loess:全局“参数”项的系数

Question

有没有一种方法可以在局部回归建模中提取全局拟合项的系数？

也许我确实误解了like函数中全局拟合项的作用，但我希望有以下几点：

# 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的值)？

Answer 1

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

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