circos 可视化手册- heatmap 篇

heatmap将一系列的数值映射到一个颜色梯度中,是最常用的图表之一。在circos中,通过plot这个block进行设置。

用法如下:

heatmap的数据和scatter plot, histogram是一样的,最少有4列内容。前3列定义染色体上的区域,第4列指定每个区域的value值,其他列可以是配置文件中的属性,示例如下

heatmapcolor属性比较特殊,在热图中,需要的是一个颜色梯度,将不同的value映射到这个颜色梯度中。所以这里的color对应的是一个list, 有两种设置方式

1. 采用内置的list

etc目录下,有很多内置的颜色配置,其中就有很多内置的颜色梯度, 示例如下

etc/colors,brewer.conf

blues-3-seqblues-3-seq-rev就是两个颜色提取的名字,对应有3种颜色

  • blues-3-seq-1

  • blues-3-seq-2

  • blues-3-seq-3


配置文件的写法如下

还可以将多个颜色梯度结合起来使用

2. 采用自定义的颜色

写法如下

在热图中,最重要的是理解valuecolor之间的对应关系。在heatmap中,共有两种映射关系

1.线性映射

默认情况下,valuecolor之间就是线性映射关系。color 的值是一个list, 对于每个value, 通过计算对应的下标color_index,就可以得到具体的颜色。

下标的计算公式如下:

从公式可以看出,二者之间是一个线性映射关系。当我们没有在配置文件中指定maxmin的值时,circos会自动根据file文件的第4列计算这两个值。当然,我们也可以手动在配置文件中指定。

当我们指定的minmax与数据不符合, 比如数据的范围是[2, 6], 配置文件中指定的是[1,5], 通过以上公式计算时,会出现下标越界的情况,此时,对于小于min的值,其颜色定义为list中的第一个颜色值;对于大于max的值,其颜色定义为list中的最后一个颜色值。

具体的是通过color_mapping参数控制线性映射的属性,共有3种取值

  1. color_mapping = 0

  2. color_mapping = 1

  3. color_mapping = 2


默认情况下,color_mapping = 0color_mapping的值设置的是一个偏移量,对于color设置的n种颜色,每种颜色均分成两份,从0开始计数,00..nncolor_mapping = 1 时,向后偏移一位,变成011..nn; color_mapping = 2时,向后偏移两位,变成112233

反应到颜色上,具体的示意图如下:

2. 非线性映射

在配置文件中添加下面的参数,开启log映射

当数据分布不均匀时,采用log映射更加合适。

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

推荐阅读更多精彩内容