Print a Nice "summary" of Robustness Weights
Print a nice “summary” about a numeric vector of robustness weights. Observations with weights around zero are marked as outliers.
summarizeRobWeights(w, digits = getOption("digits"),
header = "Robustness weights:",
eps = 0.1 / length(w), eps1 = 1e-3, ...)w |
numeric vector of robustness weigths. |
digits |
digits to be used for |
header |
string to be printed as header line. |
eps |
numeric tolerance eps: values of |
eps1 |
numeric tolerance: values of |
... |
potential further arguments, passed to
|
none; the function is used for its side effect of printing.
Martin Maechler
Our methods for weights(),
weights.lmrob(*, type="robustness") and
weights.glmrob(*, type="robustness").
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>) outputsPlease choose more modern alternatives, such as Google Chrome or Mozilla Firefox.