如何使用Circos-5:2D Data Tracks-histogram 简单配置

说明,该系列原文写于2016年3月。

如何使用Circos(2) 我们曾经说过 2D Data Tracks的文件格式基本一致,同样,就配置参数而言各类2D Data Tracks也基本相同。

在 **简单配置 **这一部分内容我们只举柱状图和热度两例,其他类型的参数配置我们留在后面的 **进阶配置 **中再做进一步介绍。

例1:histogram

# 首先让我们来复习一下前面的内容

karyotype = data/karyotype/karyotype.human.txt
chromosomes_units = 1000000

chromosomes_display_default = no
chromosomes = /hs[1234]$/
chromosomes_radius = hs4:0.9r

<colors>
chr1* = red
chr2* = orange
chr3* = green
chr4* = blue
</colors>

chromosomes_reverse = /hs[234]/
chromosomes_scale = hs1=0.5r,/hs[234]/=0.5rn

# 再来复习一下原始数据格式
#
# #chr start end value [options]
# ...
# hs3 196000000 197999999 71.0000
# hs3 198000000 199999999 57.0000
# hs4 0 1999999 28.0000
# hs4 2000000 3999999 40.0000
# hs4 4000000 5999999 59.0000
# ...
#
# 每一个柱状图被定义在一个 <plot> block中, <plots> block可以包含若干个<blaock>

<plots>
<plot>
type = histogram
#首先要说明类型,比如histogram或者scatter(散点图)等
file = data/5/segdup.hs1234.hist.txt
#之后写清楚引用原始文件位置

# 每一个track 被定义在 r0/r1 的两个半径之间,需要说明的是,当你使用r这个相对位置时,它的位置时相对于ideogram而言的
#在之前内容中我们曾经提到过一次染色体显示比例的相对位置,那时的r是相对于整个图形的周长而言的。
r1 = 0.88r
r0 = 0.81r

# 柱状图既可以有填充fill 也可以有轮廓 outline. 默认设置其 outline 为 1px thick black.
fill_color = vdgrey
#vdgrey的意思是very deep,这一部分内容会在和颜色相关的部分进行介绍。
#如果不想要轮廓,可以将outline的thickness设置为0.

#thickness = 0p
# Do not join histogram bins that do not abut .
extend_bin = no

# 和link相似,这里也可以加入相关的rule
<rule>
condition = on(hs1)
show = no
</rule>
# 不显示hs1的任何数据

</plot>

#在这个例子中,作者还介绍了一个堆叠柱状图的概念,即同一个位置可以展示不同的value,
#然后给每一个value赋予不同的颜色以区别显示。

#原始文件格式
#chr start end value,value,value,... [options]
# ...
# hs3 196000000 197999999 0.0000,7.0000,64.0000,0.0000
# hs3 198000000 199999999 21.0000,6.0000,18.0000,12.0000
# ...
#当 type=histogram 并且原始文件中有多个value时,circos会自动绘制堆叠的柱状图。

<plot>
type = histogram
file = data/5/segdup.hs1234.stacked.txt
r1 = 0.99r
r0 = 0.92r

#绘制上述堆叠状柱状图时,最好采用一系列的颜色,这里使用的是默认设置中给定的人染色体颜色。
# hs3 198000000 199999999 21.0000,6.0000,18.0000,12.0000
# the colors for the bins will be hs1 (21.0000), hs2 (6.0000), hs3 (18.0000), hs4 (12.0000).

fill_color = hs1,hs2,hs3,hs4
extend_bin = no

# 柱状图的延伸方向可以单独进行设置。
orientation = in
#修改为no则向外,具体图示可以参考example4中的circos1.png

图例


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

推荐阅读更多精彩内容