主讲老师第六十五课:ceRNA 的桑基图绘制
桑基图(Sankey diagram),当然,也有人根据其形状称为“决策树”, 或者因其流向称为“流程图”。桑基图由于其独特的特征,能够很好的反映组间变量的相互作用关系:1). 起始流量和结束流量相同,也就是我们常说的能量守恒;2).在图形内部,不同的线条代表了不同的分流情况,线条的粗细与其占比成正比例。目前,除了相关性分析外,桑基图也常常用于 ceRNA 的可视化方式。下面,我们来看下如何基于 ceRNA 调控网络进行桑基图的绘制(ceRNA_plot.zip)
1.读取数据
首先,分别读取分析得到的 miRNA_mRNA 和 lncRNA_miRNA 之间的交互关系
2.整合数据
通过 miRNA,使用 merge() 函数将两个交互网络进行合并。
随后,将宽数据转化为长数据。
3.桑基图的绘制
3.1 指定颜色
首先,我们需要对每个基因内容分配颜色,这里使用 rainbow() 彩虹颜色函数。
3.2 绘制桑基图
首先,加载绘图使用的 R 包,在此,我们使用 ggalluvial 包来绘制桑基图。
基于 ggplot()函数,使用 geom_stratum()来指定桑基图中的堆叠柱形图在文章里面,习惯了 Cytoscape 的结果后,偶尔用用这样的结果图也是不错的