在 lmer 输出上运行 car Anova 后列表类型错误答案

【问题标题】：List type error after running car Anova on lmer output在 lmer 输出上运行 car Anova 后列表类型错误
【发布时间】：2017-09-01 02:28:37
【问题描述】：

我正在尝试在 R 中做一个随机效应模型。我想在输出上运行 anova，我查看了一些教程并复制了他们的示例，但我的没有工作。我在 Anova 上不断收到以下错误：

(1:length(names))[-which.term][sapply(names[-which.term] 中的错误， function(term2) is.relative(term, : 无效的下标类型'list'

我正在运行的代码如下，任何人对问题所在有任何想法将不胜感激！

library(lme4)
library(car)

div1<-rep(1,5)
div2<-rep(2,3)
div3<-rep(3,3)
div6<-rep(6,6)
div<-c(div1,div2,div3,div6)
res<-c(3.8082479,7.7819745,3.3792467,7.2288647,3.4564646,1.8043898,5.1443293,3.9467614,2.5922306,1.9996585,4.2004104,0.7290807,2.1854365,3.4118980,3.2464388,2.9607496,1.9993038)
df<-data.frame(div,res)

randeff<-lmer(res~1+(1|div),data=df,REML=FALSE)
summary(randeff)
Anova(randeff)

【问题讨论】：

标签： r lme4 anova

【解决方案1】：

car::Anova 仅用于分析固定效果。你试图分析一个没有随机效应的模型，所以它变得不开心（不幸的是，错误消息是难以理解的）。要检验这个假设，请在您的数据集中添加一个随机列并为模型添加一个固定效应：

> df$y <- rnorm(nrow(df))
> randeff2 <- lmer(res~y+(1|div),data=df,REML=FALSE)
> Anova(randeff2)
Analysis of Deviance Table (Type II Wald chisquare tests)

Response: res
   Chisq Df Pr(>Chisq)
y 0.0418  1      0.838

【讨论】：