| multiecdf {geneplotter} | R Documentation |
Plot multiple empirical cumulative distribution functions (ecdf)
and densities with user interface similar to that of boxplot.
multiecdf(x, ...)
## S3 method for class 'formula':
multiecdf(formula, data = NULL, ..., na.action = NULL)
## Default S3 method:
multiecdf(x, xlim, col, do.points=FALSE, subsample=TRUE, ...)
multidensity(x, ...)
## S3 method for class 'formula':
multidensity(formula, data = NULL, main, xlab, ..., na.action = NULL)
## Default S3 method:
multidensity(x, bw="nrd0", xlim, ylim,
col = brewer.pal(9, "Set1"),
main = sprintf("multidensity.default(%s)", deparse(substitute(x))),
xlab = deparse(substitute(x)),
lty = 1L , ...)
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. |
The usefulness of multidensity can vary,
depending on the data and because of smoothing
artifacts. multiecdf will in many cases be preferable.
For the multidensity functions, a list of density objects.
Wolfgang Huber http://www.ebi.ac.uk/huber
f = 1 + (runif(1000)>0.5) x = rnorm(length(f), mean=f, sd=f) multiecdf(x~f) multidensity(x~f)