Не сходящийся glmmTMB - PullRequest
Новая
       17

Не сходящийся glmmTMB

0 голосов
/ 17 мая 2019

Я работаю с многомерной моделью с гамма-распределением, и я хотел бы использовать синтаксис lme4, развернутый в glmmTMB, однако я заметил кое-что странное в моей модели.Видимо, модель легко сходится, когда я использую stats::glm, но, похоже, выдает ошибку в glmmTMB framework.Вот воспроизводимый пример: 

d <- data.frame(gamlss.data::plasma) # Sample dataset

m4.1 <- glm(calories ~ fat*fiber, family = Gamma(link = "log"), data = d)    # Dos parámetros con interacción
m4.2 <- glmmTMB(calories ~ fat*fiber, family = Gamma(link = "log"), data = d)    # Dos parámetros con interacción

>Warning message:
In fitTMB(TMBStruc) :
  Model convergence problem; non-positive-definite Hessian matrix. See vignette('troubleshooting')

Полагаю, решение может лежать на параметрах управления, но, посмотрев виньетку устранения неполадок, я не уверен, с чего начать.

...