在lmer输出

Question

我试图在R中做一个随机效果模型。我想在输出上运行一个anova，我已经浏览了几个教程并复制了他们的例子，但是我没有＆＃ 39;工作。我在Anova上一直收到以下错误：

  

（1：length（names））中的错误[ - which.term] [sapply（names [-which.term]，   function（term2）is.relative（term，：invalid subscript type＆＃39; list＆＃39;

我正在运行的代码如下，任何人对问题的看法都会非常感激！

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)

Answer 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