

   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

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

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

        `mad' provides a scale estimate based on the median
        absolute deviation.  The actual value calculated is
        `constant * (median(abs(x - center)))' with the default
        value of `center' being `median(x)'.

        The default `constant = 1.4826' is pprox Phi^{-1}(ac 3
        4)= 1.4826 ~=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::

        `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

