[FORGED] Re: Error from glmmTMB() --- no joy.
On 19/03/20 3:52 pm, Ben Bolker wrote:
<SNIP>
I made a small change that appears to allow the
original model to run correctly. If you want to try it out, and have
development tools (compilers etc.) and the remotes package installed,
remotes::install_github("glmmTMB/glmmTMB/glmmTMB at cloglog_fix")
should install the patched version (I think).
This will presumably be merged with the master branch sometime soon but
I want to add some tests etc.
I've finally got around to trying out your patch. No problems with the install_github() process, but after doing the install I executed:
library(glmmTMB)
fmla <- cbind(Dead, Alive) ~ (Trt + 0)/Dose + (Dose | Rep)
X <- dget("X.txt")
fit <- glmmTMB(fmla, data = X, family = betabinomial(link = "cloglog"),
dispformula = ~1
and got:
Error in (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN gradient evaluation In addition: Warning messages: 1: In (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN function evaluation 2: In (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN function evaluation 3: In (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN function evaluation 4: In (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN function evaluation 5: In (function (start, objective, gradient = NULL, hessian = NULL, : NA/NaN function evaluation Timing stopped at: 0.572 0 0.572
So my Sadim touch persists in its effect. cheers, Rolf P. S. With the logit link I got warnings
1: In (function (start, objective, gradient = NULL, hessian = NULL, :
NA/NaN function evaluation
2: In fitTMB(TMBStruc) :
Model convergence problem; non-positive-definite Hessian matrix. See vignette('troubleshooting')
but no errors as such. R.
Honorary Research Fellow Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276