基因表达模式聚类以及可视化

最近在使用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<