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

