R语法简单斜率MEM

Question

关于R上混合效应模型的语法的问题。

我已经运行了以下代码来检查简单的斜率，以确定我的一个变量(可变性)在另一个变量(模糊度)中的影响：

lmer.E1.v2 <- lmer(logRT ~ Variability.c / Ambiguity.c + (Variability.c + Ambiguity.c|ID),
                   data=data %>% filter(Experiment == "E1"),
                   control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=2e5)))
summary(lmer.E1.v2)

当我颠倒这两个变量时，代码看起来像这样：

lmer.E1.v2 <- lmer(logRT ~ Ambiguity.c / Variability.c + (Ambiguity.c + Variability.c|ID),
                   data=data %>% filter(Experiment == "E1"),
                   control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=2e5)))
summary(lmer.E1.v2)

。。我在第一段代码中得到的输出与第二段代码中的输出不同。在语法中颠倒两个变量的顺序在解释上有什么不同？

Answer 1

主要问题是/运算符是不可交换的(即a/b != b/a)：a/b扩展为a + a:b，而b/a扩展为b + a:b。您应该获得相同的总体拟合(预测、似然等)，至少达到一定程度的数值模糊，但模型参数化将有所不同。

确实存在这样的情况，(a+b|g)给出了与(b+a|g)不同的答案(参见here，但这并不常见)。