Anova函数无法在mblogit模型上使用,原因是mblogit模型不是包含单一响应变量的线性模型。因此,可以尝试使用多元响应变量线性模型,然后使用Anova函数来进行方差分析。
以下是一个使用multinom函数拟合多元响应变量线性模型,并使用Anova函数进行方差分析的示例:
library(nnet)
library(car)
data("Fishing", package="car")
# 建立多元响应变量线性模型
model <- multinom(cbind(never,less1,1.more1) ~ ., data=Fishing)
# 进行方差分析
Anova(model, type="III")
这里使用了Fishing数据集中的3个响应变量:never、less1和1.more1,分别表示钓鱼次数从未、1个或更少和1个或更多的人数。此外,Car包中的Anova函数支持type选项,可以选择使用type III 方差分析,这与SPSS中的分析相同。