R-lmer相互作用效应

我有一些关于lme4软件包的R中线性混合模型分析中的3向相互作用效应的问题

我的测试结果是“MMSE”，这是一种认知测试。我主要感兴趣的是MMSE随时间（时间变量=“fuclin”）的下降（按组（“ATgroups”1,2,3,4）和胆固醇（“BL_-CHOL”0,1）。这些分析还针对几个协变量（年龄、性别、教育程度和APOE dich）进行了校正。数据集称为“ib”

这就是我的模型的样子：

总结（lmer（MMSE~apoe|dich+EDUC+性别+年龄分+numfactors+fuclin*BL|CHOL*ATgroups+（1+fuclin | ID），ib））

这就是我的输出结果：

所以在这个输出中，它总是将交互作用的最低组（ATgroups=1和BL_CHOL=0）作为参考组。但是如果我使用冒号（：）将交互添加到模型中，我会得到不同的输出。所以我的代码是：

总结（lmer（MMSE~apoe|dich+EDUC+性别+年龄分+numfactors+fuclin:BL|CHOL:ATgroups+（1+fuclin | ID），ib））

然后我的输出如下所示：

在这种情况下，我不确定从哪个比较中生成p值。我怎样才能得到与输出1中的对比度相似的对比度，但有我想要的比较（而不仅仅是与最低类别的比较）

---

谢谢

也许您可以在运行分析之前设置变量的参考级别。是一个很好的例子。也许您可以在运行分析之前设置变量的引用级别。这是一个很好的例子。