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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
\name{tukeyPsi1}
\alias{robustbase-deprecated}
\alias{tukeyPsi1}
\alias{tukeyChi}
\title{Tukey's Bi-square Score (Psi) and "Chi" (Rho) Functions and Derivatives}
\description{
These are \bold{deprecated}, replaced by
\code{\link{Mchi}(*, psi="tukey")}, \code{\link{Mpsi}(*, psi="tukey")}
\code{tukeyPsi1()} computes Tukey's bi-square score (psi) function, its first
derivative or it's integral/\dQuote{principal function}. This is
scaled such that \eqn{\psi'(0) = 1}, i.e.,
\eqn{\psi(x) \approx x}{psi(x) ~= x} around 0.
\code{tukeyChi()} computes Tukey's bi-square loss function,
\code{chi(x)} and its first two derivatives. Note that in the general
context of \eqn{M}-estimators, these loss functions are called
\eqn{\rho (rho)}{rho}-functions.
}
\usage{
tukeyPsi1(x, cc, deriv = 0)
tukeyChi (x, cc, deriv = 0)
}
\arguments{
\item{x}{numeric vector.}
\item{cc}{ tuning constant }
\item{deriv}{integer in \eqn{\{-1,0,1,2\}} specifying the order of the
derivative; the default, \code{deriv = 0} computes the psi-, or
chi- ("rho"-)function.}
}
\value{
a numeric vector of the same length as \code{x}.
}
\note{\code{tukeyPsi1(x, d)} and \code{\link{tukeyChi}(x, d+1)} are just
re-scaled versions of each other (for \code{d in -1:1}), i.e.,
\deqn{\chi^{(\nu)}(x, c) = (6/c^2) \psi^{(\nu-1)}(x,c),}
for \eqn{\nu = 0,1,2}.
We use the name \sQuote{tukeyPsi\bold{1}}, because \code{tukeyPsi} is
reserved for a future \dQuote{Psi Function} class object, see
\code{\link{psiFunc}}.
}
\seealso{
\code{\link{lmrob}} and \code{\link{Mpsi}}; further
\code{\link{anova.lmrob}} which needs the \code{deriv = -1}.
}
\author{Matias Salibian-Barrera, Martin Maechler and Andreas Ruckstuhl}
\examples{
\dontshow{oop <- options(warn = -1) # against the "deprecated" warnings}
op <- par(mfrow = c(3,1), oma = c(0,0, 2, 0),
mgp = c(1.5, 0.6, 0), mar= .1+c(3,4,3,2))
x <- seq(-2.5, 2.5, length = 201)
cc <- 1.55 # as set by default in lmrob.control()
plot. <- function(...) { plot(...); abline(h=0,v=0, col="gray", lty=3)}
plot.(x, tukeyChi(x, cc), type = "l", col = 2)
plot.(x, tukeyChi(x, cc, deriv = 1), type = "l", col = 2)
plot.(x, tukeyChi(x, cc, deriv = 2), type = "l", col = 2)
% \ is escape for Rd
mtext(sprintf("tukeyChi(x, c = \%g, deriv), deriv = 0,1,2", cc),
outer = TRUE, font = par("font.main"), cex = par("cex.main"))
par(op)
op <- par(mfrow = c(3,1), oma = c(0,0, 2, 0),
mgp = c(1.5, 0.6, 0), mar= .1+c(3,4,1,1))
x <- seq(-5, 5, length = 201)
cc <- 4.69 # as set by default in lmrob.control()
plot. <- function(...) { plot(..., asp = 1); abline(h=0,v=0, col="gray", lty=3)}
plot.(x, tukeyPsi1(x, cc), type = "l", col = 2)
abline(0:1, lty = 3, col = "light blue")
plot.(x, tukeyPsi1(x, cc, deriv = -1), type = "l", col = 2)
plot.(x, tukeyPsi1(x, cc, deriv = 1), type = "l", col = 2); abline(h=1,lty=3)
% \ is escape for Rd
mtext(sprintf("tukeyPsi1(x, c = \%g, deriv), deriv = 0, -1, 1", cc),
outer = TRUE, font = par("font.main"), cex = par("cex.main"))
par(op)
\dontshow{options(oop)}
}
\keyword{robust}
|
