CubicHeatMap - 重磅推出 - 魔方热图 - 伪三维


Emmm, 今天的主图,

写在前面

博士毕业论文中,使用的样品包括:

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

TBtools中多种热图可视化方式

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

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


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

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

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

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


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

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

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

https://upload-images.jianshu.io/upload_images/10518391-a3309b0df0063b21.gif

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

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

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


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

Emmm.. 忽略配色,毕竟不一定能拿到原始的配色方案。这张图片,包括了 16个基因,11个组织和16个时期。

如何实现?

前述,在公众号推文中已经提及,JIGplot 引擎中有类似图层的概念,也就是 SubPanel。虽然 CubicHeatMap 看起来是立方体。不过我们完全可以拆解。

EJAR~BB)IE8QARMH_S}@DEK.png

JIGplot 除了掰弯,其实还有一个功能切变。后者的实现不是很健全,不过处理 Rectangle 完全没问题。直接画三个正常的热图,然后分别对 UpPanel 和 RightPanel 做 Shear 即可。

如何使用 CubicHeatMap ?

虽然是顺手写,不过还是打成GUI,方便自己后续使用【毕竟Java码不是R码,随取随用 】。于是,TBtools又增加了一个功能。

打开之后,可以看到,只有两个输入:


1. 基因表达矩阵,这个比较简单,跟其他热图完全一样


2. 样品分组信息,这个说起来复杂,其实也简单。一共三列:
第一列,样品名字,即表达矩阵对应的样品名字【注:没有格式和内容要求,对应起来就行】
第二列,样品的第一维度分类,比如 属于某个组织
第三列,样品的第二维度分类,比如 处于某个时期
如下:

两个文件放进去,自己设置下喜欢的配色,随后点击 Start 即可


于是得到


当然,我们不能忘了,JIGplot 是可交互的

礼物?

Emmm... 突然想起来,TBtools 的 Logo,其实就像一个热图。


既然有了CubicHeatMap, 那么我们完全可以做一组数据,就画 TBtools 的 Logo

写在后面

Emmm...
练手完毕,继续....
相信今年的 2月14 是不少人可能在很多年以后 ,仍然觉得特殊的日子。喜欢过节的,在往常可能会花一点时间去准备礼物,但是今年应该就不行了。
但这一切,应是没有所谓。毕竟
始终如一是最好的相处,而
陪伴就是最好的礼物。

祝大家身体健康了。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,657评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,662评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,143评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,732评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,837评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,036评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,126评论 3 410
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,868评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,315评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,641评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,773评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,470评论 4 333
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,126评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,859评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,095评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,584评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,676评论 2 351