obtaining means/SD after fitting a mixed model

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
2 messages Options
Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

obtaining means/SD after fitting a mixed model

Pablo Inchausti-2
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

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: obtaining means/SD after fitting a mixed model

Spencer Graves
          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
>
> ______________________________________________
> [hidden email] mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Loading...