[会员专享] nature communications图表复现之个性化地图绘制

R语言数据分析指南

2023-03-11

欢迎关注R语言数据分析指南

❝
本节来复现「nature communications」上一篇文章中地图的绘制方法，下面来进行复现过程，由于未找到作者提供的数据信息，小编自己构建了绘图文件因此结果与原图有出入仅供参考。整张图均使用R代码进行绘制，「数据+代码已经上传2023VIP群，加群的观众老爷请自行下载」
❞

2023-VIP购买方式

❝
有需要购买2023年文档的朋友欢迎到小编的「淘宝店铺」 「R语言数据分析指南」下单购买初始价格99元，微信发小编订单截图即可邀请进会员交流群，小编的文档为按年售卖，只包含当年度的「除系列课程外」的文档，有需要往年文档的朋友也可下单购买，需要了解更多信息的朋友欢迎交流咨询，「添加小编微信请备注，以提高沟通效率」
❞

论文原图

Land-use diversity predicts regional bird taxonomic and functional richness worldwide

复现图

加载R包

# install.packages("rnaturalearthdata")

# install.packages("tidygeocoder")
library(tidygeocoder)
library(tidyverse)
library(sf)
library(camcorder)
library(scico)
library(rnaturalearth)
library(terra)
library(tidyterra)
library(geodata)
library(cowplot)
library(ggsci)

绘制点图


df1 <- read_tsv("data.xls") %>% filter(type=="Taxonomic richness") %>% 
  select(5:9) %>% 
  group_by(REALM) %>% 
  slice_head(n=1)

df1$REALM <- factor(df$REALM,levels = c("Nearctic","Palearctic","Indomalayan","Neotropic","Afrotropic","Australasia"))  

plot1 <- df1 %>% ggplot(aes(y = fct_rev(REALM))) +
  theme_bw()+
  geom_errorbarh(aes(xmin=Lower_ci,xmax=Upper_ci),height=0.1) +
  geom_point(aes(x=visregFit,color=REALM),fill="black",size=3,show.legend = F) +
  labs(x="Taxonomic richness",y=NULL)+
  scale_color_npg()+
  theme(axis.ticks.y= element_blank(),
        axis.title.y= element_blank(),
        axis.title.x = element_text(color="black",size=8,face="bold"),
        axis.text.y=element_text(color="black",size=8,face="bold"),
        axis.text.x=element_text(color="black",size=8,face="bold"))


df2 <- read_tsv("data.xls") %>%
  filter(type=="Functional richness") %>% 
  select(5:9) %>% 
  group_by(REALM) %>% 
  slice_head(n=1)

df2$REALM <- factor(df$REALM,levels = c("Nearctic","Palearctic","Indomalayan","Neotropic","Afrotropic","Australasia"))  

plot2 <- df2 %>% ggplot(aes(y = fct_rev(REALM))) +
  theme_bw()+
  geom_errorbarh(aes(xmin=Lower_ci,xmax=Upper_ci),height=0.1) +
  geom_point(aes(x=visregFit,color=REALM),fill="black",size=3,show.legend = F) +
  labs(x="Functional richness",y=NULL)+
  scale_color_npg()+
  theme(axis.ticks.y= element_blank(),
        axis.title.y= element_blank(),
        axis.title.x = element_text(color="black",size=8,face="bold"),
        axis.text.y=element_text(color="black",size=8,face="bold"),
        axis.text.x=element_text(color="black",size=8,face="bold"))