[R] obtaining means/SD after fitting a mixed model

Spencer Graves spencer.graves at pdf.com
Thu Mar 2 04:08:55 CET 2006

	  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

More information about the R-help mailing list