由于mblogit模型属于多层次模型,因此不能直接使用Anova函数进行分析。需要使用multcomp包中的glht函数来执行post hoc比较。以下是一个示例代码:
library(mlogit) data("Fishing", package = "mlogit") Fishing$id <- 1:nrow(Fishing) fish.mlogit <- mlogit(mode ~ price | income, data = Fishing) fish.mblogit <- mlogit(mode ~ price | income | 1, data = Fishing) summary(fish.mblogit) library(multcomp) fish.glht <- glht(fish.mblogit, linfct = mcp(mode = "Tukey")) summary(fish.glht)
这将计算出mblogit模型的概要信息以及post hoc比较的结果。