ggplot2-绘图-箱线图

前面的啰嗦

本教程对绘图所需数据和绘图过程进行了一步步的解释,稍显啰嗦,没时间细看或不想一步步看的朋友可以直接点击目录最后一节:
箱线图完整版代码

绘图数据说明

用到的数据是mpg数据集,这是ggplot2包自带的数据集,可以使用以下代码查看该数据集的基本情况:

library(ggplot2)
?mpg # 或者help(mpg)

从帮助文档中可以看出这是一个234x11的数据框,它记录了美国1999年和2008年部分汽车的制造厂商,型号,类别,驱动程序和耗油量等信息。
下面看一下这11列的情况:

library(ggplot2)
# str查看数据集的基本情况
str(mpg)
## tibble [234 x 11] (S3: tbl_df/tbl/data.frame)
##  $ manufacturer: chr [1:234] "audi" "audi" "audi" "audi" ...
##  $ model       : chr [1:234] "a4" "a4" "a4" "a4" ...
##  $ displ       : num [1:234] 1.8 1.8 2 2 2.8 2.8 3.1 1.8 1.8 2 ...
##  $ year        : int [1:234] 1999 1999 2008 2008 1999 1999 2008 1999 1999 2008 ...
##  $ cyl         : int [1:234] 4 4 4 4 6 6 6 4 4 4 ...
##  $ trans       : chr [1:234] "auto(l5)" "manual(m5)" "manual(m6)" "auto(av)" ...
##  $ drv         : chr [1:234] "f" "f" "f" "f" ...
##  $ cty         : int [1:234] 18 21 20 21 16 18 18 18 16 20 ...
##  $ hwy         : int [1:234] 29 29 31 30 26 26 27 26 25 28 ...
##  $ fl          : chr [1:234] "p" "p" "p" "p" ...
##  $ class       : chr [1:234] "compact" "compact" "compact" "compact" ...

每一列的解释如下:

  • model: model name
  • displ: 发动机排量(engine displacement, in litres)
  • year: 发行时间,仅有1999和2008
  • cyl: 汽缸数目(number of cylinders)
  • trans: 变速器类型(type of transmission)
  • drv: 表示驱动程序,是一个分类变量,详细解释为:the type of drive
    train, where f = front-wheel drive, r = rear wheel drive, 4 = 4wd
  • cty: city miles per gallon,每加仑城市里程,耗油量的一种计量方式。
  • hwy: highway miles per
    gallon,每加仑公路行驶里程,也是耗油量的一种计量方式。

基本箱线图

下面对hwy绘制箱线图,查看这些汽车的耗油量分布情况。

library(ggplot2)
ggplot(data=mpg, aes(x="", y=hwy)) +
  geom_boxplot()
image

可以看到耗油量的中位数在25附近,其中有两个离群点,这两辆汽车的耗油量很大。

这里绘制的图是一个横向的箱线图,如果要绘制纵向的箱线图(即上图顺时针旋转90度),只需要改一下代码中aes的参数,设置为x=hwy即可,如下所示:

ggplot(data=mpg, aes(x=hwy, y="")) +
  geom_boxplot()
image

分组箱线图

下面绘制一下分组箱线图,看一下不同驱动程序的汽车的耗油量分布,其中驱动程序是mpg数据集中的drv变量,这是一个分类变量

class(mpg$drv)
## [1] "character"
table(mpg$drv)
## 
##   4   f   r 
## 103 106  25
ggplot(data=mpg, aes(x=factor(drv), y=hwy, fill=factor(drv))) +
  geom_boxplot()
image

由于mpg$drv的类别是character,而ggplot绘图进行分组所需要的变量类型是factor,因此需要将drv变量转为factor类别,但这里存在一个问题,就是横坐标轴标题出现了factor字样,图例的标题中也出现了factor字样,不是很美观,这个问题的解决方法有两个:

  1. 直接对数据框的drv变量进行修改,将其改为factor类别:
mpg_new <- mpg
mpg_new$drv <- as.factor(mpg_new$drv)
ggplot(data=mpg_new, aes(x=drv, y=hwy, fill=drv)) +
  geom_boxplot()
image
  1. 修改坐标轴与图例的标题,将其修改为drv:
ggplot(data=mpg, aes(x=factor(drv), y=hwy, fill=factor(drv))) +
  geom_boxplot() + 
  scale_fill_discrete(name='drv') + 
  labs(x="drv") # x表示横坐标轴的标题
image
  # scale_一系列的函数都是用来配置图例的,scale_fill_discrete指的是修改fill的图例,和aes中的fill相对应。

美化分组箱线图

上面的箱线图存在一些地方不是很好看,主要有以下几点:

  1. 图例的颜色与横坐标轴的三个标签信息是一致的,存在冗余,因此可以去掉图例,去掉图例的代码为guides(fill=FALSE)
  2. 没有添加标题,可以在labs函数中添加标题
  3. 所用主题可以优化,可以使用theme_系列函数

下面使用mpg_new数据框进行操作,避免不停的对drv列进行操作。

mpg_new <- mpg
mpg_new$drv <- as.factor(mpg_new$drv)
ggplot(data=mpg_new, aes(x=drv, y=hwy, fill=drv)) +
  geom_boxplot() +
  guides(fill=FALSE) + # fill指的是图例中的fill,和aes中对应
  labs(title="Boxplot of highway miles per gallon") +
  theme_minimal() +
  theme(plot.title=element_text(hjust=0.5)) # 这一句让标题居中,ggplot2默认的标题是左对齐的
## Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> =
## "none")` instead.
image

这里有一个warning,表示guides(<scale> = FALSE)这种用法已经废弃,这种用法是旧版本ggplot2的用法,新版本的应使用guides(<scale> = "none")(我的ggplot2版本是3.3.5),修改后的代码如下:

ggplot(data=mpg_new, aes(x=drv, y=hwy, fill=drv)) +
  geom_boxplot() +
  guides(fill="none") + # fill指的是图例中的fill,和aes中对应
  labs(title="Boxplot of highway miles per gallon") +
  theme_minimal() +
  theme(plot.title=element_text(hjust=0.5))  #这一句让标题居中,ggplot2默认的标题是左对齐的
image

这里有一点需要注意的是theme_minimaltheme两个函数的顺序,由于theme_minimal是一种简化式的主题设置,ggplot2遵循图层叠加的原理,如果将theme_minimal放在theme后面,那么图标标题仍然是左对齐,而不是居中显示。

箱线图完整版代码

基本箱线图的完整代码

  1. 横向箱线图
library(ggplot2)
ggplot(data=mpg, aes(x="", y=hwy)) +
  geom_boxplot() +
  labs(x="") # 去除横坐标轴的标签
image
  1. 纵向箱线图
library(ggplot2)
ggplot(data=mpg, aes(x=hwy, y="")) +
  geom_boxplot() +
  labs(y="") # 去除纵坐标轴的标签
image

分组箱线图的完整代码

# 首先将数据框中的drv修改为factor类型
library(ggplot2)
mpg_new <- mpg
mpg_new$drv <- as.factor(mpg_new$drv)

# 其次绘图
ggplot(data=mpg_new, aes(x=drv, y=hwy, fill=drv)) +
  geom_boxplot() +
  guides(fill="none") + # fill指的是图例中的fill,和aes中对应,fill=none表示不显示图例
  labs(title="Boxplot of highway miles per gallon") + # title是图表标题
  theme_minimal() +
  theme(plot.title=element_text(hjust=0.5))  #这一句让图表标题居中,ggplot2默认的标题是左对齐的
image
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,636评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,890评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,680评论 0 330
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,766评论 1 271
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,665评论 5 359
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,045评论 1 276
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,515评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,182评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,334评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,274评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,319评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,002评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,599评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,675评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,917评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,309评论 2 345
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,885评论 2 341

推荐阅读更多精彩内容