Hello,
I am dealing (non-balanced) mixed models (using libraries lme4 and MASS) 
such as :
model1<-glmmPQL(Y~grouping variable+ covariate, random=~1|yr/month, 
data=dat, family= gaussian or poisson)
Besides the usual output in terms of the statistical significance of the 
fitted coefficients, I'd like to obtain the adjusted means and standard 
deviations as one obtains after fitting a lm model.
My problem is that I cannot find a command equivalent to "model.table 
(model, grouping variable, mean/sd)" that would apply to mixed models.
Any suggestion would be appreciated.
Thanks in advance

Pablo 
	[[alternative HTML version deleted]]


