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

我试图在一个线性模型中报告所有的相互作用，该模型如下：

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

如果我去掉截距，在末尾加-1，我得到：

编写此代码的最佳方法是什么

---

感谢

您之所以能看到您现在的输出，是因为

治疗

的一个因子水平成为参考水平。解释模型时，系数变为“与参考水平的有效差异”。只要模型包含一个截距，这是必要的，因此获得所有系数解释的唯一方法是删除截距，如下所示

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

---

非常感谢。我知道消除截距可以避免引用，但只会对整个组产生影响。我的兴趣是互动。A:治疗，B:治疗，C:治疗。这有意义吗？能做到吗？嗨，莉莉，你说的话听起来就像是删除了截获。您所寻找的效果与

summary（mod1.lme）

删除截距后显示的效果有何不同？为了更准确，我对文本进行了编辑：）thanksAh很有意义。您希望将“治疗”改为“治疗：A组”，以明确其仍受研究组的影响。在这里，你会撞到一堵墙。这是完全可能的，但格式化系数名称不是

帮助（summary.lm）

中的选项，也不是它的打印方法（同一页）。最简单的方法是将摘要存储在变量中，然后更改名称，或者使用提供更多选项的包（例如）。太棒了，谢谢！

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")

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