前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >ggplot2多维分面多图层对应规则

ggplot2多维分面多图层对应规则

作者头像
数据小磨坊
发布2018-04-11 17:28:05
1K0
发布2018-04-11 17:28:05
举报
文章被收录于专栏:数据小魔方

今天只给大家讲一个知识点,是属于ggplot2高阶用法中的分面与多图层关系如何对应,这个用法之前困扰我很久,也是最近帮朋友做东西才发现这个漏洞,于是感觉分享给大家。

ggplot2的多维分面系统非常完美,可以让我们非常方便的将一个多维度的复杂图形按照某个维度的类别进行矩阵化,使得单个类别的信息更加清晰明了,数据呈现直观易懂。

通常我们可以通过如下步骤做出一个基于地理位置的分面图来。

library(ggplot2) library(plyr) library(maptools) library(Cairo) library(RColorBrewer) library(dplyr)

CHN_adm2 <- readShapePoly("D:/R/rstudy/CHN_adm/CHN_adm2.shp") CHN_adm2_1 <- fortify(CHN_adm2) data1 <- CHN_adm2@data data2 <- data.frame(id=row.names(data1),data1) china_map_data <- join(CHN_adm2_1,data2, type = "full") dongsansheng <-subset(china_map_data,NAME_1==c("Heilongjiang","Jilin","Liaoning")) dongsansheng$NAME_1<-as.character(dongsansheng$NAME_1)

mydata<-read.csv("D:/R/rstudy/State/dongsansheng.csv",header=T)

dongsansheng<-within(dongsansheng,{ NAME_1[NAME_1=="Heilongjiang"]="黑龙江" NAME_1[NAME_1=="Jilin"]="吉林" NAME_1[NAME_1=="Liaoning"]="辽宁" })

ggplot()+ geom_polygon(data=dongsansheng,aes(x=long,y=lat,group=group),colour="grey40",fill="white") + facet_grid(.~NAME_1,scales = "free")+ coord_map("polyconic")+ theme_void(base_size=18)

然而问题来了,以上图形仅仅基于同一个图层进行维度分面,倘若我有如下需求,不仅要对地图进行分面,而且要在单个区域地图上呈现一些点信息、线条的信息,这就意味着我们需要在保持分面的基础上,叠加图层,那么我们给分面函数指定的分面规则是否能够作用于第二个图层呢,或者说想要让分面函数同事控制所有图层应该 如何进行参数设定呢,分面参数的控制权限到底有多高呢?

带着以下疑问,我们先尝试着在原始分面图层的基础上叠加一个散点图层。

midpos <- function(x) mean(range(x,na.rm=TRUE)) centres <- ddply(dongsansheng,.(NAME_2),colwise(midpos,.(long,lat))) mydata<-mydata%>%merge(centres) ggplot() + geom_polygon(data=dongsansheng,aes(x=long,y=lat,group=group),colour="grey40",fill="white") + geom_point(data=mydata,aes(x=long,y=lat,size=zhibiao),colour="red",alpha=.5)+ facet_grid(.~NAME_1,scales = "free")+ coord_map("polyconic")+ scale_size_area(max_size=8)+ theme_void(base_size=18)

然而遗憾的是,我们得到的结果是这样的,分面函数仅仅控制了第一个图层(也就是地图的图层),却对第二个图层(散点图层没有任何影响),这不是我们想要的结果,我们想要的是这个分面参数同事完成地图和散点图的对应区域分割、匹配。

实际上以上结果并不难解释,因为我们在分面参数设定是,参数设置的依据是第一个图层的数据源中的NAME_1字段,但是我们并没有保证第二个图层中有这个同属性,同名称的字段。

事实上为了进行接下来的案例讲解,我确实在散点图的数据源中设定了一个跟多边形(也就是第一个图层)数据源的NAME_1同属性的字段(里面记录的都是类别相同的省份名称),为了对比效果暂时命名为Province。但是分面函数只能在以上两个数据框中找到第一个图层数据源中含有NAME_1变量,而第二个图层的数据源中尽管有同性质的变量,但是名称不同,分面函数是无法识别的,因为忽略了对图层二的分面操作。

接下来我将图层二中的省份名称变量更改为更图层一中名称相同,再看下结果:

mydata<-mydata%>%rename(NAME_1=Province) ggplot() + geom_polygon(data=dongsansheng,aes(x=long,y=lat,group=group),colour="grey40",fill="white") + geom_point(data=mydata,aes(x=long,y=lat,size=zhibiao),colour="red",alpha=.5)+ facet_grid(.~NAME_1,scales = "free")+ coord_map("polyconic")+ scale_size_area(max_size=8)+ theme_void(base_size=18)

这下效果一目了然,分面函数很顺利的识别了通过设定的分面依据字段NAME_1,在图层一、图层二中均探测到了同名字段,而且字段结构类别均一致,分面操作成功。

这个知识点相对难理解,属于ggplot2高阶用法中比较深奥的部分,与此相同的还有分面参数的控制权限范围问题,即分面参数的控制权限范围到底有多高,是否可以控制ggplot父函数内设定的同名参数,关于这一点儿,我先挖个坑,以后有时间再填,如果你感兴趣,可以自己通过以上提示,使用今天的案例数据修改参数,自己探索,相信在实践中提升技能,你学到的才是硬实力。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2017-08-26,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 数据小魔方 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档