[ADMB Users] Error in a negative binomial, zero-inflated glmmADMB with random effect

Holly G hollygoy at hotmail.com
Wed Mar 30 17:35:46 PDT 2011

Hi, Hans,
thanks for the pointer, unfortunately, when I use a variable that has already been converted to an integer, I get the same error in all cases: "Incompatible array bounds". This occurs whether zeroInflation is TRUE or FALSE. Using the simplest model, with only one fixed variable, both glm.nb{MASS} and zeroinfl{pscl} converge, whereas glmm.ADMB does not. Therefore, I am inclined to ask if, in glmmADMB:
1. is there a way to designate an initial value of theta (to produce a geometric rather than a negative binomial distribution)?
2. are there constraints in picking a grouping factor?
Thank you for the advice,

glmmadmbZINBtotalCRATdensInt1fixGrpSrvy <- glmm.admb(
data= Data110217noFlNaRVfallsub, 

# Incompatible array bounds in dvar_vector  operator * (_CONST dvar_matrix& m,_CONST dvar_vector& x)


glmnb1totalCRATdensInt1fix <- glm.nb(
data= Data110217noFlNaRVfallsub,
link = log,
control=glm.control(maxit = 25, trace = TRUE))


psclZINB1totalCRATdensInt1fix <- zeroinfl(
formula= (TotalCRATdensInt~Dist),

> Date: Wed, 30 Mar 2011 22:32:58 +0200
> Subject: Re: [ADMB Users] Error in a negative binomial, zero-inflated glmmADMB with random effect
> From: hskaug at gmail.com
> To: hollygoy at hotmail.com
> CC: users at admb-project.org
> Hi Holly,
> It seems that "round" cannot be used in the formula. Try to add
> a new variable to the dataframe, and use this as the response.
> Hans
> On Tue, Mar 29, 2011 at 5:44 PM, Holly G <hollygoy at hotmail.com> wrote:
> > Hello list,
> >
> > I have been working with a dataset that is large, overdispersed,
> > non-normally distributed, non-linear, spatially autocorrelated, and
> > zero-inflated. So far, I have been able to account for everything except for
> > the zero-inflation using GAMMs in R. To my knowledge, there is no general
> > software code available yet to incorporate zero- inflation into negative
> > binomial GAMMs, please correct me if I am wrong. As far as GLMMs, I have
> > modeled the data with glmmPQL (accounting for random effects and
> > autocorrelation structure), and zeroinfl{pscl} (accounting for
> > zero-inflation). For a zero-inflated model that includes random effects, I
> > am attempting glmm.admb (the example from the documentation ran smoothly),
> > but have encountered a few errors, all a variation of “The function
> > maximizer failed” (below).  I am not sure if these errors are due to the
> > nature of the data or to a mistake that I overlooked, so I provide details
> > below. I would appreciate any advice. Apparently, the temporary nbmm.dat and
> > nbmm.pin files are not being created correctly, but I am unsure as to why or
> > how to fix it.
> >
> > Thank you,
> > Holly
> >
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20110331/c3872cbc/attachment.html>

More information about the Users mailing list