glmmTMB ZINB混合效应模型拟合不收敛问题求助

问题背景 #

• 拟拟合含混合效应的ZINB模型，模型包含三个表示嵌套分层分组结构的变量的随机截距，以及含165个时间点的var2变量的随机斜率
• 仅含随机截距的模型可正常输出结果，但加入随机斜率后模型运行报错，无法生成有效输出

模型公式 #

var9 ~ var13 + var14 + var19 + var20 + var34 + var46 + var47 + var50 + var51 + var52 + var56 + var57 + var59 + var60 + var62 + var65 + var66 + var89 + var90 + var91 + var92 + var93 + var94 + (var2 | var4/var3/var1)

已尝试操作及结果 #

• 按包开发者建议，将随机斜率中的时间变量var2指定为连续型后重新拟合：
  • 在约为原数据集1%的测试子集上可返回模型结果，但AIC、BIC、loglik、deviance均为NA值
  • 在全量数据集上运行失败，报错信息如下：

    In (function (start, objective, gradient = NULL, hessian = NULL,  :
      NA/NaN function evaluation
    Error in ff(x, order = 1) : 
      inner newton optimization failed during gradient calculation
    outer mgc:  NaN 
    Error in (function (start, objective, gradient = NULL, hessian = NULL,  : 
      NA/NaN gradient evaluation
    Calls: glmmTMB ... with.default -> eval -> eval -> do.call -> <Anonymous>
    Timing stopped at: 5258 879.3 6540
    Execution halted
    

数据集基本属性 #

• 总观测数量：2429283条
• var4的水平数：53
• var4与var3的交互项水平数：375
• var4、var3、var1三者的交互项水平数：61125

求助需求 #

恳请指导如何让模型收敛并获得有效输出，可提供测试数据集及其他所需细节，感谢帮助。

内容的提问来源于stack exchange，提问作者StatisticsFanBoy