Hi all, I am running the following model:
> glm89.nb <- glm.nb(AvGUD ~ Year*Trt*Micro) where Year has 3 levels, Trt has 2 levels and Micro has 3 levels. However when I run it has a zero inflated negative binomial (as I have lots of zeros) I get the below error message: > Zinb <- zeroinfl(AvGUD ~ Year*Trt*Micro |1, data = AvGUD89, dist = "negbin") Error in optim(fn = loglikfun, gr = gradfun, par = c(start$count, start$zero, : non-finite value supplied by optim For what I have read I think the problem is that for Year level 3 there is no Trt 1 and for Year level 1 there is no Micro level 3. I cannot find a solution to this problem, is there any way I can solve this so I can run the zero inflated model? Thank you for your time Rebecca [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.