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 ? 

Thank you very much for your help,

Marie
	[[alternative HTML version deleted]]

