我正在运行零膨胀的glmmTMB模型,其计数数据可用here。当我的完整模型运行时没有任何错误或警告时,子模型“ zm13”会生成错误消息和2条其他警告消息。该错误似乎表明参数设置过高,但由于完整模型尚可,因此我不确定在解决问题方面下一步该怎么做。
我的模型(完整的和产生子集的错误消息)如下:
fullEVCounts1 <- glmmTMB(ev ~ z.jdate + z.start + z.tourist + fp + offset(log(duration)) + (1|name), zi = ~ 1, family = poisson(), data = EVCountsSTD, na.action = "na.fail")
zm13<- glmmTMB(ev ~ z.start + fp + offset(log(duration)) + (1|name), zi = ~ 1, family = poisson(), data = EVCountsSTD, na.action = "na.fail")
完整模型和变量的所有可能组合(构成子集模型)运行时没有错误或警告消息,而“ zm13”产生以下消息:
Error in optimHess(par.fixed, obj$fn, obj$gr) :
gradient in optim evaluated to length 1 not 5
In addition: Warning messages:
1: In (function (start, objective, gradient = NULL, hessian = NULL, :
NA/NaN function evaluation
2: In (function (start, objective, gradient = NULL, hessian = NULL, :
NA/NaN function evaluation
对于此模型子集为什么会弹出此错误消息和相关的警告消息,有任何想法吗?预先感谢您提供的任何指导!