这一讲,我们将介绍如何用ggplot2软件包中的函数geom_boxplot()绘制复杂箱式图。1. 概论
函数geom_boxplot()简化格式为:
geom_boxplot(outlier.colour="black", outlier.shape=16,outlier.size=2, notch=FALSE)
·Outlier.colour, outlier.shape, outlier.size:分别表示离群点的颜色,形状和大小。
·notch:逻辑值。如果为TRUE,则绘制一个缺口箱式图。缺口在中位数及其置信区间,该置信区间通常基于中位数+/- 1.58 * IQR / sqrt(n)。缺口用于比较组;如果两个箱式的缺口不重叠,则有力证明中位数不同。
使用R内置的ToothGrowth数据集:
# 将可变剂量从数字转换为因子变量ToothGrowth$doseas.factor(ToothGrowth$dose)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
确保使用上述R脚本将变量剂量转换为因子变量。
3. 基本箱式图library(ggplot2)# 基本箱式图,将基本箱式图保存在p中p<-ggplot(ToothGrowth, aes(x=dose, y=len))+geom_boxplot() #表示对ggplot()函数中绘图内容进行geom_boxplot()函数运算p #打印基本箱式图,即输入p# 转置箱式图p+coord_flip()# 绘制缺口箱式图ggplot(ToothGrowth, aes(x=dose, y=len))+geom_boxplot(notch=TRUE)# 更改离群值的颜色、形状和大小ggplot(ToothGrowth, aes(x=dose, y=len))+geom_boxplot(outlier.colour="red", outlier.shape=8,outlier.size=4)