基本简介
小提琴图 (Violin Plot) 用于显示数据分布及其概率密度，该图类是箱型图和核密度图的结合，主要用来显示数据的分布状况。中间的黑色粗条表示四分位数范围，从其延伸的细黑线代表数据范围，两端为最大值和最小值，而白点则为中位数。
基本用法
vioplot包：vioplot()函数
ggplot2包：plot()函数
示例
vioplot包
vioplot(
x, #数据矩阵
...,
data = NULL,
range = 1.5,
h = NULL,
xlim = NULL,
ylim = NULL,
names = NULL,
horizontal = FALSE,
col = "grey50", #颜色
border = par()$fg, #边界
lty = 1,
lwd = 1
)
# 创建一个数据集
treatment <- c(rep("A", 40) , rep("B", 40) , rep("C", 40) )
value <- c( sample(2:5, 40 , replace=T) , sample(c(1:5,12:17), 40 , replace=T), sample(1:7, 40 , replace=T) )
data <- data.frame(treatment,value)
#绘图
with(data , vioplot(
value[treatment=="A"] , value[treatment=="B"], value[treatment=="C"],
col=rgb(0.1,0.4,0.7,0.7) , names=c("A","B","C")
))
ggplot2包绘制
#加载包
{library(ggplot2)
library(dplyr)}
#创建数据集
data <- data.frame(
name=c( rep("A",500), rep("B",500), rep("B",500), rep("C",20), rep('D', 100) ),
value=c( rnorm(500, 10, 5), rnorm(500, 13, 1), rnorm(500, 18, 1), rnorm(20, 25, 4), rnorm(100, 12, 1) )
)
#基本绘图
ggplot(data, aes(x=name, y=value, fill=name)) + # fill=name allow to automatically dedicate a color for each group
geom_violin()+
theme_bw()
#mpg数据集绘制
mpg$class = with(mpg, reorder(class, hwy, median))
mpg %>%
ggplot( aes(x=class, y=hwy, fill=class)) +
geom_violin() +
xlab("class") +
theme(legend.position="none") +
xlab("")+
theme_bw()
参考文献
[1] https://www.r-graph-gallery.com/violin.html
[2] https://www.jianshu.com/p/e8e69c561f4f