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
|
\name{radarImage}
\alias{radarImage}
\docType{data}
\title{Satellite Radar Image Data from near Munich}
\description{
The data were supplied by A. Frery. They are a part of a synthetic aperture
satellite radar image corresponding to a suburb of Munich. Provided
are coordinates and values corresponding to three frequency bands for
each of 1573 pixels.
}
\usage{data(radarImage, package="robustbase")}
\format{
A data frame with 1573 observations on the following 5 variables.
\describe{
\item{\code{X.coord}}{a numeric vector}
\item{\code{Y.coord}}{a numeric vector}
\item{\code{Band.1}}{a numeric vector}
\item{\code{Band.2}}{a numeric vector}
\item{\code{Band.3}}{a numeric vector}
}
}
% \details{
% }
\source{
The website accompanying the MMY-book:
\url{https://www.wiley.com/legacy/wileychi/robust_statistics/}
}
\examples{
data(radarImage)
plot(Y.coord ~ X.coord, data = radarImage)
## The 8 "clear" outliers (see also below)
ii8 <- c(1548:1549, 1553:1554, 1565:1566, 1570:1571)
outF <- 1+(seq_len(nrow(radarImage)) \%in\% ii8)
pairs(radarImage[, 3:5], main = "radarImage (n = 1573)",
col = outF, pch=outF)
## Finding outliers -----------------------------------------
set.seed(1)
system.time(cc.ri <- covMcd(radarImage))# ~ 0.1 sec
## check for covMcd() consistency:
iiO <- as.integer(
c(262, 450:451, 480:481, 509, 535, 542, 597, 643, 669, 697, 803:804, 832:834,
862:864, 892, 989, 1123, 1145, 1223:1224, 1232:1233, 1249:1250, 1267, 1303,
1347, 1357, 1375, 1411, 1419:1420, 1443, 1453, 1504, 1510:1512,
1518:1521, 1525:1526, 1543:1544, 1546:1555, 1557:1558, 1561:1562, 1564:1566,
1569:1571, 1573))
length(iiO) # 73 -- other seeds sometimes give 72, rarely 71 "outliers"
table(isO <- cc.ri$mcd.wt == 0) # 2023-05: 118
stopifnot(exprs = {
## identical(iiO, which(isO)) -- TRUE before 2023-05 covMcd() change
ii8 \%in\% which(isO) # ii8 is subset of isO
identical(ii8, which(cc.ri$mah > 200))
length(intersect(cc.ri$best, iiO)) == 0
})
cc <- c(adjustcolor("black", 0.4), adjustcolor("tomato", 0.8))
pairs(radarImage, main = "radarImage (n = 1573) + Outliers", gap=0,
col = cc[1+isO], pch = c(1,8)[1+isO], cex = 0.8)
}
\keyword{datasets}
|
