估算MGCV的特定价值

Question

我正在尝试运行GAM模型，并根据来自this thread的包mgcv的建议，在给定值下获得精确估计。使用预测函数时，结果与GAM的估计值不匹配。

这是我感兴趣的模型：

mod <- gam(y~s(a, by=b), data = dat)

当我使用以下方式获得结果时

plot(mod)

结果与我尝试使用以下方法估计变量a上的单个绘图点时不同：

pdat <- with(dat,
             data.frame(a = 0,
                        b = mean(b,na.rm = TRUE)))
pred <- predict(mod, pdat, type = "response", se.fit = TRUE)
pdat <- transform(pdat, fitted = pred$fit)
pdat <- transform(pdat, upper = fitted + (1.96 * pred$se.fit),
                  lower = fitted - (1.96 * pred$se.fit))

有没有人对我可能做错了什么有任何建议？

Answer 1

通过反复试验，我找到了问题的答案。它接缝，而不是值的平均值，因为我的变量我需要用1替换b - 这将获得图中的确切值。

也许这将有助于其他人。