本文首发于基迪奥生物公众号,如需转载,请先留言~
火山图的本质是散点图,筛选数据后用Excel其实也能画,不过可能需要花点时间。而用Origin去完成数据的生成和颜色的映射会更方便一点。
今天主要是想借绘制火山图这一案例,为大家介绍Origin数据的计算生成、“速度模式”和参考线的添加。
数据准备
本例所用的数据是来自周老师ggplot2视频教程中的示例数据(点击 阅读原文 下载),如下。数据有47430行,首先要做的是删除一些不要的列(数据量较大,只需保留log2FC、FDR、significance这3列),然后另存为Excel的专用格式(xlsx)。
然后用Origin直接导入外部数据,方法如下。
对导入的数据稍作整理,然后新增1列空列,方法如下。
然后选中这一列,通过Analysis/Mathematics/Set Column Values,给新列计算相应的-log10FDR值。
在计算窗口中输入公式:log(Col("FDR"))*-1,如下图,点OK即可完成计算,非常方便。
计算的结果如下图,给这列新数据自定义列名(Long Name),我这里简单命名为y,数据的准备工作就完成啦。
图表绘制
通过Plot/Scatter/Color Mapped绘制散点图,如下图。
在Plot Setup窗口,给X轴、Y轴、映射颜色选择要绘制的数据,方法如下,点OK按钮完成绘制。
由于数据量较大(接近5万多个散点),Origin会提示自动进入Speed Mode(也就是仅显示一部分数据点,便于反复自定义调整而不必每次都进行大量的运算),初始的结果如下:
双击图表,在Plot Details窗口可调整散点的形状,大小,透明度等,如下。
还可以在Color List选项下自定义散点的颜色,如下图。
改完颜色后,右键图例,通过Legend\Reconstruct Legend更新一下图例,并改一下y轴的标题,得到的初始效果如下图。
添加参考线
双击坐标轴,在坐标轴调节窗口的Reference Lines选项卡下,先添加垂直的参考线,Reference Linesat Value输入框中输入-1 1(空格分隔),如下图,点Insert按钮,再点Apply按钮完成数据垂直方向参考线的添加。
点Details按钮可对参考线的颜色、线的类型、粗细等进行自定以设置。这里可同时选中-1、1两条线进行同时设置,如下图。
同样的方法,也可以添加水平的参考线,我这里输入1.301(LOG10(0.05)*-1)得到水平的参考线。同时,我也调整了下y轴和x轴的显示范围,结果如下图。调整满意后,点右侧的“速度模式”开关按钮(红色箭头所指),关掉Speed mode,让更多的点显示出来。
最终的效果如下:
更多的Origin视频教程可到Omicshare课堂观看,今天的内容就到这里啦~
为了更好的阅读体验,建议关注基迪奥生物公众号~
网友评论