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