R中的边际方差分析

Question

我试图使用以下代码运行一个边际(也称为type-III)的ANOVA。不幸的是，我得到了“错误：$操作符对于原子向量无效”。谷歌搜索发现了另一个人收到了同样的错误(请看这里)，但不幸的是没有解决方案。

rm(list = ls())

data(iris)

iris.lm <- lm(formula = Sepal.Length ~ Sepal.Width + Petal.Length , data = iris)
print(anova(object = iris.lm))
print(anova(object = iris.lm , type = "marginal"))

Answer 1

作为以后的参考，来自包Anova的函数car接受type参数(如"III"或3)。此函数不应与标准anova库中的函数R混淆，后者不接受type参数。

Answer 2

使用traceback，您可以本地化错误：

print(anova(object = iris.lm , type = "marginal"))
Error: $ operator is invalid for atomic vectors
> traceback()
7: deparse(x$terms[[2L]])
6: FUN(X[[2L]], ...)
5: lapply(objects, function(x) deparse(x$terms[[2L]]))
4: anova.lmlist(object, ...)
3: anova.lm(object = iris.lm, type = "marginal")
2: anova(object = iris.lm, type = "marginal")
1: print(anova(object = iris.lm, type = "marginal"))

因此，您会得到错误，因为您试图在不支持它的对象上调用terms。这将复制错误：

 lapply(list(iris.lm,type='marginal'),terms)
 Error: $ operator is invalid for atomic vectors

或者只是：

terms('marginal')

但你为什么认为这会奏效呢？也许我错过了什么，但我认为没有提到阿诺瓦可以把type作为一个论点。