R中GAM模型异常：分类变量首个类别未参与建模

这不是mgcv的bug，是R里因子变量建模的默认参照组机制：

• 因子的第一个水平会被自动设为「参照组」，模型里其他水平的系数都是和这个参照组的差值，所以参照组不会单独出现在系数表和绘图里——它的效应被当作基准（0值），所有其他类别都是相对于它的偏移量。

解决方法 #

1. 手动更换参照组 #

如果想让某个特定类别作为参照（比如把原来的非首个类别设为参照，让首个类别显示差值），用relevel()函数调整：

# 把Bottom_typef的参照组改成Seagrass
data$Bottom_typef <- relevel(data$Bottom_typef, ref = "Seagrass")
# 把Tidal_stagef的参照组改成Outgoing
data$Tidal_stagef <- relevel(data$Tidal_stagef, ref = "Outgoing")

# 重新拟合模型
test1 <- gam(Depth ~ Bottom_typef + Tidal_stagef, data = data, method = "REML")
summary(test1)
plot(test1, all.terms = TRUE, page = 1)

这样原来的首个类别（Artificial Reef、Incoming）就会以新参照组为基准，出现在系数表和绘图里。

2. 使用总和对比（显示所有类别系数） #

如果想让因子的所有水平都显示单独的系数（没有参照组，所有系数总和为0），可以修改因子的对比方式为contr.sum：

# 对单个因子设置总和对比
contrasts(data$Bottom_typef) <- contr.sum(nlevels(data$Bottom_typef))
contrasts(data$Tidal_stagef) <- contr.sum(nlevels(data$Tidal_stagef))

# 或者全局设置（会影响后续所有模型的因子对比）
options(contrasts = c("contr.sum", "contr.poly"))

# 重新拟合模型
test1 <- gam(Depth ~ Bottom_typef + Tidal_stagef, data = data, method = "REML")

此时summary里会列出每个因子的所有水平系数，绘图也会展示每个类别的效应值。

3. 可视化所有类别的实际预测值 #

如果只是想直观看到每个类别对应的Depth预测值，而不是相对差值，可以用emmeans包计算并绘图：

library(emmeans)
# 计算每个因子组合的预测均值
emm <- emmeans(test1, ~ Bottom_typef + Tidal_stagef)
# 可视化
plot(emm)

这样能直接看到每个类别对应的实际Depth估计值，更直观。

内容的提问来源于stack exchange，提问作者Zesra