| Multilocation {SASmixed} | R Documentation |
The Multilocation data frame has 108 rows and 7 columns.
This data frame contains the following columns:
B < D < E < I < G <
A < C < F < H
1 to 3
1 to 4
B/1 < B/2 < B/3 < D/1 <
D/2 < D/3 < E/1 < E/2 <
E/3 < I/1 < I/2 < I/3 <
G/1 < G/2 < G/3 < A/1 <
A/2 < A/3 < C/1 < C/2 <
C/3 < F/1 < F/2 < F/3 <
H/1 < H/2 < H/3
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 2.8.1).
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly")) str(Multilocation) ### Create a Block Multilocation$Grp <- with(Multilocation, Block:Location) (fm1Mult <- lmer(Adj ~ Location * Trt + (1|Grp), Multilocation)) anova(fm1Mult) (fm2Mult <- lmer(Adj ~ Location + Trt + (1|Grp), Multilocation)) fm3Mult <- lmer(Adj ~ Location + (1|Grp), Multilocation) fm4Mult <- lmer(Adj ~ Trt + (1|Grp), Multilocation) fm5Mult <- lmer(Adj ~ 1 + (1|Grp), Multilocation) anova(fm2Mult) anova(fm1Mult, fm2Mult, fm3Mult, fm4Mult, fm5Mult) Multilocation$LocTrt <- with(Multilocation, Location:Trt) ### Treating the location as a random effect ## Not run: (fm1MultR <- lmer(Adj ~ Trt + (1|LocTrt) + (1|Location) + (1|Grp), Multilocation)) #intervals(fm1MultR) anova(fm1MultR) fm2MultR <- lmer(Adj ~ Trt + (Trt - 1|Location) + (1|Block), Multilocation) anova(fm1MultR, fm2MultR) ## End(Not run)