lmm 中的事后分析:未执行

Posthoc analysis in lmm: unexecuted

我正在尝试解决错误:

if ((misc$estType == "pairs") && (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")

但是returns出现了上面的错误。我试图更改代码,但我没有成功。 请问有人可以帮我解决吗?

我想我知道发生了什么。我相信 X1X2 都是数字预测变量。因此,您的 post_hoc 只列出了一种情况 - X1X2 均处于平均值的情况。试试这个:

post_hoc

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

我敢打赌只列出了一个案例,但在这个网站上投注可能会让人皱眉。你不能比较一件事,这就是你得到错误的原因。

如果这些预测变量真的应该是因素,您需要拟合一个识别它们的模型,例如,

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

你的model1几乎可以肯定是没用的,因为它符合两个数值变量的双曲曲面方程。

我看看能否让 contrast()pairs() 显示更详细的错误消息。