r-lmm中的事后分析:未执行

Question

我正在尝试解决错误：

if((misc$estType=="pairs"(&amp；(paste(c("，by(，collapse="(！=：缺少需要TRUE/FALSE的值

我的代码：

model1 <- lmer (Y ~ X1 * X2 + (1|"FACTOR of 13 levels"), data=data, REML = FALSE)
library (emmeans)
post_hoc <- emmeans (model1, ~ X1*X2)
pairs (post_hoc, adjust="tukey")

但它返回上述错误。我试图更改代码，但没有成功。有人能帮我解决吗？

Answer 1

我想我知道发生了什么。我相信X1和X2都是数字预测因子。因此，您的post_hoc只列出了一种情况，即X1和X2分别处于其平均值的情况。试试这个：

post_hoc

(只需键入对象的名称并点击回车。这将显示它。(

我敢打赌，只有一个案例列出，但在这个网站上下注可能是不受欢迎的。你无法比较一件事，这就是为什么你会出错的原因。

如果这些预测因子真的应该是因子，那么你需要拟合一个识别它们的模型，例如

model2 <- lmer(factor(X1) * factor(X2) + (1|subject), 
data = data, REML = TRUE)

你的model1几乎肯定是无用的，因为它适用于两个数值变量中的双曲面方程。

我会看看是否可以让contrast()和pairs()显示一条信息更丰富的错误消息。