Description Usage Arguments Details Value Note Author(s) References See Also Examples
Plots beans to compare the distributions of different groups; it draws one bean per group of data. A bean consists of a one-dimensional scatter plot, its distribution as a density shape and an average line for the distribution. Next to that, an overall average for the whole plot is drawn per default.
Refine from beanplot
source code by using multiple color
schemes from R package RColorBrewer
.
1 2 3 4 5 6 7 8 | beanplot2(..., bw = "SJ-dpi", kernel = "gaussian", cut = 3, cutmin = -Inf,
cutmax = Inf, grownage = 10, what = c(FALSE, TRUE, TRUE, TRUE),
add = FALSE, col, axes = TRUE, log = "auto", handlelog = NA,
ll = 0.16, wd = NA, maxwidth = 0.8, maxstripline = 0.96,
method = "stack", names, overallline = "mean", beanlines = overallline,
horizontal = FALSE, side = "no", jitter = NULL, beanlinewd = 2,
frame.plot = TRUE, border = NA, innerborder = NA, at = NULL,
boxwex = 0.8, ylim = NULL, xlim = NULL, show.names = NA)
|
... |
data which to perform the beanplot on. This data can consist of dataframes, vectors and/or formulas. For each formula, a dataset can be specified with data=[dataset], and a subset can be specified with subset=[subset]. If subset/data arguments are passed, but there are not enough subset/data arguments, they are reused. Additionally, na.action, drop.unused.levels and xlev can be passed to |
bw |
the bandwidth (method) being used, used by |
kernel |
see |
cut |
the beans are cut beyond |
cutmin |
the low-ends of the beans are cut below |
cutmax |
the high-ends of the beans are cut beyond |
grownage |
the width of a bean grows linearly with the count of points, until |
what |
a vector of four booleans describing what to plot. In the following order, these booleans stand for the total average line, the beans, the bean average, and the beanlines. For example, |
add |
if true, do not start a new plot |
col |
the colors to be used. A vector of up to four colors can be used. In the following order, these colors stand for the area of the beans (without the border, use |
axes |
if false, no axes are drawn. |
log |
use |
handlelog |
if |
ll |
the length of the beanline per point found. |
wd |
the linear transformation that determines the width of the beans. By default determined using |
maxwidth |
the maximum width of a bean. |
maxstripline |
the maximum length of a beanline. |
method |
the method used when two points on a bean are the same. |
names |
a vector of names for the groups. |
overallline |
the method used for determining the overall line. Defaults to |
beanlines |
the method used for determining the average bean line(s). Defaults to |
horizontal |
if true, the beanplot is horizontal |
side |
the side on which the beans are plot. Default is |
jitter |
passed to |
beanlinewd |
the width used for the average bean line |
frame.plot |
if true, plots a frame by calling |
border |
the color for the border around a bean. |
innerborder |
a color (vector) for the border inside the bean(s). Especially useful if |
at |
the positions at which a bean should be drawn. |
boxwex |
a scale factor applied to all beans. Compatible with |
ylim |
the range to plot. |
xlim |
the range to plot the beans at. |
show.names |
if true, plots the names as axis labels |
Most parameters are compatible with boxplot
and
stripchart
. For compatibility, arguments with the
name "formula" or "x" are used as data. However, data or formulas
do not need to be named "x" or "formula". The function handles
(combinations of) dataframes, vectors and/or formulas.
bw |
The bandwith (bw) used. |
wd |
The bean width (wd) used. |
In case of more than 5000 values per bean, the autodetection of log fails. In such cases, use
log=""
or log="y"
.
Peter Kampstra <pkampst=beanplot@cs.vu.nl>
Kampstra, P. (2008) Beanplot: A Boxplot Alternative for Visual Comparison of Distributions. Journal of Statistical Software, Code Snippets, 28(1), 1-9. URL http://www.jstatsoft.org/v28/c01/
boxplot
, stripchart
, density
,
rug
, vioplot
in package vioplot
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 | ###############Examples for beanplot2, which provided by Xiangchun LI##############
A=rnorm(100, 1)
B=rnorm(100, 1.3)
C=rnorm(100, 1.5)
D=rnorm(100, 1.6)
E=rnorm(100, 1.8)
F=rnorm(100, 0.8)
data=data.frame(A, B, C, D, E, F)
beanplot2(data)
par(lwd=0.4);beanplot2(data) # Control lwd
beanplot2(data, axes=FALSE, frame.plot=FALSE)
axis(side=1, at=1:6, labels=colnames(data), xaxs='i', lwd=0.4, padj=-1)
axis(side=2, lwd=0.4, yaxs='i', las=1, hadj=0)
box(lwd=0.4)
# If the length of each variables are not equal, use a list as input, e.g.
G=rnorm(200, 1.9)
data=list(A=A,B=B,C=C,D=D,E=E,F=F,G=G)
beanplot2(data)
func.color <- function(color.name) brewer.pal(4,color.name)
col <- lapply(rownames(brewer.pal.info), func.color)
beanplot2(data, col=col[13:18], axes=T, frame.plot=T)
###############End of examples for beanplot2, which provided by Xiangchun LI##############
beanplot(rnorm(22),rnorm(22),rnorm(22),main="Test!",rnorm(3))
#mostly examples taken from boxplot:
par(mfrow = c(1,2))
boxplot(count ~ spray, data = InsectSprays, col = "lightgray")
beanplot(count ~ spray, data = InsectSprays, col = "lightgray", border = "grey", cutmin = 0)
boxplot(count ~ spray, data = InsectSprays, col = "lightgray")
beanplot(count ~ spray, data = InsectSprays, col = "lightgray", border = "grey", overallline = "median")
boxplot(decrease ~ treatment, data = OrchardSprays,
log = "y", col = "bisque", ylim = c(1,200))
beanplot(decrease ~ treatment, data = OrchardSprays,
col = "bisque", ylim = c(1,200))
par(mfrow = c(2,1))
mat <- cbind(Uni05 = (1:100)/21, Norm = rnorm(100),
T5 = rt(100, df = 5), Gam2 = rgamma(100, shape = 2))
par(las=1)# all axis labels horizontal
boxplot(data.frame(mat), main = "boxplot(*, horizontal = TRUE)",
horizontal = TRUE, ylim = c(-5,8))
beanplot(data.frame(mat), main = "beanplot(*, horizontal = TRUE)",
horizontal = TRUE, ylim = c(-5,8))
par(mfrow = c(1,2))
boxplot(len ~ dose, data = ToothGrowth,
boxwex = 0.25, at = 1:3 - 0.2,
subset = supp == "VC", col = "yellow",
main = "Guinea Pigs' Tooth Growth",
xlab = "Vitamin C dose mg",
ylab = "tooth length", ylim = c(-1, 40), yaxs = "i")
boxplot(len ~ dose, data = ToothGrowth, add = TRUE,
boxwex = 0.25, at = 1:3 + 0.2,
subset = supp == "OJ", col = "orange")
legend("bottomright", bty="n", c("Ascorbic acid", "Orange juice"),
fill = c("yellow", "orange"))
allplot <- beanplot(len ~ dose+supp, data = ToothGrowth,
what=c(TRUE,FALSE,FALSE,FALSE),show.names=FALSE,ylim=c(-1,40), yaxs = "i")
beanplot(len ~ dose, data = ToothGrowth, add=TRUE,
boxwex = 0.6, at = 1:3*2 - 0.9,
subset = supp == "VC", col = "yellow",border="yellow2",
main = "Guinea Pigs' Tooth Growth",
xlab = "Vitamin C dose mg",
ylab = "tooth length", ylim = c(3, 40), yaxs = "i",
bw = allplot$bw, wd = allplot$wd, what = c(FALSE,TRUE,TRUE,TRUE))
beanplot(len ~ dose, data = ToothGrowth, add = TRUE,
boxwex = 0.6, at = 1:3*2-0.1,
subset = supp == "OJ", col = "orange",border="darkorange",
bw = allplot$bw, wd = allplot$wd, what = c(FALSE,TRUE,TRUE,TRUE))
legend("bottomright", bty="n", c("Ascorbic acid", "Orange juice"),
fill = c("yellow", "orange"))
par(mfrow = c(1,2))
boxplot(len ~ dose, data = ToothGrowth,
boxwex = 0.25, at = 1:3 - 0.2,
subset = supp == "VC", col = "yellow",
main = "Guinea Pigs' Tooth Growth",
xlab = "Vitamin C dose mg",
ylab = "tooth length", ylim = c(-1, 40), yaxs = "i")
boxplot(len ~ dose, data = ToothGrowth, add = TRUE,
boxwex = 0.25, at = 1:3 + 0.2,
subset = supp == "OJ", col = "orange")
legend("bottomright", bty="n",c("Ascorbic acid", "Orange juice"),
fill = c("yellow", "orange"))
beanplot(len ~ reorder(supp, len, mean) * dose, ToothGrowth,
side = "b", col = list("yellow", "orange"), border = c("yellow2",
"darkorange"), main = "Guinea Pigs' Tooth Growth",
xlab = "Vitamin C dose mg", ylab = "tooth length", ylim = c(-1,
40), yaxs = "i")
legend("bottomright", bty="n",c("Ascorbic acid", "Orange juice"),
fill = c("yellow", "orange"))
#Example with multiple vectors and/or formulas
par(mfrow = c(2,1))
beanplot(list(all = ToothGrowth$len), len ~ supp, ToothGrowth, len ~ dose)
title("Tooth growth length (beanplot)")
#Trick using internal functions to do this with other functions:
mboxplot <- function(...){
graphics::boxplot(beanplot:::getgroupsfromarguments(), ...)
}
mstripchart <- function(..., method = "overplot", jitter = 0.1, offset = 1/3,
vertical = TRUE, group.names, add = FALSE,
at = NULL, xlim = NULL, ylim = NULL,
ylab = NULL, xlab=NULL, dlab = "", glab = "",
log = "", pch = 0, col = par("fg"), cex = par("cex"),
axes = TRUE, frame.plot = axes) {
graphics::stripchart(beanplot:::getgroupsfromarguments(),
method, jitter, offset, vertical, group.names, add,
at, xlim, ylim, ylab, xlab, dlab, glab, log, pch, col, cex,
axes, frame.plot)
}
mstripchart(list(all = ToothGrowth$len), len ~ supp, ToothGrowth, len ~ dose,
xlim = c(0.5,6.5))
title("Tooth growth length (stripchart)")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.