Hello!
If I run this example:
library(nlme)
fm1 - lme(distance ~ age+Sex, Orthodont, random = ~ age + Sex| Subject)
If I run:
summary(fm1)
then I can see the fixed effects for age and sex (17.7 for intercept,
0.66 for age, and -1.66 for SexFemale)
If I run:
ranef(fm1)
Then it looks like it's
Dmitri:
Not quite sure what you mean by easier ... fixef() and ranef() will
both give coefficients which can be easily manipulated to produce the
results for all subjects.
However, note that there are numerous built-in lme
functions(especially for graphics) that do this internally to produce,
Thank you very much, but not I am not sure now - does ranef(fm1) give
the (total) slope and
intercept values directly for each group or not?
Thanks a lot for clarifying - because I might well have been wrong.
Dimitri
On Mon, Oct 18, 2010 at 5:57 PM, Bert Gunter gunter.ber...@gene.com wrote:
Oh -- I get your question (I think). Not the total, just the random
effects. You have to add them to the fixed effects.
See e.g. p. 39 of Bates and Pinheiro.
-- Bert
On Mon, Oct 18, 2010 at 3:00 PM, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
Thank you very much, but not I am
Yes, sorry for the confusion. Maybe I should have used a different term.
So, I guess, I was right - it gives only the random effects that I
have to add to the fixed effects.
And there is no way to get it done by R (not that I can't do it myself)?
Dimitri
On Mon, Oct 18, 2010 at 6:24 PM, Bert
5 matches
Mail list logo