Displaying 1 result from an estimated 1 matches for "bushid".
Did you mean:
bus_id
2009 Oct 26
1
GLMMPQL and negbinomial: trouble with the X-axis in PREDICT
...therefore not to
scale.
I did try creating a new vector to plot the data, but I got an error
message (see below).
Does anyone have suggestions on where to look to solve this problem?
Erin Latham,
M.GIS Candidate, Geography
University of Calgary, AB, Canada
> m0 <- glm.nb(MANUAL ~ BERRIES + BUSHID , data=a)
>theta<-m0$theta
> m1 <- glmmPQL(MANUAL ~ BERRIES, random = ~ 1|BUSHID, data=a, family=negative.binomial(theta))
> summary(m1)
Linear mixed-effects model fit by maximum likelihood
Data: a
AIC BIC logLik
NA NA NA
Random effects:
Formula: ~1 | BUSHID
(Int...