我正在尝试使用emmeans通过MuMin包中的model.avg从平均clmm模型生成和可视化概率。但是,我始终收到一条错误消息,指出“ vcov0 [i,i] <-v中的错误: 下标的作业中不允许使用NA。”
我尝试了以下方法提供的解决方法:Plot model fit for discrete variable, from average model 并在以下位置查看了第I组-多个模型(通过推算或平均)下提供的信息: https://cran.r-project.org/web/packages/emmeans/vignettes/models.html
但是这些都不能让我解决问题。这些建议该对象必须具有“ modelList”属性,该属性通过挖掘我创建的全局模型和子设置结果。我还确保了每个模型都将数据作为命名参数,并且这些都不能解决问题。
model1<-clmm(fairness~scenario+number_meals+tlu+gender+land_woned+bns_score+(1|original_order_id)+(1|interviewer_no),data=fairness_workbook,Hess = TRUE,na.action = "na.fail")
d <- dredge(model1,evaluate=TRUE,rank = "AICc",trace=2)
selected = get.models(d, delta<2)
AV = model.avg(selected)
attr(AV, "modelList")
em=emmeans(AV, ~ fairness | scenario, mode = "prob",data=fairness_workbook)
em
vcov0 [i,i]中的错误<-v:下标分配中不允许使用NA