compareDistributions: Compare two density distributions side by side

Description Usage Arguments Value Author(s) Examples

Description

This function will produce a plot of two density functions displayed side by side

Usage

1
2
3
compareDistributions(left, right, add.spread = TRUE,
  print.stats = TRUE, xlim = NULL, ylim = NULL, clrs = c("purple",
  "lightblue"), xlab = "density", ylab = "value", ...)

Arguments

left

numeric vector.

right

numeric vector.

add.spread

logical, whether to plot the spread (q25 to q75 and the median).

print.stats

logical, whether to print summary statistics for each distribution.

xlim, ylim

numeric axis limits, see xyplot.

clrs

A character vector of length 2 specifying the colors for the filled density regions.

xlab, ylab

character axis labels, see plot.

...

Additional arguments passed to density.

Value

A trellis object.

Author(s)

Tim Appelhans

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
compareDistributions(rnorm(1000, 2, 3), rnorm(1000, -5, 1))
compareDistributions(rnorm(1000, 2, 3), rnorm(1000, -5, 1), 
                     add.spread = FALSE)
compareDistributions(rnorm(1000, 2, 3), rnorm(1000, -5, 1), 
                     add.spread = TRUE, clrs = c("red", "brown"))
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = FALSE, add.spread = FALSE)
                     
## pass additional parameters to density() 
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = FALSE, add.spread = FALSE, bw = 5)
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = FALSE, add.spread = FALSE, bw = 8,
                     kernel = "rectangular")
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = FALSE, add.spread = TRUE, bw = 8,
                     n = 3)
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = TRUE, add.spread = TRUE, bw = 0.1)
compareDistributions(rnorm(1000, 2, 5), rnorm(1000, -5, 4), 
                     print.stats = TRUE, add.spread = TRUE, bw = 0.5)                    

Orcs documentation built on May 2, 2019, 2:37 a.m.