[R-sig-ME] Find Hessian in MixCat and glmmTMB packages
Sun, John
j@un20 @end|ng |rom @|b@ny@edu
Sat Sep 24 18:25:39 CEST 2022
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
More information about the R-sig-mixed-models
mailing list