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

[Holly G]

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,
Holly

glmmadmbZINBtotalCRATdensInt1fixGrpSrvy <- glmm.admb(
fixed=TotalCRATdensInt~
Dist,
group="fSurveyID",
data= Data110217noFlNaRVfallsub, 
zeroInflation=FALSE,
family="nbinom")

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

library(MASS)

glmnb1totalCRATdensInt1fix <- glm.nb(
TotalCRATdensInt~
Dist,
data= Data110217noFlNaRVfallsub,
init.theta=1,
link = log,
control=glm.control(maxit = 25, trace = TRUE))

library(pscl)

psclZINB1totalCRATdensInt1fix <- zeroinfl(
formula= (TotalCRATdensInt~Dist),
dist="geometric",
link="logit",
data=Data110217noFlNaRVfallsub)

> 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>