

mad {base}                                   R Documentation

_M_e_d_i_a_n _A_b_s_o_l_u_t_e _D_e_v_i_a_t_i_o_n

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

     Compute a scale estimate based on the median absolute
     deviation.

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

     mad(x, center, constant = 1.4826, na.rm = FALSE)

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

       x: a numeric vector.

  center: Optionally, the centre: defauls to the median.

constant: scale factor.

   na.rm: if `TRUE' then `NA' values are stripped from `x'
          before computation takes place.

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

     The actual value calculated is `constant *
     (median(abs(x - center)))' with the default value of
     `center' being `median(x)'.

     The default `constant = 1.4826' (approximately 1/
     Phi^(-1)(3/4) = `1/qnorm(3/4)') ensures consistency,
     i.e.,

                   E[mad(X_1,...,X_n)] = sigma

     for X_i distributed as N(mu,sigma^2) and large n.

     If `na.rm' is `TRUE' then `NA' values are stripped from
     `x' before computation takes place.  If this is not
     done then an `NA' value in `x' will cause `mad' to
     return `NA'.

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

     `IQR' which is simpler but less robust, `median',
     `var'.

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

     mad(c(1:9))
     print(mad(c(1:9),     constant=1)) ==
           mad(c(1:8,100), constant=1)       # = 2 ; TRUE

