Raw Message

Message-ID: <BL0PR04MB4563778B131E342C305020ACB9509@BL0PR04MB4563.namprd04.prod.outlook.com>
Date: 2022-09-24T16:25:39Z
From: Sun, John
Subject: Find Hessian in MixCat and glmmTMB packages

Dear All,

I am writing to ask how to find the hessian from these fitted objects to calculate the condition number of the hessian to diagnose loss of precision.
Please let me know if you know any information.

library(mixcat)
data(schizo)
attach(schizo)
a<-npmlt(y~trt*sqrt(wk),formula.npo=~trt,random=~1+trt,id=id,k=2,EB=FALSE)

library(glmmTMB)
(m1 <- glmmTMB(count ~ mined + (1|site),
  zi=~mined,
  family=poisson, data=Salamanders))

I cannot find how to find the hessian from these fitted objects.

Best regards,
John