语言绘制基因表达热图_绘制火山图&热图

在上一篇文章中，我们已经对基因进行了差异分析，接下来我们根据结果中的FDR值和FC值筛选出上调基因和下调基因(上调基因：基因转录成mRNA时受到正向调控，促进表达；下调基因：转录成mRNA时受到抑制，表达量减少)，并绘制成火山图与热图。

所用工具：R语言；

所需要包：ggplot2、pheatmap。

第一部分：火山图

首先，加载所需的包并导入数据：

library(ggplot2)diff_stat 

其次，筛选上调趋势数据和下调趋势数据，对于Fold Change值和p值阈值的选择，还需在实际的分析中视情况而定，本文以|log2FC| ≥2以及FDR p-value < 0.05作为差异OTUs的判断依据：

diff_stat[which(diff_stat$FDR < 0.05 & diff_stat$logFC >= 2),'diff'] diff_stat[which(diff_stat$FDR < 0.05 & diff_stat$lo