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出现了上面的错误。我试图更改代码,但我没有成功。
请问有人可以帮我解决吗?
我想我知道发生了什么。我相信 X1
和 X2
都是数字预测变量。因此,您的 post_hoc
只列出了一种情况 - X1
和 X2
均处于平均值的情况。试试这个:
post_hoc
(只需键入对象的名称并点击 return。这将显示它。)
我敢打赌只列出了一个案例,但在这个网站上投注可能会让人皱眉。你不能比较一件事,这就是你得到错误的原因。
如果这些预测变量真的应该是因素,您需要拟合一个识别它们的模型,例如,
model2 <- lmer(factor(X1) * factor(X2) + (1|subject),
data = data, REML = TRUE)
你的model1
几乎可以肯定是没用的,因为它符合两个数值变量的双曲曲面方程。
我看看能否让 contrast()
和 pairs()
显示更详细的错误消息。
我正在尝试解决错误:
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出现了上面的错误。我试图更改代码,但我没有成功。 请问有人可以帮我解决吗?
我想我知道发生了什么。我相信 X1
和 X2
都是数字预测变量。因此,您的 post_hoc
只列出了一种情况 - X1
和 X2
均处于平均值的情况。试试这个:
post_hoc
(只需键入对象的名称并点击 return。这将显示它。)
我敢打赌只列出了一个案例,但在这个网站上投注可能会让人皱眉。你不能比较一件事,这就是你得到错误的原因。
如果这些预测变量真的应该是因素,您需要拟合一个识别它们的模型,例如,
model2 <- lmer(factor(X1) * factor(X2) + (1|subject),
data = data, REML = TRUE)
你的model1
几乎可以肯定是没用的,因为它符合两个数值变量的双曲曲面方程。
我看看能否让 contrast()
和 pairs()
显示更详细的错误消息。