R-ggplot2-箱圖系列（1） basic - 簡書 (jianshu.com)
R-ggplot2-箱圖系列（2） 注釋P值與組間比較 - 簡書 (jianshu.com)
https://www.leansigmacorporation.com/box-plot-with-minitab/
0、加載包和示例數據
library(ggplot2)
library(patchwork)
#組別名最好是字符型；如果是數值類型，最好轉為因子化
ToothGrowth$dose = factor(ToothGrowth$dose)
summary(ToothGrowth)
# len supp dose
# Min. : 4.20 OJ:30 0.5:20
# 1st Qu.:13.07 VC:30 1 :20
# Median :19.25 2 :20
# Mean :18.81
# 3rd Qu.:25.27
# Max. :33.90
head(ToothGrowth)
# len supp dose
# 1 4.2 VC 0.5
# 2 11.5 VC 0.5
# 3 7.3 VC 0.5
# 4 5.8 VC 0.5
# 5 6.4 VC 0.5
# 6 10.0 VC 0.5
Each animal received one of three dose levels of vitamin C (0.5, 1, and 2 mg/day) by one of two delivery methods, orange juice or ascorbic acid
1、基礎繪圖
# x 指定組名列；
# y 指定值的列
p1 = ggplot(ToothGrowth, aes(x=dose, y=len, fill=dose)) +
geom_boxplot()
p2 = ggplot(ToothGrowth, aes(x=dose, y=len, fill=supp)) +
geom_boxplot()
p1 | p2
2、離群點相關
?geom_boxplot
p1 = ggplot(ToothGrowth, aes(x=dose, y=len)) +
geom_boxplot(outlier.color = "red",
outlier.size = 0.5)
p2 = ggplot(ToothGrowth, aes(x=dose, y=len)) +
geom_boxplot(outlier.alpha = 0) #透明度為0，相當于不繪制離群點
p1 + p2
3、隨機抖動點
ggplot(ToothGrowth, aes(x=dose, y=len)) +
geom_boxplot(outlier.color = "red") +
geom_jitter(color = "black", size=0.8)
4、給箱圖添加whisker須線
ggplot(ToothGrowth, aes(x=dose, y=len)) +
stat_boxplot(geom = "errorbar", width = 0.2) +
geom_boxplot()
5、小提琴圖
p1 = ggplot(ToothGrowth, aes(x=dose, y=len, fill=dose)) +
geom_violin()
p2 = ggplot(ToothGrowth, aes(x=dose, y=len, fill=dose)) +
geom_violin() +
geom_boxplot(width=0.2, color="black", alpha=0.5)
p1 | p2
6、排序
- 按照每組的均值等依次降序/升序排列箱圖的順序
參考之前筆記：http://www.lxweimin.com/p/552bee6119ad -
scale_x_discrete()自定義排列順序
ggplot(ToothGrowth, aes(x=dose, y=len)) +
geom_boxplot() +
scale_x_discrete(limits=c("1","2","0.5"))
7、一鍵生成復雜箱圖
library(ggstatsplot)
ggbetweenstats(
data = ToothGrowth,
x = dose,
y = len
)
在下一小節，會學習如何使用
ggpubr包對箱圖組間比較注釋p值結果。
