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).
str(Multilocation) if (require("lme4", quietly = TRUE, character = TRUE)) { options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly")) ### Create a Block Multilocation$Grp <- with(Multilocation, Block:Location) print(fm1Mult <- lmer(Adj ~ Location * Trt + (1|Grp), Multilocation)) print(anova(fm1Mult)) print(fm2Mult <- lmer(Adj ~ Location + Trt + (1|Grp), Multilocation)) print(fm3Mult <- lmer(Adj ~ Location + (1|Grp), Multilocation)) print(fm4Mult <- lmer(Adj ~ Trt + (1|Grp), Multilocation)) print(fm5Mult <- lmer(Adj ~ 1 + (1|Grp), Multilocation)) print(anova(fm2Mult)) print(anova(fm1Mult, fm2Mult, fm3Mult, fm4Mult, fm5Mult)) ### Treating the location as a random effect print(fm1MultR <- lmer(Adj ~ Trt + (1|Location/Trt) + (1|Grp), Multilocation)) print(anova(fm1MultR)) print(fm2MultR <- lmer(Adj ~ Trt + (Trt - 1|Location) + (1|Block), Multilocation)) print(anova(fm1MultR, fm2MultR)) ## Not run: intervals(fm1MultR) ## End(Not run) }