replications              package:stats              R Documentation

_N_u_m_b_e_r _o_f _R_e_p_l_i_c_a_t_i_o_n_s _o_f _T_e_r_m_s

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

     Returns a vector or a list of the number of replicates for each
     term in the formula.

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

     replications(formula, data=NULL, na.action)

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

 formula: a formula or a terms object or a data frame.

    data: a data frame used  to  find  the  objects in 'formula'.

na.action: function for handling missing values.  Defaults to a
          'na.action' attribute of 'data', then a setting of the option
          'na.action', or 'na.fail' if that is not set.

_D_e_t_a_i_l_s:

     If 'formula' is a data frame and 'data' is missing, 'formula' is
     used for 'data' with the formula '~ .'.

_V_a_l_u_e:

     A vector or list with one entry for each term in the formula
     giving the number(s) of replications for each level. If all levels
     are balanced (have the same number of replications) the result is
     a vector, otherwise it is a list with a component for each terms,
     as a vector, matrix or array as required.

     A test for balance is '!is.list(replications(formula,data))'.

_A_u_t_h_o_r(_s):

     The design was inspired by the S function of the same name
     described in Chambers _et al._ (1992).

_R_e_f_e_r_e_n_c_e_s:

     Chambers, J. M., Freeny, A and Heiberger, R. M. (1992) _Analysis
     of variance; designed experiments._ Chapter 5 of _Statistical
     Models in S_ eds J. M. Chambers and T. J. Hastie, Wadsworth &
     Brooks/Cole.

_S_e_e _A_l_s_o:

     'model.tables'

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

     ## From Venables and Ripley (2002) p.165.
     N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0)
     P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0)
     K <- c(1,0,0,1,0,1,1,0,0,1,0,1,0,1,1,0,0,0,1,1,1,0,1,0)
     yield <- c(49.5,62.8,46.8,57.0,59.8,58.5,55.5,56.0,62.8,55.8,69.5,
     55.0, 62.0,48.8,45.5,44.2,52.0,51.5,49.8,48.8,57.2,59.0,53.2,56.0)

     npk <- data.frame(block=gl(6,4), N=factor(N), P=factor(P),
                       K=factor(K), yield=yield)
     replications(~ . - yield, npk)

