Description Usage Arguments Details Value Author(s) References See Also Examples
Draws a contour plot of Tukey's depth function.
1 2 3 4 |
x |
Bivariate data as a matrix, data frame or list. If it is a matrix or data frame, then each row is viewed as one bivariate observation. If it is a list, both components must be numerical vectors of equal length (coordinates of observations). |
dpth |
Vector of positive integers. Numbers 1,2,… refer to contours of depth 1/n, 2/n,…, where n is the number of observations. Useful to draw particular contours. Default |
output |
Logical. Default FALSE produces a contour plot; otherwise a list of contour vertices. |
twodim |
Logical. |
mustdith |
Logical. Should dithering be applied? Used when data set is not in general position or a numerical problem is encountered. |
maxdith |
Positive integer. Maximum number of dithering steps. |
dithfactor |
Scaling factor used for horizontal and vertical dithering. |
trace.errors |
Logical. Should all contours be considered? Used when a numerical problem is encountered for some inner contours. Default |
eps |
Error tolerance to control the calculation. |
factor |
Proportion (0 to 1) of outermost contours computed according to a version of the algorithm ISODEPTH of Rousseeuw and Ruts (1998); remaining contours are derived from an algorithm in Rousseeuw et al. (1999). |
xlab |
Title for x-axis. Must be a character string. |
ylab |
Title for y-axis. Must be a character string. |
zlab |
Title for z-axis. Used jointly with |
colcontours |
Vector of color names of some or all of the contours. Recycling is used when necessary. Colors can be specified in different ways, see color specification in |
... |
Any additional graphical parameters (see |
Tukey's depth and dimension 2 only. Contours are computed according to algorithm
ISODEPTH by Ruts and Rousseeuw (1996) or, more exactly, revised versions of this
algorithm which appear in Rousseeuw and Ruts (1998) and Rousseeuw et al.
(1999). Argument factor
determines which version to use. If n is
the number of observations, contours of depth ≤ factor
n/2 are obtained from the 1998 version, while the remaining contours are derived from the 1999 version.
When the data set is not in general position, dithering can be used in the sense that random noise is added to each component of each observation. Random noise takes the form eps
times dithfactor
times U for the horizontal component and eps
times dithfactor
times V for the vertical component, where U, V are independent uniform on [-.5, 5.]. This is done in a number of consecutive steps applying independent U's and V's.
Default output = FALSE
yields a contour plot. If not, the function returns a list of m components, where m is the number of contours and component i is a matrix whose rows are the vertices of contour i.
Jean-Claude Masse and Jean-Francois Plante, based on Fortran code by Rousseeuw, Ruts and Struyf from University of Antwerp.
Ruts, I. and Rousseeuw, P.J. (1996), Computing depth contours of bivariate point clouds, Comput. Stat. Data An., 23, 153–168.
Rousseeuw, P.J. and Ruts, I. (1998), Constructing the bivariate Tukey median, Stat. Sinica, 8, 828–839.
Rousseeuw, P.J., Ruts, I., and Tukey, J.W. (1999), The Bagplot: A Bivariate Boxplot, The Am. Stat., 53, 382–387.
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 | ## exact contour plot with 10 contours
set.seed(601) ; x = matrix(rnorm(48), nc = 2)
isodepth(x)
## exact colored contours
set.seed(159); library(MASS)
mu1 <- c(0,0); mu2 <- c(6,0); sigma <- matrix(c(1,0,0,1), nc = 2)
mixbivnorm <- rbind(mvrnorm(80, mu1 ,sigma), mvrnorm(20, mu2, sigma))
isodepth(mixbivnorm, dpth = c(35,5), col = rainbow(2))
## vertices of each contour
set.seed(601)
x <- matrix(rnorm(48), nc = 2)
isodepth(x, output = TRUE)
## data set not in general position
data(starsCYG, package = "robustbase")
isodepth(starsCYG, mustdith = TRUE)
## colored contours
set.seed(601)
x <- matrix(rnorm(48), nc = 2)
isodepth(x, colcontours= rainbow(10))
# perspective plot
library(rgl)
set.seed(601)
x <- matrix(rnorm(48), nc = 2)
isodepth(x, twodim = FALSE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.