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]]
______________________________________________
[email protected] 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.