R help - Feb 2006 - obtaining means/SD after fitting a mixed model

Hello,
I am running (non-balanced) mixed models (using library lme4) such as :
model1<-glmmPQL(Y~grouping variable+ covariate, random=~1|yr/month, 
data=dat, family= gaussian or poisson)
and 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

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]]

Have you consulted Pinheiro and Bates (2000) Mixed-Effects Models in 
S and S-Plus (Springer)?  The MASS function "glmmPQL" translates the 
problem into approximate "gaussian" terms and passes the translation 
into "lme".  Therefore, I would expect much of what appears in
Pinheiro
and Bates to be accessible to you.

	  If you would like more help from this listserve, please submit 
another question.  Please be advised, however, that there is a fairly 
high correlation between the speed and utility of the responses on the 
one hand and the clarity of the question on the other.  If you would 
like to increase your chances of getting a useful answer quicker, I 
encourage you to PLEASE read the posting guide! 
"www.R-project.org/posting-guide.html", especially the part about 
supplying a simple, self-contained example:  Vague questions get vague 
answers (if any).  Simple, clear questions are more likely to get 
simple, clear and immediate answers.

	  spencer gravesd

Pablo Inchausti wrote:
> Hello,
> I am running (non-balanced) mixed models (using library lme4) such as :
> model1<-glmmPQL(Y~grouping variable+ covariate, random=~1|yr/month, 
> data=dat, family= gaussian or poisson)
> and 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
> 
> ______________________________________________
> 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