跟着Nature Metabolism学作图:R语言ggplot2各种各样柱形图(1)

本文涉及的产品
简介: 跟着Nature Metabolism学作图:R语言ggplot2各种各样柱形图(1)

论文

Single-cell profiling of vascular endothelial cells reveals progressive organ-specific vulnerabilities during obesity

https://www.nature.com/articles/s42255-022-00674-x#Sec58

s42255-022-00674-x.pdf

https://github.com/Osynchronika/sc_EC_obesity_atlas

大部分 作图的数据都有,可以试着用论文中提供的数据复现一下论文中的图

论文中figure2和figure3中有很多种柱形图,争取把每个种类的柱形图都复现一下

加载作图用到的R包

library(readxl)
library(tidyverse)
library(ggplot2)

首先是最普通的柱形图

figure3m

image.png

示例数据集如下

image.png

作图代码

fig3m.df<-read_excel("data/20230207/ggplot2barplot.xlsx",
                     sheet = "fig3m")
fig3m.df
ggplot(data=fig3m.df,aes(x=x,y=y))+
  geom_col(fill="#5ab033",color="black")+
  theme_classic()+
  scale_y_continuous(expand = expansion(mult = c(0,0)),
                     limits = c(0,0.5))+
  labs(x=NULL,y="log(FC)")+
  theme(panel.grid.major.y = element_line(),
        axis.text.x = element_text(angle=90,face="italic",
                                   vjust=0.5,hjust=1))

更多的时候会对数值进行排序,从小到大,或者从大到小

fig3m.df %>% 
  arrange(y) %>% 
  mutate(x=factor(x,levels = x)) %>% 
  ggplot(aes(x=x,y=y))+
  geom_col(fill="#5ab033",color="black")+
  theme_classic()+
  scale_y_continuous(expand = expansion(mult = c(0,0)),
                     limits = c(0,0.5))+
  labs(x=NULL,y="log(FC)")+
  theme(panel.grid.major.y = element_line(),
        axis.text.x = element_text(angle=90,face="italic",
                                   vjust=0.5,hjust=1)) -> p1

fig3m.df %>% 
  arrange(desc(y)) %>% 
  mutate(x=factor(x,levels = x)) %>% 
  ggplot(aes(x=x,y=y))+
  geom_col(fill="#5ab033",color="black")+
  theme_classic()+
  scale_y_continuous(expand = expansion(mult = c(0,0)),
                     limits = c(0,0.5))+
  labs(x=NULL,y="log(FC)")+
  theme(panel.grid.major.y = element_line(),
        axis.text.x = element_text(angle=90,face="italic",
                                   vjust=0.5,hjust=1)) -> p2

p1/p2

image.png

柱形图除了水平摆放,也可以垂直摆放,我们把作图代码里的x和y对调位置就行,如果数据集里的数据有正有负,那么柱子呈现的就是既有朝上的,又有朝下的

比如这个figure r s

image.png

代码


fig3r.df<-read_excel("data/20230207/ggplot2barplot.xlsx",
                     sheet = "fig3r")

fig3r.df %>% 
  mutate(x=factor(x,levels = c("Vcam1","Pecam1","Alcam","Icam1","Gja4","Gja5","F11r"))) %>% 
  ggplot(aes(x,y))+
  geom_col(fill="#ee7770",color="black")+
  geom_text(aes(y=c(-0.01,-0.01,-0.01,-0.01,-0.01,0.01,0.01),
                label=x),
            angle=90,hjust=c(1,1,1,1,1,0,0),
            color=c("#bf1818","#bf1818","#bf1818","#bf1818","black","black","blue"))+
  theme_classic()+
  theme(axis.line.x = element_blank(),
        axis.text.x = element_blank(),
        axis.ticks.x = element_blank(),
        panel.grid.major.y = element_line(),
        plot.title = element_text(hjust = 0.5))+
  scale_y_continuous(limits = c(-0.25,0.25),
                     breaks = c(-0.25,seq(-0.2,0.2,by=0.1),0.25),
                     expand = expansion(mult = c(0,0)),
                     labels = c("",seq(-0.2,0.2,by=0.1),""))+
  labs(x=NULL,y="log(FC)",title = "Art")

image.png

今天的推文就介绍这么多,明天继续

示例数据和代码可以给推文点赞,然后点击在看,最后留言获取

欢迎大家关注我的公众号

小明的数据分析笔记本

小明的数据分析笔记本 公众号 主要分享:1、R语言和python做数据分析和数据可视化的简单小例子;2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记;3、生物信息学入门学习资料及自己的学习笔记!

微信公众号好像又有改动,如果没有将这个公众号设为星标的话,会经常错过公众号的推文,个人建议将 小明的数据分析笔记本 公众号添加星标,添加方法是

点开公众号的页面,右上角有三个点

image.png

点击三个点,会跳出界面

image.png

直接点击 设为星标 就可以了

image.png

相关实践学习
基于函数计算一键部署掌上游戏机
本场景介绍如何使用阿里云计算服务命令快速搭建一个掌上游戏机。
建立 Serverless 思维
本课程包括: Serverless 应用引擎的概念, 为开发者带来的实际价值, 以及让您了解常见的 Serverless 架构模式
相关文章
|
7月前
|
存储 Go
速绘丨GO富集气泡图绘制方法,利用R语言ggplot2包快速绘制,完整脚本可重复绘图
速绘丨GO富集气泡图绘制方法,利用R语言ggplot2包快速绘制,完整脚本可重复绘图
|
7月前
|
数据采集 机器学习/深度学习 SQL
绝不可错过!R语言与ggplot2实现SCI论文数据分析神器
绝不可错过!R语言与ggplot2实现SCI论文数据分析神器
115 0
|
10月前
|
数据可视化 数据挖掘 定位技术
跟着Nature Communications学作图:R语言ggplot2画世界地图并用md语法添加文字标签
跟着Nature Communications学作图:R语言ggplot2画世界地图并用md语法添加文字标签
|
10月前
|
数据可视化 数据挖掘 Python
跟着Oncogene学作图:R语言gggenomes画桑基图
跟着Oncogene学作图:R语言gggenomes画桑基图
|
10月前
|
数据可视化 数据挖掘 Python
跟着NatureCommunications学作图:R语言ggtree根据分组给进化树上色
跟着NatureCommunications学作图:R语言ggtree根据分组给进化树上色
|
10月前
|
数据可视化 数据挖掘 Python
跟着NatureMetabolism学作图:R语言ggplot2转录组差异表达火山图
跟着NatureMetabolism学作图:R语言ggplot2转录组差异表达火山图

热门文章

最新文章