我有一个由多种类型和两种处理方法组成的数据集(示例如下)。我使用 R 中的 lm() 函数来分析它,并使用 emmeans 进行事后测试。
m <- lm(Y ~ Type*Treatment)
Anova(m, type = 3)
Run Code Online (Sandbox Code Playgroud)
如何在 R 中设置正交对比?如果我跑
emmeans(m, ~Type*Treatment)
Run Code Online (Sandbox Code Playgroud)
它比较一切。我有兴趣比较类型中的不同治疗和治疗中的不同类型,但不是不同的治疗和不同类型。换句话说,我想比较已处理的 A 与未处理的 A,或已处理的 A 与已处理的 B,但不是已处理的 A 与未处理的 B。
如果这个问题包含在另一个问题中,我很抱歉。我搜索了论坛,但没有找到任何可以应用于这种情况的内容。谢谢。