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