Good morning,
I used in R contr.sum for the contrast in a lme model:
> options(contrasts=c("contr.sum","contr.poly"))
> Septo5.lme<-lme(Septo~Variete+DateSemi,Data4.Iso,random=~1|LieuDit)
> intervals(Septo5.lme)$fixed
lower est. upper
(Intercept) 17.0644033 23.106110 29.147816
Variete1 9.5819873 17.335324 25.088661
Variete2 -3.3794907 6.816101 17.011692
Variete3 -0.5636915 8.452890 17.469472
Variete4 -22.8923812 -10.914912 1.062558
Variete5 -10.7152821 -1.865884 6.983515
Variete6 0.2743390 9.492175 18.710012
Variete7 -23.7943250 -15.070737 -6.347148
Variete8 -21.7310554 -12.380475 -3.029895
Variete9 -27.9782575 -17.480555 -6.982852
DateSemi1 -5.7903419 -1.547875 2.694592
DateSemi2 3.6571596 8.428417 13.199675
attr(,"label")
[1] "Fixed effects:"
How is it possible to obtain a return with the name of my factor's levels as
with contr.treatment ?
Thanks for you help.
--
Ghislain Vieilledent
30, rue Bernard Ortet 31 500 TOULOUSE
06 24 62 65 07
[[alternative HTML version deleted]]
______________________________________________
[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
