用 lmerTest 创建的模型的方差同质性

Question

我已经运行了模型：

model2 <- lmer(tas ~ station
              + (1|date), data = all5)

其中 station 是具有 4 个水平的分类变量。 我想检查方差同质性。 我从

lmerresid <- resid(model2)

中找到残差。 如果我运行 levene 测试：

leveneTest(lmerresid ~ all5$station)

我得到 p = 0.6572 所以它意味着方差齐性。 如果我制作箱线图，则意味着方差同质性：

boxplot(lmerresid ~ all5$station)

但是，如果我绘制残差与拟合图，那么就会出现错误： #情节与因素

par(mfrow = c(1,1))
plot(lmerresid ~ fitted(model2))

方差是否齐次？

Answer 1

我不确定你是否在混合模型上使用了 levene 测试（而不是

lm()

），尽管我在那里可能是错的。

说的

我发现测试模型有效性的最简单方法是使用

library(performance)

.

performance(model2)

“超级简单，勉强和不便”。我不会说它对您的想法来说是完美的*，但它会运行您可能需要检查模型的大部分（如果不是全部）测试。这与标准

lmer()

输出一起，应该使您能够解释您的模型。

*

您需要深入阅读手册以了解测试的所有适当用途。