[R-sig-ME] (no subject)
Renwick, A. R.
a.renwick at abdn.ac.uk
Tue Jan 27 10:51:54 CET 2009
If I understand your question correctly you need to relevel your variable:
age<-relevel(age,ref="3")
Then rerun the model and age3 should be used as the reference group
-----Original Message-----
From: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf Of Lawrence Lee
Sent: 27 January 2009 09:24
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] (no subject)
Dear all,
Does anyone know how to change the reference group within lme().
For example,
lme(fixed = obs ~ line + age - 1,random = ~1 |sire, data = lambweight)
Suppose age has 3 levels, so the code above will give me estimate of age2 and age3:
Fixed effects: obs ~ line + age - 1 Value Std.Error DF t-value p-valueline1 10.491153 0.7261881 18 14.446882 0.0000line2 12.290287 0.8242309 18 14.911219 0.0000line3 11.032864 0.7647870 18 14.426061 0.0000line4 10.276735 0.7389137 18 13.907897 0.0000line5 10.952840 0.6084109 18 18.002373 0.0000age2 -0.155435 0.7157030 38 -0.217178 0.8292age3 0.009646 0.5481034 38 0.017599 0.9861
How I can change within lme() in order to give me the estimate of age1 and age2.
Thanks,
Lawrence
_________________________________________________________________
ore_012009
[[alternative HTML version deleted]]
_______________________________________________
R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
The University of Aberdeen is a charity registered in Scotland, No SC013683.
More information about the R-sig-mixed-models
mailing list