在 r 中报告线性模型中的交互

Question

我正在尝试报告线性模型中的所有交互：

mod1.lme <- lm(volume ~ Group * Treatment + Group + Treatment, data = df)

组是一个因子变量，具有 3 个水平：A、B 和 C。

我目前得到的结果是（我自己编的）：

这两个估计值参考了治疗：A，但我想单独查看每个效果。所以我想得到的输出是：

Treatment:A
Treatment:B
Treatment:C

如果我消除在末尾添加 -1 的截距，我会得到：

最好的编码方式是什么？

谢谢

Answer 1

您看到输出的原因是Treatment 的因子水平之一成为参考水平。在解释模型时，系数变为“与参考水平的效果差异”。只要模型包含截距，这就是必要的，因此获得所有显示系数的解释的唯一方法是删除截距，如下所示。

mod1.lme <- lm(volume ~ Group * Treatment - 1, data = df)

编辑：

要更改交互效果的名称，必须手动编辑名称

sum.lm <- summary(mod1.lme)
rownames(sum.lm$coef) <- c("groupA","groupB","groupC", "groupA:Treatment", "groupB:Treatment", "groupC:Treatment")

或者使用另一个包作为摘要，例如sjPlot

library(sjPlot)
tab_model(mod1.lme, pred.labels = c("groupA","groupB","groupC", "groupA:Treatment", "groupB:Treatment", "groupC:Treatment"))