    I have fit a model with glmmPQL function in MASS library.

I fit a binomial longitudinal response variable nested in 17 stations.

I would like to know how I can obtain elements of diagnostic checks about these models in order to choose best model.

I use summary(), but can I use other functions like in lme, for example anova?

I would be thankfull for all the insights.

Fabrizio Consentino.

	[[alternative HTML version deleted]]


