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, main, xlab, ...)
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 |
Line colors. If missing, the
Set1 palette from RColorBrewer
is used. |
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)