我对R中使用lme4包的线性混合模型分析中的3向交互效应有一些疑问。
我的结果指标是" MMSE"这是一种认知测试。我主要感兴趣的是随着时间的推移MMSE的下降(时间变量=" fuclin")按组(" ATgroups" 1,2,3,4)和胆固醇(& #34; BL_CHOL" 0,1)。这些分析也针对几个协变量(age_cent,GENDER,EDUC和APOE_dich)进行了更正。数据集被称为" ib"。
这就是我的模型:
摘要(lmer(MMSE~apoe_dich + EDUC + GENDER + age_cent + numfactors + fuclin * BL_CHOL * ATgroups +(1 + fuclin | ID),ib))
这就是我的输出:
因此,在此输出中,它始终将最小的交互组(ATgroups = 1和BL_CHOL = 0)作为参考组。但是如果我使用冒号(:)向模型添加交互到模型,我会得到不同的输出。所以我的代码是:
摘要(lmer(MMSE~apoe_dich + EDUC + GENDER + age_cent + numfactors + fuclin:BL_CHOL:ATgroups +(1 + fuclin | ID),ib))
然后我的输出如下:
在这种情况下,我不确定生成p值的比较。我怎样才能得到与输出1中的对比相似的对比,但是我想要的比较(而不仅仅是与最低类别的比较)?
谢谢!