今天去地里干活了，下午7点多才回到学校。所以今天的推文就没有新的学习笔记啦，翻了翻公众号很早之前发过的推文，找出来再重新发一下。

原文地址

https://stats.biopapyrus.jp/r/graph/circos-plot.html

代码

library(circlize)
library(RColorBrewer)
col<-RColorBrewer::brewer.pal(6,"Paired")
df<-data.frame(chr=sample(c(paste0("chr",1:6)),1000,replace = T),
               x=rnorm(1000),
               y=rnorm(1000))
head(df)
circos.initialize(factors = df$chr,x=df$x)
circos.trackPlotRegion(factors = df$chr,y=df$y,
                       panel.fun = function(x,y){
                         circos.axis()
for(i in 1:6){
  highlight.sector(sector.index = paste0("chr",i),col=col[i])
  circos.text(CELL_META$xcenter, CELL_META$ycenter,
              labels = paste0("chr",i),sector.index = paste0("chr",i))
circos.trackHist(df$chr,df$x,col=col)
circos.link("chr1",0,"chr2",0)
circos.link("chr4",c(-2,0),"chr3",c(-2,2),col="blue")
circos.link("chr4", c(0, 2), "chr6", c(-2, 0), col = "orange")
circos.link("chr5", c(0, 2), "chr1", c(-2, 2), col = "lightgrey")
circos.clear()