Displaying 2 results from an estimated 2 matches for "quine2".
Did you mean:
quine
2004 May 21
1
Bug in update()? (PR#6902)
...the problem also
exists under R 1.9.0.
A reproducible example uses the following contingency table:
> library(MASS)
> data(quine)
> tmp <- with(quine, expand.grid(Eth=levels(Eth), Sex=levels(Sex),
+ Lrn=levels(Lrn), Age=levels(Age)))
> n <- nrow(quine)
> quine2 <- with(quine,
+ data.frame(tmp,
+ Count=as.vector(tapply(rep(1,n), list(Eth, Sex, Lrn, Age), sum))))
First fit a saturated model and see which term we can drop:
> fm <- glm(Count ~ .^4, quine2, family=poisson)
> drop1(fm, test="Chisq")
Single term del...
2010 Mar 03
1
Zero inflated negative binomial
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 =