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 `NA` s. 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]