summarizeRobWeights: Print a Nice "summary" of Robustness Weights

In robustbase: Basic Robust Statistics

Print a Nice "summary" of Robustness Weights

Description

Print a nice “summary” about a numeric vector of robustness weights. Observations with weights around zero are marked as outliers.

Usage

summarizeRobWeights(w, digits = getOption("digits"),
                    header = "Robustness weights:",
 	            eps = 0.1 / length(w), eps1 = 1e-3, ...)

Arguments

w	numeric vector of robustness weigths.
digits	digits to be used for printing.
header	string to be printed as header line.
eps	numeric tolerance \epsilon: values of w with \left|w_i\right| < \epsilon/n are said to be outliers.
eps1	numeric tolerance: values of w with \left|1 - w_i\right| < eps1 are said to have weight ‘~= 1’.
...	potential further arguments, passed to print().

Value

none; the function is used for its side effect of printing.

Author(s)

Martin Maechler

See Also

The summary methods for lmrob and glmrob make use of summarizeRobWeights().

Our methods for weights(), weights.lmrob(*, type="robustness") and weights.glmrob(*, type="robustness").

Examples

w <- c(1,1,1,1,0,1,1,1,1,0,1,1,.9999,.99999, .5,.6,1e-12)
summarizeRobWeights(w) # two outside ~= {0,1}
summarizeRobWeights(w, eps1 = 5e-5)# now three outside {0,1}

## See the summary(<lmrob>) outputs

robustbase documentation built on Sept. 27, 2024, 5:09 p.m.