[ADMB Users] Warning message: Convergence failed:log-likelihood of gradient=...
W Robert Long
W.R.Long at leeds.ac.uk
Fri Jan 11 13:02:30 PST 2013
Hi all
I am using glmmADMB to fit a zero-inflated negative binomial model with
random slopes and random intercepts. The model formula is
zi1 <- glmmadmb(Y ~ time*X1 + time*X2 +
(time | Subject), data=final, family="nbinom2", zeroInflation=TRUE)
There are 4 time points and 85 subjects.
The outcome has some missing values, approx 20% are missing. The model
runs fine with complete case data.
I had originally wanted to do multiple imputation, but I have yet to
find a way of multiply imputing zero-inflated data in the context of
this kind of multilevel/hierarchical model. So, I have implemented a
random hot-deck imputation, which /appears/ to have worked well,
however, around 50% of the completed datasets return this warning from
glmmADMB:
Convergence failed:log-likelihood of gradient= -XX.XXX
and the parameter estimates are not as expected. Also, occasionally
(5-10% of the time) it returns the error:
The function maximizer failed (couldn't find STD file)
The remaining imputed datasets don't cause any problem and the parameter
estimates are within the ranges expected.
All the imputed values are plausible, as you would expect from a
hot-deck imputation and so far I can't identify what is causing these
warnings and errors.
I would be grateful for any hints or advice about how to proceed.
Thanks
Robert Long
More information about the Users
mailing list