

npk(MASS)                                    R Documentation

_C_l_a_s_s_i_c_a_l _N_, _P_, _K _F_a_c_t_o_r_i_a_l _E_x_p_e_r_i_m_e_n_t

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

     A classical N, P, K (nitrogen, phosphate, potassium)
     factorial experiment on the growth of peas conducted on
     6 blocks. Each half of a fractional factorial design
     confounding the NPK interaction was used on 3 of the
     plots.

_A_r_g_u_m_e_n_t_s_:

   block: which block (label 1 to 6).

       N: indicator (0/1) for the application of nitrogen.

       P: indicator (0/1) for the application of phosphate.

       K: indicator (0/1) for the application of potassium.

   yield: Yield of peas, in pounds/plot (the plots were
          (1/70) acre).

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

     The `npk' data frame has 24 rows and 5 columns:

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

     Imperial College, London, M.Sc. exercise sheet.

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

     data(npk)
     npk.aov <- aov(yield ~ block + N*P*K, npk)
     npk.aov
     summary(npk.aov)
     alias(npk.aov)
     coef(npk.aov)
     options(contrasts=c("contr.treatment", "contr.poly"))
     npk.aov1 <- aov(yield ~ block + N + K, npk)
     summary.lm(npk.aov1)
     se.contrast(npk.aov1, list(N=="0", N=="1"), data=npk)
     model.tables(npk.aov1, type="means", se=TRUE)

