circos绘制彩色的links

对于links而言,默认情况下所有的links 都是同一种颜色,通过color属性指定。当我们想要构建出彩色的links 时,通常情况下有3种做法:

1. 拆分file

每一个link都有一个file指定的数据,我们可以人为的将原来的file根据条件拆分成多份,每一份对应一个link, 然后为每个link设置不同的颜色,就可以实现彩色的link了。

示例如下:

每个link定义一个文件,对应的color值不同,生成的图片如下

2. 定义rules

拆分数据比较繁琐,一种更加简便的方法是定义rules, 根据不同的规则指定不同的颜色。示例如下

需要理解的就是condition的写法。第一种就是var的用法,var获取指定字段的值,在links中,有以下几种字段

  1. CHRn
    表示link 对应的染色体的名字,一个link连接两个区域,第一个区域对应的染色体为chr1, 第二个区域对应的染色体为chr2, condition中对应的写法为var(chr1), var(chr2)

  2. STARTn
    表示区域的起始位置,第一个区域的起始位置为start1, 第二个区域对应的起始位置为start2, condition中对应的写法为var(start1), var(start2)

  3. ENDn
    表示区域的终止位置,第一个区域的终止位置为end1, 第二个区域对应的终止位置为end2, condition中对应的写法为var(end1), var(end2)

  4. POSn
    表示区域的中心点的位置,第一个区域的中心位置为position1, 第二个区域对应的中心位置为position2, condition中对应的写法为var(position1), var(position2)

  5. SIZEn
    表示区域的长度,第一个区域的长度为size1, 第二个区域对应的长度为size2, condition中对应的写法为var(size1), var(size2)

  6. REVn
    判断某个区域是否反向,如果起始位置大于终止位置,代表是反向的,返回值为1,否则返回值为0, condition中对应的写法为var(rev1),var(rev2)

  7. INV
    如果一个link连接的两个区域方向不同,一个正向,一个反向,返回值为1,其他情况返回值为0,condition中对应的写法为var(inv)

  8. INTERCHR
    如果一个link连接的两个区域位于两条染色体上,返回值为1,其他情况返回值为0,condition中对应的写法为var(interchr)

  9. INTRACHR
    如果一个link连接的两个区域位于同一条染色体上,返回值为1,其他情况返回值为0,condition中对应的写法为var(intrachr)

var会返回对应的值,从上面可以看出,返回值可以分成以下3类

  • 字符串

  • 数字

  • 逻辑值

对于字符串,采用perl当中的字符串操作符,示例

condition = var(chr1) eq “hs1”

对于数字,采用数字操作符, 示例

condition = var(size1) > 10mb

对于逻辑值, 可以直接使用,示例

condition = var(intrachr)

除了使用var 之外,还可以使用函数,condition中支持以下函数

  1. conf

  2. on

  3. within

  4. between

  5. fromto

  6. tofrom

  7. from

  8. to

  9. chrlen

links中,支持between的用法, 示例如下

condition = between(hs1, hs2)

rules的示例:

生成的图片如下:

3.在file中添加属性

file文件支持内置属性和自定义属性,所有的属性写在最后一列,多个属性用逗号分隔。

示例如下:

上面的文件中,最后一列包含了colorvalue两种属性,color是内置属性,value是自定义属性。通过在file中添加不同取值的color属性,可以方便的实现彩色的links。这里主要看下通过value属性的值映射到不用颜色上。

配置文件写法如下:

通过var(value)获取每个link的value 值作为数组的下标,数组由qw定义,数组中的元素为5个不同的颜色,value的取值范围为0-5,正好对应该数组的下标,通过这种方式,将value映射到颜色上;而且在原来颜色的基础上,还通过sprintf函数,为原来的颜色添加_a3后缀,这个后缀代表一定的透明度。

生成的图片如下


根据自己的数据,灵活运用以上3种方法,就可以随心所欲的构建多彩的 links了。



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

推荐阅读更多精彩内容

  • Spring Cloud为开发人员提供了快速构建分布式系统中一些常见模式的工具(例如配置管理,服务发现,断路器,智...
    卡卡罗2017阅读 134,644评论 18 139
  • 第2章 基本语法 2.1 概述 基本句法和变量 语句 JavaScript程序的执行单位为行(line),也就是一...
    悟名先生阅读 4,138评论 0 13
  • HTML 5 HTML5概述 因特网上的信息是以网页的形式展示给用户的,因此网页是网络信息传递的载体。网页文件是用...
    阿啊阿吖丁阅读 3,875评论 0 0
  • rules是circos中的一种黑魔法,可以在不改变数据的情况下,对展现形式进行调整。rules代表规则,本质上就...
    生信修炼手册阅读 1,377评论 0 3
  • 《Signal 信号》韩剧原创剧评 第二篇 – 第13-16集观感 请不要轻言放弃 “绝对不要放弃。过去是可以改变...
    隽永剧评阅读 712评论 0 6