

   mad {base}                                   R Documentation

   MMeeddiiaann AAbbssoolluuttee DDeevviiaattiioonn

   DDeessccrriippttiioonn::

        Compute a scale estimate based on the median absolute
        deviation.

   UUssaaggee::

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

   AArrgguummeenntts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Deettaaiills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Seeee AAllssoo::

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

   EExxaammpplleess::

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

