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