Circlize绘制飞镖靶

今天用到的颜色

开始作图-初始化：取消track间距，取消gap间距。

myfactor<- 1:20 
library(circlize)

## ========================================
## circlize version 0.4.4
## CRAN page: https://cran.r-project.org/package=circlize
## Github page: https://github.com/jokergoo/circlize
## Documentation: http://jokergoo.github.io/circlize_book/book/
## 
## If you use it in published research, please cite:
## Gu, Z. circlize implements and enhances circular visualization 
##   in R. Bioinformatics 2014.
## ========================================

circos.par( gap.degree = 0 , track.margin = c(0,0), start.degree = 360/40)
circos.initialize(factors = myfactor, xlim = c(0,1))

circos.clear()

加个外圈，占15%，填入文字，边界黑色

myfactor<- 1:20 
library(circlize)
circos.par( gap.degree = 0 , track.margin = c(0,0), start.degree = 360/40)
circos.initialize(factors = myfactor, xlim = c(0,1))
circos.trackPlotRegion(ylim = c(0,1) , bg.col= "black" , bg.border = 1, track.height = 0.15, panel.fun = function(x,y){circos.text(CELL_META$xcenter, CELL_META$ycenter , CELL_META$sector.index ,col = "white" , font = 2)})

circos.clear()

加内圈红绿色，占5%

myfactor<- 1:20 
library(circlize)
circos.par( gap.degree = 0 , track.margin = c(0,0), start.degree = 360/40)
circos.initialize(factors = myfactor, xlim = c(0,1))
circos.trackPlotRegion(ylim = c(0,1) , bg.col= "black" , bg.border = 1, track.height = 0.15, panel.fun = function(x,y){circos.text(CELL_META$xcenter, CELL_META$ycenter , CELL_META$sector.index ,col = "white" , font = 2)})
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#E41A1C', '#4DAF4A'),10) , track.height = 0.05)

circos.clear()

以此类推加到最后一个

myfactor<- 1:20 
library(circlize)
circos.par( gap.degree = 0 , track.margin = c(0,0), start.degree = 360/40)
circos.initialize(factors = myfactor, xlim = c(0,1))
circos.trackPlotRegion(ylim = c(0,1) , bg.col= "black" , bg.border = 1, track.height = 0.15, panel.fun = function(x,y){circos.text(CELL_META$xcenter, CELL_META$ycenter , CELL_META$sector.index ,col = "white" , font = 2)})
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#E41A1C', '#4DAF4A'),10) , track.height = 0.05)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#FFFFFF', '#000000'),10), track.height = 0.275)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#E41A1C', '#4DAF4A'),10) , track.height = 0.05)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#FFFFFF', '#000000'),10), track.height = 0.375)
circos.trackPlotRegion(ylim = c(0,1) ,bg.border = NA , bg.col = "#4DAF4A" , track.height = 0.05)

circos.clear()

最后一个空白的洞不能用circos.track,而应该用draw.sector

myfactor<- 1:20 
library(circlize)
circos.par( gap.degree = 0 , track.margin = c(0,0), start.degree = 360/40)
circos.initialize(factors = myfactor, xlim = c(0,1))
circos.trackPlotRegion(ylim = c(0,1) , bg.col= "black" , bg.border = 1, track.height = 0.15, panel.fun = function(x,y){circos.text(CELL_META$xcenter, CELL_META$ycenter , CELL_META$sector.index ,col = "white" , font = 2)})
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#E41A1C', '#4DAF4A'),10) , track.height = 0.05)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#FFFFFF', '#000000'),10), track.height = 0.275)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#E41A1C', '#4DAF4A'),10) , track.height = 0.05)
circos.trackPlotRegion(ylim = c(0,1) , bg.border = "0", bg.col = rep(c('#FFFFFF', '#000000'),10), track.height = 0.375)
circos.trackPlotRegion(ylim = c(0,1) ,bg.border = NA , bg.col = "#4DAF4A" , track.height = 0.05)
draw.sector(center = c(0, 0), start.degree = 0, end.degree = 360,
    rou1 = 0.05, col = "#E41A1C", border = "#EEEEEE")

circos.clear()

circlize的这一系列都是模仿jokergoo在github上提供的example.