生信小白7行代码画一个高颜值的世界地图

不要害怕编程，它其实像泡方便面一样简单。先在百度搜索R，然后像安装QQ一样，一直选择下一步，安装成功后，双击打开桌面R的快捷方式。

复制粘贴下面代码，敲一个回车键，就能得到下面的世界地图。如果这个图上面，恰巧能表示你选取的200份小麦全球的分布情况，岂不是心里很美。

代码

install.packages("maps") ## 安装R包
install.packages("ggplot2")
library("ggplot2") ### 加载ggplot2
library("maps") ### 加载地图包
world_map <- map_data("world") ### 导入地图数据 
head(world_map) ### 看看地图数据是啥样子
world_map$size <- 1:dim(world_map)[1] ### 生产一些随机的数据

### 画图开始 ####
ggplot(world_map, aes(x = long, y = lat, group = group, fill = size)) +
  geom_polygon(colour = "white",size=0.2) + scale_fill_distiller(palette = "YlGnBu",direction = 1)

效果图

地图数据与画图参数讲解

world_map初始样子就是下面这个样子，有6列数据，分别是经度，维度，组别，顺序，国家，地区

现在我们需要按照一些特性，给不同国家画上不同颜色，使地图更有意义。比如每个国家博士数量，GDP，国足进世界杯次数，SCI文章数量，搜集的小麦材料各个国家占比等等。

这里我们使用一个随机的数据。world_map$size <- 1:dim(world_map)[1]数值是从1到n的数据，没啥具体的含义，就为了一个例子。

我们画自己的图就只要修改这个数据，怎么修改呢？可以使用write.csv把地图数据写出来（这是个exel格式）然后用read.csv再读进去就可以了。

网上很多类似教程，其他参数不再细讲，聪明如你，一定可以。

猜你喜欢

系列教程：微生物组入门 Biostar 微生物组  宏基因组

专业技能：学术图表 高分文章 生信宝典 不可或缺的人

一文读懂：宏基因组 寄生虫益处 进化树

必备技能：提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

扩增子分析：图表解读 分析流程 统计绘图

在线工具：16S预测培养基 生信绘图

科研经验：云笔记  云协作 公众号

编程模板: Shell  R Perl

生物科普:  肠道细菌 人体上的生命 生命大跃进  细胞暗战 人体奥秘  

写在后面

为鼓励读者交流、快速解决科研困难，我们建立了“宏基因组”专业讨论群，目前己有国内外5000+ 一线科研人员加入。参与讨论，获得专业解答，欢迎分享此文至朋友圈，并扫码加主编好友带你入群，务必备注“姓名-单位-研究方向-职称/年级”。PI请明示身份，另有海内外微生物相关PI群供大佬合作交流。技术问题寻求帮助，首先阅读《如何优雅的提问》学习解决问题思路，仍未解决群内讨论，问题不私聊，帮助同行。

学习16S扩增子、宏基因组科研思路和分析实战，关注“宏基因组”

点击阅读原文，跳转最新文章目录阅读