oats                  package:MASS                  R Documentation

_D_a_t_a _f_r_o_m _a_n _O_a_t_s _F_i_e_l_d _T_r_i_a_l

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

     The yield of oats from a split-plot field trial using three
     varieties and four levels of manurial treatment.  The experiment
     was laid out in 6 blocks of 3 main plots, each split into 4
     sub-plots.  The varieties were applied to the main plots and the
     manurial treatments to the sub-plots.

_U_s_a_g_e:

     data(oats)

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

     This data frame contains the following columns:

     `_B' Blocks, levels I, II, III, IV, V and VI

     `_V' Varieties, 3 levels.

     `_N' Nitrogen (manurial) treatment, levels  0.0cwt, 0.2cwt, 0.4cwt
          and 0.6cwt, showing the application in cwt/acre.

     `_Y' Yields in 1/4lbs per sub-plot, each of area 1/80 acre.

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

     Yates, F. (1935) Complex experiments, Journal of the Royal
     Statistical Society Suppl. 2, 181-247.

     Also given in  Yates, F. (1970)  Experimental design: Selected
     papers of Frank Yates, C.B.E, F.R.S. London: Griffin.

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

     data(oats)
     oats$Nf <- ordered(oats$N, levels=sort(levels(oats$N)))
     oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = TRUE)
     summary(oats.aov)
     # next line not implemented in R
     summary(oats.aov, split=list(Nf=list(L=1, Dev=2:3)))
     par(mfrow=c(1,2), pty="s")
     plot(fitted(oats.aov[[4]]), studres(oats.aov[[4]]))
     abline(h=0, lty=2)
     oats.pr <- proj(oats.aov)
     qqnorm(oats.pr[[4]][,"Residuals"], ylab="Stratum 4 residuals")
     qqline(oats.pr[[4]][,"Residuals"])

     par(mfrow=c(1,1), pty="m")
     oats.aov2 <- aov(Y ~ N + V + Error(B/V), data = oats, qr = TRUE)
     model.tables(oats.aov2, type = "means", se = TRUE)

