如何在sjPlot/ggplot中设置LMER模型固定效应图的X轴范围与参考线

错误原因分析 #

你遇到的「离散值传入连续刻度」错误，是因为sjPlot默认生成的固定效应图（比如type="est"）中，若误操作了离散轴（比如对应自变量名称的y轴），或使用了非效应值为x轴的图类型（比如type="pred"），就会触发该报错。另外geom_vline()无效果，大概率是图层顺序问题，或是轴映射搞反了。

方案一：基于sjPlot现有图修改 #

sjPlot的plot_model()返回的是ggplot对象，直接用ggplot语法调整即可，核心是针对连续的效应值轴（默认是x轴）设置范围，并用oob_squish截断超出部分：

library(lme4)
library(sjPlot)
library(ggplot2)
library(scales)

# 示例lmer模型
data(sleepstudy)
model <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)

# 生成固定效应森林图（x轴为连续效应值）
p <- plot_model(model, type = "est")

# 修改x轴范围并截断超出部分，添加x=0垂直线
p +
  scale_x_continuous(limits = c(-1, 2), oob = oob_squish) +
  geom_vline(xintercept = 0, linetype = "dashed", color = "darkred", size = 0.8)

方案二：改用ggplot手动绘制（完全自定义） #

如果sjPlot的默认设置限制太多，直接用broom.mixed提取模型结果，手动构建ggplot更灵活：

library(broom.mixed)
library(ggplot2)
library(scales)

# 提取固定效应及置信区间
tidy_fe <- tidy(model, effects = "fixed", conf.int = TRUE)

# 手动绘制森林图
ggplot(tidy_fe, aes(x = estimate, y = term)) +
  geom_point(size = 2) +
  geom_errorbarh(aes(xmin = conf.low, xmax = conf.high), height = 0.2) +
  # 设置x轴范围并截断超出值
  scale_x_continuous(limits = c(-1, 2), oob = oob_squish) +
  # 添加x=0垂直线
  geom_vline(xintercept = 0, linetype = "dashed", color = "darkred") +
  labs(x = "固定效应估计值", y = "自变量") +
  theme_bw()

注意事项 #

• 确认调整的是连续轴：type="est"的固定效应图中，x轴是连续的效应值，y轴是离散的自变量名称，别搞反轴的设置。
• oob_squish参数的作用是把超出x轴范围的点/误差线截断到边界，而非删除数据，符合你「截断超出部分」的需求。

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