[R-sig-ME] Error in reS$ziReStruc[[i]]

[Ben Bolker]

 Sorry, hit 'send' by accident.

 ~village would treat village as a fixed effect in the zero-inflation
model; ~1|village [or ~(1|village)] would treat it as a random effect.
I can't reproduce a problem with summary(); here's a reproducible
example.  Can you construct a reproducible example that fails ?

dd <- data.frame(village=factor(rep(1:6,50)))
set.seed(101)
zp <- c(0.2,0.8,0.3,0.5,0.6,0.7)
dd$z <- rbinom(300,size=1,prob=zp[dd$village])
dd$z[dd$z==1] <- rpois(sum(dd$z),lambda=3)
library(glmmTMB)
g1 <- glmmTMB(z~1,ziformula=~1|village,data=dd,
        family=poisson)
g2 <- update(g1,ziformula=~village)
summary(g1)
summary(g2)
>
> On Mon, Aug 13, 2018 at 9:11 AM C. AMAL D. GLELE <altessedac2 using gmail.com> wrote:
> >
> > Hi, dear all.
> > I've built a glmm (model) with glmmTMB; it had fitted without,  neither
> > error message, nor warnings;
> > but, since I've added the following ziformula=~village (village is a 6
> > levels grouping factor variable), I got the warning message
> > "Model convergence problem; non-positive-definite Hessian matrix"
> > Trying to resolve this, I've used ziformula=~(1|village) instead;
> >  then, the model fits without neither error message, nor warnings; but when
> > I'd ran "summary(model)", I got the following error message "Error in
> > reS$ziReStruc[[i]] : subscript out of bounds": it's my first time to get
> > such an error message and I don't know any idea about its possible cause;
> > so, any help to solve it will be much appreciated.
> > In advance, thanks.
> > Kind regards,
> >
> >         [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-mixed-models using r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models