Displaying 1 result from an estimated 1 matches for "avgud89".
Did you mean:
avgud
2010 Mar 03
1
Zero inflated negative binomial
...l:
> 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 solu...