emmeans 模型与交互和控制变量的事后比较

问题描述

我正在使用这个混合模型，包括固定效应和随机截距和斜率：

model=lmer(volt ~ (congruency+condition+channel)^2-congruency:channel+sessorder+(1+congruency_code*(condition_code+condition_codelevel3)+sessorder_code+sessorder_codelevel3+channel_code+channel_codelevel3|pp),REML=T,data=data,control=contr1)

所有固定效应都是分类变量（congruency，2levels；condition，3levels；channel，3levels）和sessorder被用作控制变量（3levels）。

在使用 emmeans 的事后测试中，我想调查“条件*通道”之间的重要关系。为此，我使用此代码：

emms <- emmeans(model,~ condition|channel)
summary(pairs(emms),by=NULL,adjust="bonferroni")

R 中的结果表明：

Results are averaged over the levels of: congruency,sessorder
Degrees-of-freedom method: kenward-roger 
P value adjustment: bonferroni method for 9 tests 

我的问题是，这是在我的重要控制变量“sessorder”上取平均值的问题吗？我是否必须在事后测试中以某种方式考虑这个变量？

非常感谢您的意见！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）