使用allEffects生成的图表数值错误，有什么解决方法及替代绘图方案？

allEffects问题修复方案 #

• 确认变量指定逻辑：调用allEffects()时需显式指定模型中的交互项、自变量取值范围，避免函数默认插值导致数值和预期不符，参考调用示例：

# 假设你的拟合模型为mod <- lm(因变量 ~ 性别 * 核心自变量, data = 数据集)
eff_obj <- allEffects(mod,
                      # 显式指定核心自变量的取值序列，和你预期的展示范围一致
                      xlevels = list(核心自变量 = seq(min(数据集$核心自变量), max(数据集$核心自变量), length.out = 100)),
                      # 保留性别分组
                      factors = list(性别 = unique(数据集$性别)))
# 提取效应结果转为数据框，可直接核对数值是否和预估一致
eff_df <- as.data.frame(eff_obj[[1]])

• 调整绘图参数避免分面：使用effects包自带绘图函数时，添加multiline = TRUE参数即可将两组折线合并到同一张图中，无需拆分：

plot(eff_obj,
     multiline = TRUE,
     confint = TRUE,
     main = "",
     xlab = "自定义X轴标签",
     ylab = "自定义Y轴标签")

无需allEffects的替代方案（效果可控，推荐毕设使用） #

直接基于模型预测值配合ggplot2绘图，所有元素可自定义，不会出现数值或样式不符的问题，实现步骤如下：

1. 构造预测数据集并计算预测值及置信区间

library(ggplot2)
# 生成覆盖所有展示维度的预测数据集
pred_df <- expand.grid(
  核心自变量 = seq(min(数据集$核心自变量), max(数据集$核心自变量), length.out = 100),
  性别 = unique(数据集$性别)
)
# 计算预测值及95%置信区间
pred_res <- predict(mod, newdata = pred_df, se.fit = TRUE)
pred_df$预测值 <- pred_res$fit
pred_df$置信区间下限 <- pred_res$fit - 1.96 * pred_res$se.fit
pred_df$置信区间上限 <- pred_res$fit + 1.96 * pred_res$se.fit

2. 绘制目标双折线图

ggplot(pred_df, aes(x = 核心自变量, y = 预测值, color = 性别, fill = 性别)) +
  # 绘制折线
  geom_line(linewidth = 1.2) +
  # 绘制置信区间阴影
  geom_ribbon(aes(ymin = 置信区间下限, ymax = 置信区间上限), alpha = 0.2, color = NA) +
  # 自定义标签
  labs(x = "X轴标签", y = "Y轴标签", color = "性别", fill = "性别") +
  # 调整主题
  theme_classic()

你当前拆分性别的出图参考：

你提供的预估数值参考：

上述代码可直接实现你需要的目标图表样式，可自行调整配色、坐标轴刻度、文字大小等参数匹配需求：

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