help with glmmPQL
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