

medpolish {eda}                              R Documentation

_M_e_d_i_a_n _P_o_l_i_s_h _o_f _a _M_a_t_r_i_x

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

     Fits an additive model using Tukey's median polish pro-
     cedure.

_U_s_a_g_e_:

     medpolish(x, eps = 0.01, maxiter = 10, trace.iter = TRUE)

     plot(medpolish.obj)
     print(medpolish.obj)

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

       x: a numeric matrix.

     eps: real number greater than 0. A tolerance for con-
          vergence: see Details.

 maxiter: the maximum number of iterations

trace.iter: logical. Should progress in convergence be
          reported?

medpolish.obj: object of class `medpolish'.

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

     The model fitted is additive (constant + rows +
     columns). The algorithm works by alternately removing
     the row and column medians, and continues until the
     proportional reduction in the sum of absolute residuals
     is less than `eps' or until there have been `maxiter'
     iterations.  The sum of absolute residuals is printed
     at each iteration of the fitting process, if
     `trace.iter' is `TRUE'.

     `medpolish' returns an object of class `medpolish' (see
     below).  There are printing and plotting methods for
     this class, which are invoked via by the generics
     `print' and `plot'.

_V_a_l_u_e_:

     An object of class `medpolish' with the following named
     components:

 overall: the fitted constant term.

     row: the fitted row effects.

     col: the fitted column effects.

residuals: the residuals.

    name: the name of the dataset.

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

     Tukey, J. W. (1977).  Exploratory Data Analysis, Read-
     ing Massachusetts: Addison-Wesley.

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

     `median'; `aov' for a mean instead of median decomposi-
     tion.

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

     ## Deaths from sport parachuting;  from ABC of EDA, p.224:
     deaths <-
         rbind(c(14,15,14),
               c( 7, 4, 7),
               c( 8, 2,10),
               c(15, 9,10),
               c( 0, 2, 0))
     dimnames(deaths) <- list(c("1-24", "25-74", "75-199", "200++", "NA"),
                              paste(1973:1975))
     deaths
     (med.d <- medpolish(deaths))
     plot(med.d)
     ## Check decomposition:
     all(deaths == med.d$overall + outer(med.d$row,med.d$col, "+") + med.d$resid)

