DescTools (version 0.99.11)

PlotMultiDens: Plot Multiple Density Curves

Description

Multiple density curves are plotted on the same plot. The function plots the density curves in the defined colors and linetypes, after having calculated the globally appropriate xlim- and ylim-values. A legend can directly be included.

Usage

PlotMultiDens(x, ...)

## S3 method for class 'default':
PlotMultiDens(x, xlim = NULL, ylim = NULL,
              col = rainbow(length(x)),
              lty = "solid", lwd = 1, xlab = "x", ylab = "density",
              args.dens = NULL, args.legend = NULL,
              na.rm = FALSE, flipxy=FALSE, ...)

## S3 method for class 'formula':
PlotMultiDens(formula, data, subset, \dots)

Arguments

x
a list of vectors whose densities are to be plotted. Usesplit to separate a vector by groups. (See examples)
xlim, ylim
xlim, ylim of the plot.
col
colors of the lines, defaults to rainbow(1:length(x)).
lty
line type of the lines.
lwd
line widths for the lines.
xlab, ylab
a title for the x, resp. y axis. Defaults to "x" and "density".
args.dens
list of additional arguments to be passed to the density function. If set to NULL the defaults will be used. Those are n = 4096 (2^12) and kernel = "epanechnikov".
args.legend
list of additional arguments to be passed to the legend function. Use args.legend = NA if no legend should be added.
na.rm
should NAs be omitted? Defaults to FALSE.
flipxy
logical, should x- and y-axis be flipped? Defaults to FALSE.
formula
a formula of the form lhs ~ rhs where lhs gives the data values and rhs the corresponding groups.
data
an optional matrix or data frame (or similar: see model.frame) containing the variables in the formula formula. By default the variables are taken from environment(formula)
subset
an optional vector specifying a subset of observations to be used.
...
the dots are passed to plot(...).

Value

  • data.frame with 3 columns, containing the bw, n and kernel parameters used for the list elements. The number of rows correspond to the length of the list x.

Details

All style arguments, density arguments and data list elements will be recycled if necessary. The flipxy parameter flips x an y-values, so as to plot density curves on the x-axis.

See Also

PlotViolin, density

Examples

Run this code
x <- rnorm(1000,0,1)
y <- rnorm(1000,0,2)
z <- rnorm(1000,2,1.5)

# the input of the following function MUST be a numeric list
PlotMultiDens(list(x=x,y=y,z=z))


PlotMultiDens( x=split(d.pizza$delivery_min, d.pizza$driver), na.rm=TRUE
  , main="delivery time ~ driver", xlab="delivery time [min]", ylab="density"
  , lwd=1:7, lty=1:7
  , panel.first=grid())
# this example demonstrates the definition of different line types and -colors
# an is NOT thought as recommendation for good plotting practice... :-)


# the formula interface
PlotMultiDens(delivery_min ~ driver, data=d.pizza)

# recyling of the density parameters
res <- PlotMultiDens(x=split(d.pizza$temperature, d.pizza$driver), 
              args.dens = list(bw=c(5,2), kernel=c("rect","epanechnikov")), na.rm=TRUE)
res

# compare bandwidths
PlotMultiDens(x=split(d.pizza$temperature, d.pizza$driver)[1], 
                     args.dens = list(bw=c(1:5)), na.rm=TRUE,
                     args.legend=NA, main="Compare bw")
legend(x="topright", legend=gettextf("bw = %s", 1:5), fill=rainbow(5))

Run the code above in your browser using DataCamp Workspace