ave                   package:base                   R Documentation

_G_r_o_u_p _A_v_e_r_a_g_e_s _O_v_e_r _L_e_v_e_l _C_o_m_b_i_n_a_t_i_o_n_s _o_f _F_a_c_t_o_r_s

_D_e_s_c_r_i_p_t_i_o_n:

     Subsets of `x[]' are averaged, where each subset consist of those
     observations with the same factor levels.

_U_s_a_g_e:

     ave(x, ..., FUN = mean)

_A_r_g_u_m_e_n_t_s:

       x: A numeric.

     ...: Grouping variables, typically factors, all of the same
          `length' as `x'.

     FUN: Function to apply for each factor level combination.

_V_a_l_u_e:

     A numeric vector, say `y' of length `length(x)'. If `...' is
     `g1,g2', e.g., `y[i]' is equal to `FUN(x[j]', for all `j' with
     `g1[j]==g1[i]' and `g2[j]==g2[i])'.

_S_e_e _A_l_s_o:

     `mean', `median'.

_E_x_a_m_p_l_e_s:

     ave(1:3)# no grouping -> grand mean

     data(warpbreaks)
     attach(warpbreaks)
     ave(breaks, wool)
     ave(breaks, tension)
     ave(breaks, tension, FUN = function(x)mean(x, trim=.1))
     plot(breaks, main =
          "ave( Warpbreaks )  for   wool  x  tension  combinations")
     lines(ave(breaks, wool, tension            ), type='s', col = "blue")
     lines(ave(breaks, wool, tension, FUN=median), type='s', col = "green")
     legend(40,70, c("mean","median"), lty=1,col=c("blue","green"), bg="gray90")
     detach()

