multiecdf {geneplotter}R Documentation

Multiple empirical cumulative distribution functions (ecdf) and densities

Description

Plot multiple empirical cumulative distribution functions (ecdf) and densities with user interface similar to that of boxplot.

Usage

multiecdf(x, ...)
## S3 method for class 'formula':
multiecdf(formula, data = NULL, xlab, na.action = NULL, ...)
## S3 method for class 'matrix':
multiecdf(x, xlab, ...) 
## Default S3 method:
multiecdf(x, xlim,
       col=brewer.pal(9, "Set1"),
       main="multiecdf",
       xlab,
       do.points=FALSE,
       subsample=TRUE, ...)

multidensity(x, ...)
## S3 method for class 'formula':
multidensity(formula, data = NULL, xlab, na.action = NULL, ...)
## S3 method for class 'matrix':
multidensity(x, xlab, ...) 
## Default S3 method:
multidensity(x,
  bw="nrd0",
  xlim,
  ylim,
  col  = brewer.pal(9, "Set1"),
  main = "multidensity",
  xlab,
  lty  = 1L , ...)

Arguments

formula a formula, such as y ~ grp, where y is a numeric vector of data values to be split into groups according to the grouping variable grp (usually a factor).
data a data.frame (or list) from which the variables in formula should be taken.
na.action a function which indicates what should happen when the data contain NAs. The default is to ignore missing values in either the response or the group.
x a list of numeric vectors.
bw the smoothing bandwidth, see the manual page for density. If bw is a character string specifying a rule to choose the bandwidth, this rule is applied to x[[1]] and then the same numerical value of bw is used throughout.
xlim Range of the x axis. If missing, the data range is used.
ylim Range of the y axis. If missing, the range of the density estimates is used.
col, lty Line colors and line type.
main Plot title.
xlab x-axis label.
do.points logical; if TRUE, also draw points at the knot locations.
subsample logical; if TRUE, subsamples of size 1000 are used to compute and plot the ecdf for list items with many observations (>1000)
... Further arguments that get passed on to the plot functions.

Details

The usefulness of multidensity can vary, depending on the data and because of smoothing artifacts. multiecdf will in many cases be preferable.

Value

For the multidensity functions, a list of density objects.

Author(s)

Wolfgang Huber http://www.ebi.ac.uk/huber

See Also

boxplot, ecdf density

Examples

  f = 1 + (runif(1000)>0.5)
  x = rnorm(length(f), mean=f, sd=f)
  
  multiecdf(x~f)
  multidensity(x~f)

[Package geneplotter version 1.20.0 Index]