marieline gentes <mlgentes2@...> writes: > Dear list, > > I'm just starting to learn how to use glmmPQL and I > have a very basic question. I'm a bit embarrassed asking > this, but I could not find a clear answer anywhere, > and I did search quite a bit. Unfortunately there is no > one around me to ask - none of my collegues work with glmmPQL. > > Here is the fixed effects section of my output: > > Fixed effects: cbind(ColoYes.allnoF24, ColoNo.allnoF24) ~ year.coded + sex > > Value Std.Error DF t-value p-value > (Intercept) 1.4442270 0.1824173 74 7.917161 0.0000 > year.codedY2011 -0.1733713 0.1864744 74 -0.929732 0.3555 > year.codedY2012 -0.3004284 0.2027411 74 -1.481833 0.1426 > sexM 0.4403108 0.1507471 74 2.920857 0.0046 > I simply want to report the fact that year (categorical, three > levels) is not significant. In traditional linear models I would > simply report the p value for treatment "year" from an anova > table. But my understanding is that producing an anova > from glmmPQL is not possible (?). Reporting those individual p > values of my three years separately seems clumsy and out of place... > What is the proper way to do this ?
Since glmmPQL doesn't technically compute likelihoods, and its authors are purists, it doesn't provide standard likelihood-based model comparison tools. You could switch to glmer, from the lme4 package, or I suppose you could use Fisher's method or some equivalent to combine the p values: see this thread (which also suggests alternatives) http://r.789695.n4.nabble.com/combining-P-values-using-Fisher-s-method- td832459.html [URL broken] or http://finzi.psych.upenn.edu/R/library/MADAM/html/fisher.method.html good luck Ben Bolker _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology