【画图专题】火山图Volcano Plot--差异表达分析结果展示系列1

一般做完差异表达分析后，会有一堆的结果与数据。那么，怎么样让你的分析结果在paper中以更加直观形象（好看炫酷）的方式展示出来呢？前面我们已经出过【画图专题】史上最详细画热图函数pheatmap()参数讲解。下面我们教你另外一种展示方式，火山图（Volcano Plot）。

首先，我们要注意两个重要概念，

FDR：在差异表达分析过程中，采用统计学方法对原有假设检验得到的p值进行校正，并最终采用校正后的p值，即FDR（False Discovery Rate）。

FC（Fold Change）：是两样品（组）间基因表达水平的比值，是表示差异倍数的变量。一般差异表达分析中会同时控制这两个参数来筛选显著差异表达基因，如取FDR<0.05和FC>=2。

火山图，它可以直观地展现所有基因的FDR与Fold Change之间的关系，以便快速查看基因在两组样品间表达水平的差异程度和其统计学显著性。

火山图的本质其实就是散点图。

横坐标表示某一个基因在两样品中表达水平差异倍数的对数值，即log2(FC);横坐标绝对值越大，说明表达量在两样品间的表达量倍数差异越大。

纵坐标表示pvalue 的负对数值，即-log10(pvalue)。纵坐标值越大，表明差异表达越显著，筛选得到的差异表达基因越可靠。

1，数据准备

差异表达分析结果数据示例：第一列为基因名，第二列为差异表达分析算法得到的p值统计量，第三列为矫正后的p值，即FDR，第四列为log2FC值。画火山图只需要最后两列就行。

2，R代码画图

最后画出来的火山图如下，是不是比你干巴巴的在文章里面贴一个差异表达分析结果的表格好看多了！！！

此文章同步更新在

画图专题 | 火山图Volcano Plot