|撰文:莫比
之前生信大神谢老师在基迪奥生信交流群展示了一张用R语言绘制的组合图,令群内无数粉丝疯狂不已。这张图是这样子的:
对于我来说,绘制这张图的难点是如何将3张图组合到一起,而且“右图”的绘制方向与“上图”相反才能和主图(散点图)的坐标轴一致。
此外还有直方图和累计频率点线图的组合、Offset axis、全画布辅助线的添加、图例的添加等,如果不精通R的绘图函数,寸步难行!
那么,除了R语言,常规的作图软件比如Prism可以做到吗?
于是,找谢老师要了范例数据准备验证一下。顺便讲一下,如果大家对全基因组关联分析(GWAS)感兴趣可以前往Omicshare课堂观看谢老师的系列视频教程。
我的作图思路是:分别绘制“主”、“上”、“右”三个图表,然后把三个工程文件merge成1个,使用Prism的layout功能将三个图表组合在一起,最后手动添加辅助线、图例即可。下面就一起看下如何绘制吧!
|“主图”绘制
数据准备的准备很简单,只需π_poolA/π_poolB和Fst这两列数据即可,然后筛选出π_poolA/π_poolB≥1.174且Fst≥0.069的数据,以及其余数据。为了方便快速演示,我这里仅选取了3万个位点数据。如果不会筛选作图数据,可参考《如何使用Excel对数据进行筛选?》一文的方法。
然后,将数据复制粘贴到Prism的XY表格中