Question

我是R语言编程的新手，我需要为数据使用零膨胀混合效果广义线性模型。我在网上发现，使用glmmTMB软件包是最好的选择。我有代码并且可以运行，但是我缺少一些输出-它显示为“ NA”。我可以解决什么问题，使它给我值而不是“ NA”？我还不断收到警告，说“模型收敛问题；非正定Hessian矩阵”，不明白那意味着什么，有帮助吗？

我已经在故障排除页面（https://cran.r-project.org/web/packages/glmmTMB/vignettes/troubleshooting.html）上在线查询了错误代码，但解释却很麻烦。

这是我的代码的样子（有注释掉的变量，所以我们现在只能查看特定的变量）：

No.WavesModel2 <- glmmTMB(No.WavesScaled~(MudTemp+I(MudTemp^2)+
                                      #Ba+I(Ba^2)+
                                      #Cd+I(Cd^2)+
                                      Pb+I(Pb^2)+
                                      #Sr+I(Sr^2)+
                                      #Zn+I(Zn^2)+
                                      #Fe+I(Fe^2)+
                                      Mn+I(Mn^2)+
                                      #Cu+I(Cu^2)+
                                      CarIndex+I(CarIndex^2)+
                                      No.Pedestrians)+I(No.Pedestrians^2)
                            (1|Site),
                          data=MaleWavingR,
                          ziformula=~1,
                          family=poisson)
summary(No.WavesModel2)

这些是显示在我的输出中的警告消息：

1: In nlminb(start = par, objective = fn, gradient = gr, control = control$optCtrl) :
  NA/NaN function evaluation
2: In nlminb(start = par, objective = fn, gradient = gr, control = control$optCtrl) :
  NA/NaN function evaluation
3: In fitTMB(TMBStruc) :
  Model convergence problem; non-positive-definite Hessian matrix. See vignette('troubleshooting')

我在R中使用glmmTMB包在输出中不断得到“ NA”

0 个答案: