[R-sig-ME] Error in reS$ziReStruc[[i]]
bbolker @ending from gm@il@com
Mon Aug 13 20:31:15 CEST 2018
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)))
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)
g1 <- glmmTMB(z~1,ziformula=~1|village,data=dd,
g2 <- update(g1,ziformula=~village)
> 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
More information about the R-sig-mixed-models