如何绘制具有对数转换响应的 GAM 模型

Question

我使用 mgcv 创建了一个 GAM 模型。我的响应变量是对数转换的。我想使用

gratia

绘制这个模型。但是，我想确保在使用绘图函数绘制 y 轴后对 y 轴进行幂运算。有什么办法可以做到这一点吗？预先感谢！

下面是我的代码的片段。

gam.mod1 <- gam(log(sum7) ~ tap_type.fac + aq_race___1.factor + cityfac +
                  aq_race___2.factor + aq_race___3.factor + aq_race___4.factor
                + aq_race___5.factor + aq_race___99.factor + s(aq_age),family= gaussian, 
                data=d2.male)
draw(gam.mod1, residuals=TRUE)

忘了提及我的响应变量是连续的，我对其进行了对数转换以改进诊断。

Answer 1

这不适用于

residuals = TRUE

元素（但您可以添加实际的响应数据），但是，您可以使用

smooth_estimates()

及其

draw()

方法：

m <- gam.mod1 # to save typing (am doing this on a tablet)

m |>
  smooth_estimates() |>
  add_confint() |>
  draw(fun = exp, constant = model_constant(m)) &
  # add raw data on
  geom_point(data = d2.male, aes(x = aq_age, y = sum7))