Question

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

也许我误解了全局拟合术语在功能黄土中的作用，但我想要的是以下内容：

# 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)

然而，这两种型号的配合并不完全相同（可能是由于不同的控制设置？）......

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

1 个答案: