Multilocation            package:SASmixed            R Documentation

_A _m_u_l_t_i_l_o_c_a_t_i_o_n _t_r_i_a_l

_D_e_s_c_r_i_p_t_i_o_n:

     The 'Multilocation' data frame has 108 rows and 7 columns.

_F_o_r_m_a_t:

     This data frame contains the following columns:

     _o_b_s a numeric vector

     _L_o_c_a_t_i_o_n an ordered factor with levels 'B' < 'D' < 'E' < 'I' < 'G'
          < 'A' < 'C' < 'F' < 'H'

     _B_l_o_c_k a factor with levels '1' to '3' 

     _T_r_t a factor with levels '1' to '4'

     _A_d_j a numeric vector

     _F_e a numeric vector

     _G_r_p an ordered factor with levels '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'

_S_o_u_r_c_e:

     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).

_E_x_a_m_p_l_e_s:

     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)
     }

