Lucrez cu un model multivariat cu o distribuție Gamma și aș dori să folosesc sintaxa lme4 implementată în glmmTMB, cu toate acestea, am observat ceva ciudat cu modelul meu. Aparent, modelul converge ușor atunci când folosesc stats: glm, dar pare să dea o eroare în cadrul glmmTMB. Iată un exemplu reproductibil:

overflow

Cred că soluția ar putea sta pe parametrii de control, dar după ce am analizat vigneta de depanare, nu sunt sigur de unde să încep.

1 Răspuns 1

O soluție poate fi scalarea variabilelor (atâta timp cât sunt numerice).

Aici, cel de-al doilea model converge bine, în timp ce nu mai făcuse înainte. Cu toate acestea, rețineți diferența în estimările parametrilor între cele două modele:

Acest lucru se datorează faptului că estimările se bazează pe parametri la scară, deci trebuie să fie interpretați cu precauție sau „nescalate”. A se vedea: Înțelegerea „scalei” în R pentru o înțelegere a funcției scale () și a se vedea: interpretarea coeficienților de regresie scalată. pentru o înțelegere mai aprofundată a ceea ce înseamnă acest lucru într-un model

Ca o ultimă notă, faptul că modelele converg nu înseamnă că sunt potrivite.