depthContour: Approximate depth contours

Description Usage Arguments Details Author(s) See Also Examples

Description

Draws an approximate contours of depth for bivariate data.

Usage

1
2
3
4
depthContour(x, xlim = extendrange(x[, 1], f = 0.1), ylim = extendrange(x[,
  2], f = 0.1), n = 50, pmean = TRUE, mcol = "blue", pdmedian = TRUE,
  mecol = "brown", legend = TRUE, points = FALSE, colors = heat_hcl,
  levels = 10, ...)

Arguments

x

Bivariate data

xlim

Determines the width of x-axis.

ylim

Determines the width of y-axis.

n

Number of points in each coordinate direction to be used in contour plot.

pmean

Logical. If TRUE mean will be marked.

mcol

Determines the color of lines describing the mean.

pdmedian

Logical. If TRUE depth median will be marked.

mecol

Determines the color of lines describing the depth median.

legend

Logical. If TRUE legend for mean and depth median will be drawn.

points

Logical. If TRUE points from matrix x will be drawn.

colors

function for colors pallete (e.g. gray.colors).

levels

number of levels for color scale.

...

Any additional parameters for function depth (such as method) or graphical parameters (e.g. lwd, lty, main).

Details

The set of all points that have depth at least α is called α - trimmed region. The α - trimmed region w.r.t. F is denoted by {D}_{α }(F) , i.e.,

{D}_{α }(F)=≤ft\{ z\in {{{R}}^{d}}:D(z,F)≥ α \right\}

.

Author(s)

Daniel Kosiorowski, Mateusz Bocian, Anna Wegrzynkiewicz and Zygmunt Zawadzki from Cracow University of Economics.

See Also

depthPersp

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# EXAMPLE 1
x = mvrnorm(1000,c(0,0),diag(2))
depthContour(x, colors = gray.colors)
#with points
depthContour(x, points = TRUE)
depthContour(x, points = FALSE, levels = 10)

 #EXAMPLE 2
 data(inf.mort,maesles.imm)
 data1990=na.omit(cbind(inf.mort[,1],maesles.imm[,1]))
 depthContour(data1990, n = 50, pmean = TRUE, mcol = "blue",
 pdmedian = TRUE, mecol = "brown", legend = TRUE, points = TRUE,
 xlab="infant mortality rate per 1000 live birth",
 ylab="against masles immunized percentage", main='L2 depth,
 UN Fourth Goal 2011 year',method = "LP")

DepthProc documentation built on May 2, 2019, 6:22 p.m.