我的数据是零膨胀的，所以我正在使用 glmmamdb 运行一个零膨胀模型：
Model3z <- glmmadmb(Count3 ~ Light3 + (1|Site3), zeroInflation = T, family= "poisson", data = dframe3)
但是，当我尝试对此模型中的不同灯光类型进行成对比较时pwcs3 <- lsmeans(Model3z, "Light")我收到错误消息：
ref_grid(object, ...) 中的错误：无法处理“glmmadmb”类的对象 使用 help("models", package = "emmeans") 获取有关受支持模型的信息。
当我访问 emmeans 软件包网站时，它说 glmmadmb 不再受支持。
我已切换到 pscl 和 zeroinfl 函数，但不确定如何重组我的代码以适应 pscl 格式。键入P <- zeroinfl(Count3 ~ Light3 + (1|Site3), family = poisson, data = dframe3)得到错误消息：
contrasts<-( , value = contr.funs *tmp*[1 + isOF[nn]]) 中的错误：对比只能应用于具有 2 个或更多级别的因子此外：警告消息：在 Ops.factor(1, Site3) 中：'|' 对因子没有意义
是否有另一种方法可以将 glmmadmb 与 lsmeans 一起使用？如果没有，有谁知道 pscl 中的零膨胀模型代码应该是什么样子？谢谢。