一、beeswarm作为一维散点图包R包,可以生成点不重复的图,与stripchart的区别就是等值点不会重叠到一起,下图展示了stripchart与beeswarm图的区别:

beeswarm-蜜蜂图-LMLPHP

stripchart(decrease ~ treatment, data = OrchardSprays, vertical = TRUE, log = "y", method = 'jitter', jitter = 0.2, cex = 1,pch = 16, col = rainbow(8),main = 'stripchart')

beeswarm(decrease ~ treatment, data = OrchardSprays, log = TRUE, pch = 16, col = rainbow(8),main = 'beeswarm')

  OrchardSprays是R自带的数据,decrease 和 treatment必须是数据的header

  pch:描点的样式

  jitter:把横坐标抖散(jitter),使本来都拥有相同坐标的点的横坐标稍有不同。jitter是基类函数{base},无需调用任何包

  

二、不同的beeswarm图:

  1.The color of individual points can be specified

install.packages(“beeswarm”)
library(beeswarm)
data(breast)#必须先library(beeswarm),才能导入breast数据集
head(breast)
getwd()
setwd('/home/name/R')#自己设置路径
png('diff_color.png',width=400*2,height=400*2,res=72*2,type='cairo-png')
beeswarm(time_survival ~ ER, data=breast,pch=16, pwcol=1+as.numeric(event_survival),xlab="",ylab="Follow-up time (months)",labels = c("ER neg", "ER pos"))
legend("topright",legend=c("yes","no"),title="Censored",pch=16,col=1:2)
dev.off()

beeswarm-蜜蜂图-LMLPHP

  2.Compare the four methods for arranging points#点的组织方式可以多种,如("swarm", "center", "hex", "square")

## Generate some random data
set.seed(123)
distro <- list(runif = runif(100, min = -3, max = 3), rnorm = rnorm(100)) png("arrange_points.png",width=600*2,height=600*2,res=72*2,type='cairo-png')
par(mfrow=c(2,2))#两行两列画布
for (m in c("swarm", "center", "hex", "square")) {
beeswarm(distro, col = 2:3, pch = 16,method = m,
main = paste('method = "', m, '"', sep = ''))
} dev.off()

beeswarm-蜜蜂图-LMLPHP

  3.Combining beeswarms with "boxplot" or "bxplot"

  

png("box_plus_beeswarm.png",width=600*2,height=600*2,res=72*2,type='cairo-png')
par(mfrow=c(1,2))
boxplot(len ~ dose, data = ToothGrowth,outline = FALSE,main = 'boxplot + beeswarm')
#outline = FALSE :avoid double-plotting outliers, if any
beeswarm(len ~ dose, data = ToothGrowth,col = 4, pch = 16, add = TRUE) beeswarm(len ~ dose, data = ToothGrowth,col = 4, pch = 16,main = 'beeswarm + bxplot')
bxplot(len ~ dose, data = ToothGrowth, add = TRUE)
dev.off()

beeswarm-蜜蜂图-LMLPHP

  4.New arguments "side" and "priority"

png("side_priority.png",width=600*2,height=600*2,res=72*2,type='cairo-png')
par(mfrow=c(2,3))
beeswarm(distro, col = 2:4, main = 'Default') beeswarm(distro, col = 2:4, side = -1, main = 'side = -1') beeswarm(distro, col = 2:4, side = 1, main = 'side = 1') beeswarm(distro, col = 2:4, priority = "descending", main = 'priority = "descending"') beeswarm(distro, col = 2:4, priority = "random", main = 'priority = "random"') beeswarm(distro, col = 2:4, priority = "density", main = 'priority = "density"') dev.off()

beeswarm-蜜蜂图-LMLPHP

  参考:http://www.cbs.dtu.dk/~eklund/beeswarm/

  

05-04 03:47