本文首发于微信公众号
“ 生信补给站 ”
,期待您的关注!!!
“一张统计图形就是从数据到几何对象(geometric object, 缩写为geom, 包括点、线、条形等)的图形属性(aesthetic attributes, 缩写为aes, 包括颜色、形状、大小等)的一个映射。此外, 图形中还可能包含数据的统计变换(statistical transformation, 缩写为stats), 最后绘制在某个特定的坐标系(coordinate system, 缩写为coord)中, 而分面(facet, 指将绘图窗口划分为若干个子窗口)则可以用来生成数据中不同子集的图形。”
—– Hadley Wickham
一 ggplot2 背景介绍
ggplot2
是由
Hadley Wickham
创建的一个十分强大的可视化R包。按照
ggplot2
的绘图理念,Plot(图)= data(数据集)+ Aesthetics(美学映射)+ Geometry(几何对象)。本文将从ggplot2的八大基本要素逐步介绍这个强大的R可视化包。
- 数据(Data)和映射(Mapping)
- 几何对象(Geometric)
- 标尺(Scale)
- 统计变换(Statistics)
- 坐标系统(Coordinante)
- 图层(Layer)
- 分面(Facet)
- 主题(Theme)
二 数据(data) 和 映射(Mapping)
数据:
用于绘制图形的数据,本文主要使用经典的mtcars数据集和diamonds数据集子集为例来画图。
#install.packages("ggplot2")
library(ggplot2)
data(diamonds)
set.seed(1234)
diamond <- diamonds[sample(nrow(diamonds), 2000), ]
head(diamond)
映射:aes()
函数是ggplot2中的映射函数, 所谓的映射即为数据集中的数据关联到相应的图形属性过程中一种对应关系, 图形的颜色,形状,分组等都可以通过通过数据集中的变量映射。
#使用diamonds的数据子集作为绘图数据,克拉(carat)数为X轴变量,价格(price)为Y轴变量。
p <- ggplot(data = diamond, mapping = aes(x = carat, y = price))
#将钻石的颜色(color)映射
颜色
属性:
p <- ggplot(data=diamond, mapping=aes(x=carat, y=price, shape=cut, colour=color))p+geom_point() #绘制点图
#将钻石的切工(cut)映射到
形状
属性:
p <- ggplot(data=diamond, mapping=aes(x=carat, y=price, shape=cut))p+geom_point() #绘制点图
#将钻石的切工(cut)映射到
分组
属性:
#默认分组设置, 即group=1p + geom_boxplot()
#分组(group)也是ggplot2种映射关系的一种, 如果需要把观测点按额外的离散变量进行分组处理, 必须修改默认的分组设置。
p1 <- ggplot(data=diamond, mapping=aes(x=carat, y=price, group=factor(cut)))
p1 + geom_boxplot()
注意:不同的几何对象,要求的属性会有些不同,这些属性也可以在几何对象映射时提供,以下语法与上面的aes中是一样的。
#结果与上述一致
ggplot(data = diamond)+geom_point(aes(x=carat, y=price, colour=color))
ggplot(data = diamond) +geom_point(aes(x=carat, y=price, shape=cut))
ggplot(data = diamond) +geom_boxplot(aes(x=carat, y=price, group=factor(cut)))
ggplot(data = diamond) +geom_point(aes(x=carat, y=price, colour=color,shape=cut))
注:ggplot2支持图层,可以把不同的图层中共用的映射提供给ggplot函数,而某一几何对象才需要的映射参数提供给geom_xxx函数。
三 几何对象(Geometric)
几何对象
代表我们在图中实际看到的图形元素,如点、线、多边形等。数据与映射部分介绍了ggplot函数执行各种属性映射,只需要添加不同的几何对象图层,即可绘制出相应的图形。
此处介绍几种常用的几何对象,geom_histogram用于直方图,geom_bar用于画柱状图,geom_boxplot用于画箱式图等。
直方图
单变量连续变量
:可绘制直方图展示,提供一个连续变量,画出数据的分布。
#以价格(price)变量为例,且按照不同的切工填充颜色
ggplot(diamond)+geom_histogram(aes(x=price, fill=cut))
#设置position=”dodge”,side-by-side地画直方图
ggplot(diamond)+geom_histogram(aes(x=price, fill=cut), position="dodge")
#设置使用position=”fill”,按相对比例画直方图
ggplot(diamond)+geom_histogram(aes(x=price, fill=cut), position="fill")
柱状图
单变量分类变量
:可使用柱状图展示,提供一个x分类变量,画出数据的分布。
#以透明度(clarity)变量为例,且按照不同的切工填充颜色,柱子的高度即为此分类下的数目。
ggplot(diamond)+geom_bar(aes(x=clarity, fill=cut))
注:ggplot2会通过x变量自动计算各个分类的数目。
#直接指定个数,需要通过stat参数,指定geom_bar按特定高度画图
ggplot()+geom_bar(aes(x=c(LETTERS[1:5]),y=1:5), stat="identity")
区分与联系:
直方图把连续型的数据按照一个个等长的分区(bin)切分,然后计数画柱形图。
柱状图是把分类数据,按类别计数。