For better or worse, it's holidays in the states. Very amusing for me
being in a non-Thanksgiving celebrating country.
In addition, it's not a problem. The complaint is valid. Probably no
one has coded up the right solution yet for comparison.
I can't recall if one would want those statistics for a binomial
random effects model, but I do recall some issues with model
comparison in that setting, though they are a bit dated (say, 2 years
or so).
On Fri, 26 Nov 2004 09:31:40 +0700, Andrew Criswell
<r-stats at arcriswell.com> wrote:> Hello:
>
> Will someone PLEASE help me with this problem. This is the third time
> I've posted it.
>
> When I appply anova() to two equations estimated using glmmPQL, I get a
> complaint,
>
> > anova(fm1, fm2)
>
> Error in anova.lme(fm1, fm2) : Objects must inherit from classes
"gls",
> "gnls" "lm","lmList",
"lme","nlme","nlsList", or "nls"
>
> >
>
> The two equations I estimated are these:
>
> > fm1 <- glmmPQL(choice ~ day + stereotypy,
>
> + random = ~ 1 | bear, data = learning, family = binomial)
>
> > fm2 <- glmmPQL(choice ~ day + envir + stereotypy,
>
> + random = ~ 1 | bear, data = learning, family = binomial)
>
> Individually, I get results from anova():
>
> > anova(fm1)
>
> numDF denDF F-value p-value
> (Intercept) 1 2032 7.95709 0.0048
> day 1 2032 213.98391 <.0001
> stereotypy 1 2032 0.42810 0.5130
>
> >
> > anova(fm2)
>
> numDF denDF F-value p-value
> (Intercept) 1 2031 5.70343 0.0170
> day 1 2031 213.21673 <.0001
> envir 1 2031 12.50388 0.0004
> stereotypy 1 2031 0.27256 0.6017
>
> >
>
> I did look through the archives but didn't finding anything relevant to
> my problem.
>
> Hope someone can help.
>
> ANDREW
> ____________________________
> _
> platform i586-mandrake-linux-gnu
> arch i586
> os linux-gnu
> system i586, linux-gnu
> status
> major 2
> minor 0.0
> year 2004
> month 10
> day 04
> language R
>
> --
> Andrew R. Criswell, Ph.D.
> Graduate School, Bangkok University
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
>
--
best,
-tony
---
A.J. Rossini
blindglobe at gmail.com