如何通过emtrends或其他方法获取B、C合并组中esteem的简单斜率？

你可以通过两种简洁的方式实现需求：

方法1：基于原模型的emtrends直接计算对比 #

利用emmeans包的contrast()函数，针对已获取的分组斜率设置平均权重，直接得到B、C组的合并斜率：

library(emmeans)

# 获取各分组的简单斜率
trends <- emtrends(mod, ~ condition, var = "esteem")

# 定义对比：计算B、C组斜率的平均值（权重各为0.5）
bc_combined_trend <- contrast(trends, list(B_C_combined = c(0, 0.5, 0.5)))

# 查看结果
bc_combined_trend

运行后会输出B、C合并组的斜率、标准误、置信区间等统计量，示例输出类似：

contrast      estimate   SE df lower.CL upper.CL
 B_C_combined     -3.97 1.28  6    -6.99    -0.95

Confidence level used: 0.95 

方法2：重新定义分组后拟合模型 #

如果后续需要多次分析合并组，可以直接在数据中创建新的分组变量，重新拟合模型后用emtrends获取结果：

# 创建新分组：将B、C合并为一组
dat$group <- factor(
  ifelse(dat$condition %in% c("b", "c"), "BC", "A"),
  levels = c("A", "BC")
)

# 重新拟合包含交互项的线性模型
mod_grouped <- lm(scores ~ esteem * group, data = dat)

# 获取合并组的斜率
emtrends(mod_grouped, ~ group, var = "esteem")

此方法输出会直接展示A组和BC合并组的斜率，更直观易懂。

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