「plot.gam」置信區間是如何計算的？

Question

如果模型是使用mgcv安裝然後順利術語繪製，

m <- gam(y ~ s(x)) 
plot(m, shade = TRUE) 

，那麼你得到一個置信區間的曲線的曲線圖。我猜想，這些是逐點置信區間。他們是如何計算的？

我試着寫

object <- plot(m, shade = true) 
object[[1]]$fit +- 2*object[[1]]$se 

，以便提取使用標準誤差和2的倍數的上限和下限，但是當我繪製它，它看起來有點比繪製的置信區間不同通過plot.gam？

那麼，這些如何計算？

我不使用seWithMean = true或類似的東西。

Answer 1

它是1的標準差。

oo <- plot.gam(m) 
oo <- oo[[1]] 
points(oo$x, oo$fit, pch = 20) 
points(oo$x, oo$fit - oo$se, pch = 20) 

重複的例子：

x <- seq(0, 2 * pi, length = 100) 
y <- x * sin(x) + rnorm(100, 0, 0.5) 
m <- gam(y ~ s(x))