最近在使用RNA_seq数据做些分析,结果得到了大量差异表达以及共表达的基因,如何合理展示这些基因也是一件不简单的事情。除了常见的热图(heatmap)展现形式,今天在推荐另外一种展示方式(下图C)。
需要R包TCseq或者Mfuzz。我这里给出的代码是基于TCseq。
输入文件就是一个基因表达量的矩阵,如下图。
代码也很简单,见下图,也请点击阅读原文查看代码。
library(TCseq)
svg(filename="/stress_analysis/drought_and_heat/heat.svg",width=3.5,height=15)
data <- read<