如何計算“ 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))