CubicHeatMap - 重磅推出 - 魔方热图

Emmm, 今天的主图，

写在前面

博士毕业论文中，使用的样品包括：

不同的类型 / 或者说不同的组织
不同的时间 / 或者说时间序列
所以，是一个多X多的设计。对于这类实验设计所产生的数据，往往会有较高的分析难度。在分析上，现在多数人直接使用的是共表达网络分析，最常用的莫过于 WGCNA 。而在可视化上，单单就基因表达热图绘制上，没有较好的形式，尤其是在多基因X多组织X多时期的数据上。

TBtools 中多种热图可视化方式

TBtools中，我们提供过数种热图的展示方式，

第一种，常规的热图，适用于少基因X少样品的数据

这种热图，四四方方，其实在TBtools里面，可以有不少变化，比如不要用颜色对应数值，直接用圈的大小来对应数值

当然还有其他，比如三角形的，或者前述我们推送过的对角线热图的热图组合方式。无论怎么可视化，这类热图能展示的表达矩阵，仍然很少。所以在基因数达到数十或者上百之后，TBtools中有提供了掰弯的热图或者说，环形热图。

第二种，环形的热图，适用于多基因X少样品的数据

前前后后，我在公众号上推环形热图有超过一年时间。目前市面上，我自认为暂时就TBtools一个工具，也是第一个，是以热图工具为中心开发的环形热图，如果说之前已经有，那么那些都是以进化树注释为中心开发的。进化树注释功能，TBtools也有，对应的也有简单版本的热图Track【注：功能明显有限】。如果说，还有其他工具可以做了，那么他们都只是追随者。环形热图，最大的好处，就是可以容纳比较多的基因，同时基本不占版面。
大体如下，

或者，换个配色可能好看点

第三种，卡通热图，适用于少基因X多样品的数据

虽然掰弯的热图，其实把基因和样品做一个转置【注：TBtools的热图功能内置矩阵转置功能，也有专门的表格转置功能】，就可以处理少基因X多样品的情况。但是，当样品很多，或者比较复杂时。我们人类的直观视觉可以帮我们更好的分析数据。为此，在多年前，有一个浏览器框架叫eFP Browser，针对芯片数据开发。在前述，TBtools推出了基于FPKM的，支持矢量图导入导出的FancyHeatmap。从比较多的方面有所超越。大体如下：

这种热图，最适合的是可视化某个基因在多组织多时期的表达变化。

CubicHeatMap 可视化多基因X多组织X多时期

前几天跟朋友在讨论的课题的时候【注：大家都是做植物发育生物学的】，讨论到几年前的一篇Ncomm文章，其中有图如下：

一个非常不错的可视化方式，虽然也带有缺陷【注：只有一个基因的表达量可以为完全展示】。不过这个缺点不影响他的实用性，尤其是我们一系列共表达量基因的时候。很不幸，这个工具，似乎仍然是一个浏览器框架，Tomato Expression Altas。事实上，我也不知道这个浏览器框架到底有没有另外的人用过，不过我是不会用的。
想想很久没有Coding，手和脑子都会生疏。于是，一般看娃，干脆就一边实现这个热图功能。结果图如下：