不要害怕编程，它其实像泡方便面一样简单。先在百度搜索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再读进去就可以了。

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

猜你喜欢

10000+：菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature Cell专刊 肠道指挥大脑

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

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

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

必备技能：提问 搜索  Endnote

文献阅读 热心肠 SemanticScholar Geenmedical

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

16S功能预测   PICRUSt  FAPROTAX  Bugbase Tax4Fun

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

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

编程模板: Shell  R Perl

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

写在后面

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

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

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