[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