1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
\name{summarizeRobWeights}
\alias{summarizeRobWeights}
\title{Print a Nice "summary" of Robustness Weights}
\description{
Print a nice \dQuote{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, \dots)
}
\arguments{
\item{w}{numeric vector of robustness weigths.}
\item{digits}{digits to be used for \code{\link{print}}ing.}
\item{header}{string to be printed as header line.}
\item{eps}{numeric tolerance \eqn{\epsilon}{eps}: values of \code{w}
with \eqn{\left|w_i\right| < \epsilon/n}{|w[i]| < eps/n} are said to
be outliers.}
\item{eps1}{numeric tolerance: values of \code{w} with
\eqn{\left|1 - w_i\right| < eps1}{|1 - w[i]| < eps1} are said to
have weight \sQuote{\code{~= 1}}.}
\item{\dots}{potential further arguments, passed to
\code{\link{print}()}.}
}
\seealso{
The \code{\link{summary}} methods for \code{\link{lmrob}}
and \code{\link{glmrob}} make use of \code{summarizeRobWeights()}.
Our methods for \code{\link{weights}()},
\code{\link{weights.lmrob}(*, type="robustness")} and
\code{\link{weights.glmrob}(*, type="robustness")}.
}
\value{
none; the function is used for its side effect of printing.
}
\author{Martin Maechler}
\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
}
\keyword{utilities}
|
