R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Approximate depth contours

depthContour

R Documentation

Approximate depth contours

Description

Draws an approximate contours of depth for bivariate data.

Usage

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{ zin {{{R}}^{d}}:D(z,F)≥ α ight}

Author(s)

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

Examples

# 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")

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(DepthProc)
Loading required package: ggplot2
Loading required package: Rcpp
Loading required package: rrcov
Loading required package: robustbase
Scalable Robust Estimators with High Breakdown Point (version 1.3-11)

Loading required package: MASS
Loading required package: np
Nonparametric Kernel Methods for Mixed Datatypes (version 0.60-2)
[vignette("np_faq",package="np") provides answers to frequently asked questions]
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/DepthProc/depthContour.Rd_%03d_medium.png", width=480, height=480)
> ### Name: depthContour
> ### Title: Approximate depth contours
> ### Aliases: depthContour
> ### Keywords: contour depth
> 
> ### ** Examples
> 
> # 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")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>