首页/文章/ 详情

分享18种特色科研图表的绘制方法

6月前浏览8793

1.彩色标签条形图


   


条形图,展示分析结果直观明了,绘制方法非常简单,在科研数据的展示中占有非常重要的地位,差不多每篇学术论文都可以看到它。可是,条形图看似简单,但有些条形图的绘制还是有一点难度的,比如下图这样的。


我们可发现上图最大特点(也是绘制难点)是“两级”刻度标签的添加,并且使坐标轴(Y轴)刻度标签的颜色是与“条形”的颜色是完全一致的,这里其实涉及到数字“0”和空值的妙用。关于这类图表的绘制,比如气泡图,箱型图,小提琴图等,其实利用ggplot2包的一个小小的“Bug”就可以完美实现。

#读入数据;
dt1<-read.table("koenrich.xls",sep = "\t",header = T)
#指定纵轴标签顺序,按照输入文件的顺序排序,否则默认按照首字母顺序,同时逆序绘制,保持与表格顺序一致;
dt1$KEGG_A_Class<-factor(dt1$KEGG_A_Class,
                         levels = rev(unique(dt1$KEGG_A_Class)),
                         ordered = TRUE)
dt1$KEGG_B_Class<-factor(dt1$KEGG_B_Class,
                         levels = rev(unique(dt1$KEGG_B_Class)),
                         ordered = TRUE)
# 加载ggplot2包;
library(ggplot2)
#建立数据(Genes.Number, KEGG_B_Class)与图形(点)的映射关系;
p1<-ggplot(dt1, aes(x=Genes_Number,
                    y=KEGG_B_Class,
                    fill=KEGG_A_Class,na.rm = FALSE))+
  geom_bar(stat="identity",na.rm = FALSE)+
 geom_text(aes(x=Genes_Number,y=KEGG_B_Class,label=label),size=2.5,hjust="left",nudge_x=0.1)+
  scale_x_continuous(limits = c(0, 25),expand=expansion(mult = c(0, .1)))+
  labs(x="Number of Gene",y="",title="KEGG pathway anotation")
p1

#获取颜色;
g <- ggplot_build(p1)
mycol<-g$data[[1]]["fill"]
col<-rev(mycol[,1])
#将A Class对应的颜色设为深黑色;
num <- rev(dt1$Genes_Number)
index <- which(num==0)
col[index] <- "grey10"





#自定义图表主题,对图表做精细调整;
top.mar=0.2
right.mar=0.2
bottom.mar=0.2
left.mar=0.2
mytheme1<-theme(plot.title = element_text(size = rel(1),hjust = 0.5,face = "bold"),
                          axis.title = element_text(size = rel(1)),
                          axis.text.y = element_text(size=rel(0.85),
                                                     colour =col,face = "bold"),
                          legend.position = "none",
                          plot.margin=unit(x=c(top.mar,right.mar,
                                               bottom.mar,left.mar),
                                           units="inches"))
#查看绘图效果;
p1+mytheme1




详细教程链接:
《如何绘制好看的条形图?标签颜色相对应的这种。。。》


   

   

     

2.发散PCA散点图


   




PCA (Principal Component Analysis) 散点图,是科研文章中常见的一类散点图。下面就以PCA散点图为例,为大家介绍如何使用ggplot2绘制好看的散点图。

#读入数据;
dt <- read.csv("PCA_data.csv")
#载入ggplot2绘图包;
library(ggplot2)
library(ggh4x)

#绘制实心散点图;
p1 <- ggplot(dt,aes(x=PC1,y=PC2,fill=Diagnosis))+
  stat_centroid(aes(xend = PC1, yend = PC2, colour = Diagnosis),
                geom = "segment", crop_other = F,
                alpha=0.3,size = 1,show.legend = F)+
  geom_point(size=3,alpha=0.7,
             color="white",shape = 21,show.legend = T)+
  scale_color_manual(name="",
                     values = c("#FF9999","#c77cff"))+
  scale_fill_manual(name="",
                     values = c("#FF9999","#c77cff"))+
  scale_x_continuous(expand=expansion(add = c(0.7,0.7)),
                     limits=c(-10,5))+
  scale_y_continuous(expand=expansion(add = c(0.5,0.5)),
                     limits=c(-7.5,5))+
  guides(x = "axis_truncated",y = "axis_truncated")

p1




详细教程链接:
《不就是PCA图么,还能有什么新花样?》


   

   

     

3.漂亮的玫瑰图


   




条形图(或柱状图)是科研文章中比较常见的数据展示方式。除了墨守常规的样式,其实你也可以试试南丁格尔玫瑰图,特别是当你的数据具有周期性的时候,比如四季变化、昼夜节律等,非常适合。

#读入数据;
dt <- read.csv("test_data.csv")
#实用month.abb[]将月份转换为英文缩写;
dt$Month <- month.abb[dt$Month]
#载入相关的R包;
library(dplyr)
library(ggplot2)
#转成tibble格式;
df <- as_tibble(dt)
#将月份转换成因子,固定顺序;
df$Month <- factor(df$Month,levels = unique(df$Month),ordered = T)
#将季度数值转成字符;
df$Quarter <- as.character(df$Quarter)

#提取作图数据;
df1 <- filter(df,Group=="Day")
#自定义颜色;
subcol<-rainbow(12)
mycol1<-colorRampPalette(subcol[1:4])(12)
mycol2<-colorRampPalette(c("#A5CC26","yellow","orange","tomato"))(12)

#绘制相互叠合的玫瑰图;
p1 <- ggplot(df1,aes(x = Month, y=Length,fill=Month))+
  geom_col(width = 1.2,color=NA,alpha=1)+
  geom_text(aes(label = Length),nudge_y=-2,
            colour="white",size=3)+
  ylim(-1.2,17)+
  scale_fill_manual(values = mycol1)+
  coord_polar(start = 0)+
  theme_void()
p1




详细教程链接:
《遥遥领先!如何绘制漂亮的玫瑰图?》


   

   

     

4.富集分析圈图


   




OmicShare的GO、KEGG富集分析工具非常强大,只需上传目的基因文件,即可进行常见物种的GO、KEGG富集分析。分析结果中,直接包含数据可视化结果,比如柱状图、圈图、网络图等。

工具链接:
https://www.omicshare.com/tools/Home/Soft/gogseasenior
https://www.omicshare.com/tools/Home/Soft/pathwaygseasenior

详细教程链接:
《如何绘制GO、KEGG富集圈图》


   

   

     

5.相关性网络热图


   




网络热图,应用于多组学关联分析,可展示组学内(热图)和组学间(网络图)的相关性结果。其中,组学内相关性热图,用于展示每两个元素(环境因子、代谢物、性状等)之间相关系数;组学间相关性网络图,可展示16S-OTU、宏基因组-物种、基因等丰度信息与相关性数据文件每个元素之间的相关性。
工具链接:
https://www.omicshare.com/tools

详细教程链接:
《推荐一个绘制Science网络热图的在线工具》


   

   

     

6.漂亮的桑基图


   




桑基图(SanKey Plot)非常适合展示数据的流动变化,输入关联分析数据或其他数据绘制桑基图,查看数据的流向,或展示数据间的关联。
工具链接:
https://www.omicshare.com/tools/home/report/report_sankey.html

详细教程链接:
《你的论文中还缺一张“NB闪闪”的美图?》


   

   

     

7.多组差异散点图


   




多组差异散点图,以组间差异倍数的对数值log2FC为纵坐标,以比较组的名称为横坐标,可一次性展示多个比较组的差异基因。一般而言,对于传统的火山图,一张图一次只能展示一个比较组的差异基因,而多组差异散点图可视作多组“火山图”,可一次性展示多个比较组的差异基因。


工具链接:


https://www.omicshare.com/tools/home/report/report_diff_scatter_chart.html

详细教程链接:
《分享一种“过于先进,不宜展示”的图表绘制方法》


   

   

     

8.自定义标签热图


   




热图,主要使用渐变颜色来展现数据在不同样本中的变化规律,工具除了支持为热图添加分组注释条,也可以指定展示特定的文字标签,即使绘制上百个基因的热图,基因标签也不会重叠在一起。

工具链接:
https://www.omicshare.com/tools/Home/Soft/heatmap

详细教程链接:
《如何绘制“春日风格”的个性标签热图?》


   

   

     

9.动态交互韦恩图


   




韦恩图(Venn diagram)用于展示不同数据集之间的交集情况,在韦恩图中一般以椭圆、正圆等图形代表不同的数据集,以图形之间的重叠区域代表数据集之间的交集情况。下面为大家推荐一个强大的韦恩图在线工具,可以轻松完成多分组韦恩图的绘制。

工具链接:
https://www.omicshare.com/tools/home/report/reportvenn.html

该工具支持通过“筛选后”的元素数据文件或“筛选前”的丰度数据文件绘制韦恩图,并支持查看和提取子集部分的元素信息。


详细教程链接:
《如何绘制漂亮的韦恩图,可以提取子集元素的那种?》


   

   

     

10.环状网络图


   




环状网络绘制方法其实比较简单,除了使用Cytoscape软件,使用ggraph包可实现网络图的绘制。

#载入所需的绘图包;
library(ggraph)
library(tidygraph)
library(igraph)
library(ggplot2)
#读入边文件和节点信息文件;
net <- read.csv("edge_data.csv")
info <- read.csv("node_info.csv")
#预览数据;
head(net)
head(info)

#创建网络图对象;
g <- tbl_graph(nodes = info,edges = net,directed = F)
#计算网络图的degree属性;
V(g)$degree <- degree(g)
#自定义颜色;
mycol <- c("#FF8901","#00C5FF","#FF5485")

#绘制环状网络图(曲线连线);
p1 <- ggraph(g, layout = 'linear',circular = TRUE)+
  geom_edge_arc(colour="grey50",width=1,alpha=0.3)+
  geom_node_point(aes(color=class),size=6,alpha=0.8)+
  scale_colour_manual(values = mycol)+
  geom_node_text(aes(x = 1.05 * x,y = 1.05 * y,
                     angle = -((-node_angle(x, y) + 90) %% 180) + 90,
                     label=name),
                 nudge_y = 0,
                 hjust = 'outward',
                 repel = F,
                 size=2.5)+
  coord_fixed(clip = "off")+
  theme_graph()
p1




详细教程链接:
《网络图画成这样,是什么水平?》


   

   

     

11.星号标记热图


   




一般在绘制常规的相关性热图后,根据pvalue数值继续在热图上添加显著性标记(星号),如下图。

详细教程链接:
《如何绘制带有显著性“星号”的热图?》


   

   

     

12.正负号标记热图


   




除了在热图上添加显著性标记,我们也可以使用“正负号”标记正相关和负相关,如下图,或者仅显示某一范围的相关性系数值。

详细教程链接:
《如何绘制带有“正负号”标记的热图?》


   

   

     

13.环形热图


   




环状热图可与环形聚类树结合使用,看起来很整洁美观,特别适合展示不同分组基因的表达数据。

详细教程链接:
《热图也开始“内卷”——环形热图如何绘制?》


   

   

     

14.网络韦恩图


   




网络韦恩图,用以展示不同分组之间共有或特有的基因,特别是当分组元素数量较少时,数据展示效果更加直观。

详细教程链接:
《如何绘制各种高颜值的网络韦恩图?》


   

   

     

15.好看的火山图


   




火山图(Volcano plot)是一种比较“远古”的一种散点图,广泛应用于转录组、蛋白组等组间差异分析结果的展示。图表的横轴一般展示差异倍数的变化,而纵轴表示差异分析结果的可靠性。


我这里对传统的火山图做了进一步的优化,比如通过调整点的大小突出展示感兴趣基因对应的点,并将感兴趣基因的名称以标签的方式展示出来。那么,如何绘制这般个性又好看的火山图呢?我这里主要用到的还是最常见的ggplot2包。

详细教程链接:


《如何绘制新颖好看的火山图?》




   

   

     

16. t-SNE散点图


   




在单细胞转录组数据分析中,t-SNE散点图是展示细胞亚群图谱的重要手段,而Seurat的标准绘图效果并不能做更多调整,不过我们可以提取数据使用ggplot2重新绘制!


详细教程链接:
《如何提取Seurat数据,绘制更加个性化的tSNE/UMAP图?》


   

   

     

17.双坐标轴图表


   




转录组测序后,往往需要用qPCR结果去验证,我们常常需要绘制双坐标轴图表进行比较,如下图,左侧的坐标轴对应Q-PCR的结果(柱状图),右侧的坐标轴对应RNA-seq的测序结果(折线图)。


详细教程链接:


《验证组学数据?试试绘制双坐标轴图表》


   

   

     

18.上下对称组合图


   




上下对称图表适合展示两组量纲不同却具有一定相关性的数据,如下图。图表的下部分表示一段时间内的降雨量(mm)变化,而上部分则表示相应时段对应的蒸发量(m3/s)。



来源:现代石油人
ACTTosca断裂非线性化学电子油气MATLABUG岩土UM裂纹理论材料控制曲面
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-05-07
最近编辑:6月前
现代石油人
博士 签名征集中
获赞 26粉丝 63文章 814课程 1
点赞
收藏
作者推荐
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈