R中glm与sjPlot绘图的Odds Ratio结果不一致问题求助

核心矛盾根源：变量编码/参考组不一致 #

这是逻辑回归OR值可视化与模型输出矛盾的最常见原因，具体分两类：

• 因变量（MAE）的事件方向反转
  逻辑回归中，glm默认将因变量的最后一个水平作为拟合的事件（即计算P(Y=最后一个水平)的概率）。如果你在glm中实际关注的是“MAE高”的概率，但plot_model默认采用了相反的事件方向（比如拟合P(MAE=低)），那么OR值就会变成模型计算值的倒数——原本>1的结果就会显示为<1。

• 自变量（Plt）的参考组设置差异
  若Plt是分类变量：glm默认以因子的第一个水平作为参考组，如果你在绘图前无意间修改了Plt的水平顺序，或者plot_model被配置了不同的参考组，也会导致OR的解释方向反转。

验证与修正步骤 #

1. 核对因变量编码

   • 先查看因变量的水平顺序：table(your_data$MAE)，再看summary(model)的输出，确认glm拟合的是哪个事件的概率。
   • 显式指定事件方向后重新拟合并绘图：

     # 将MAE的目标事件（比如"high"）设为对应水平，统一模型与绘图的事件方向
     your_data$MAE <- factor(your_data$MAE, levels = c("low", "high"))
     model <- glm(MAE ~ Plt, data = your_data, family = binomial)
     # 绘图时显示数值方便核对
     plot_model(model, type = "odds", show.values = TRUE)
     

2. 检查自变量参考组

   • 若Plt是分类变量，运行levels(your_data$Plt)确认参考组，对比summary(model)中系数对应的对比组，再看绘图的x轴标签是否一致。
   • 统一参考组后重新建模：

     your_data$Plt <- relevel(factor(your_data$Plt), ref = "你的参考组水平")
     model <- glm(MAE ~ Plt, data = your_data, family = binomial)
     plot_model(model, type = "odds")
     

3. 手动验证OR值
   直接计算模型的OR值：exp(coef(model))，将结果与plot_model图中的数值对比。如果两者互为倒数，可100%确认是事件方向或参考组的问题。

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