Skip to content
Prev 18250 / 20628 Next
Original Raw

Error from glmmTMB().

Ben Bolker 
I did a little more poking around with this.

  It works (i.e. fits with no warnings with binomial() instead of
betabinomial().

  It works, sort of (i.e. gives 1 warning about NaN evaluation and a
non-pos-def Hessian warning, probably from an effectively singular
model, but gives sensible fixed-effect estimates and std devs) with
link=logit instead of cloglog.

   Debugging my way through glmmTMB:::fitTMB(), even the first
evaluation of the objective function gives NaN, and the gradient is a
vector of NaN values.  I suspect we have a failure to clamp the cloglog
function somewhere, and may know how to fix it.

  This is a very nice, simple example that may help us improve a bunch
of stuff (it will also help us put some tests in useful places to get
*slightly* less opaque error messages).

  I know this doesn't solve your problem immediately - but if I'm lucky
and have a bit of time I may be able to implement something soon.

  thanks Rolf!

  Ben
On 2020-03-17 9:50 a.m., Mollie Brooks wrote:

Thread (12 messages)

Rolf Turner Error from glmmTMB(). Mar 14 Mollie Brooks Error from glmmTMB(). Mar 17 Dimitris Rizopoulos Error from glmmTMB(). Mar 18 Rolf Turner Error from glmmTMB(). Mar 18 Ben Bolker Error from glmmTMB(). Mar 18 Ben Bolker Error from glmmTMB(). Mar 18 Rolf Turner Error from glmmTMB(). Mar 19 Rolf Turner Error from glmmTMB(). Mar 19 Rolf Turner Error from glmmTMB() --- no joy. Mar 19 Ben Bolker Error from glmmTMB() --- no joy. Mar 19 Rolf Turner Error from glmmTMB() --- no joy. Mar 19 Ben Bolker Error from glmmTMB() --- no joy. Mar 19