R语言绘图--高级图形之Circos

Circos图：

R包：Rcircos

支持Circos 2D数据轨道绘图，例如散点图、线、柱状图、连接器等

实现：

library(RCircos)
library(Cairo)

outF<-'/Output/Circos.pdf'
data(UCSC.HG19.Human.CytoBandIdeogram)
data(RCircos.Scatter.Data)
data(RCircos.Histogram.Data)
data(RCircos.Link.Data)
data(RCircos.Ribbon.Data)
RCircos.Scatter.Data$chromosome<-paste('chr',RCircos.Scatter.Data$chromosome,sep = '')

###########
#chr.exclude<-NULL #设置不显示的染色体，如 c(1,3)
#tracks.inside <-10 #设置内部track 个数
#tracks.outside<-0 #设置外部track 个数
RCircos.Set.Core.Components(cyto.info=UCSC.HG19.Human.CytoBandIdeogram,chr.exclude=NULL,tracks.inside=5,tracks.outside=0)

CairoPDF(outF,width = 8,height = 8)
RCircos.Set.Plot.Area()
RCircos.Chromosome.Ideogram.Plot()
RCircos.Scatter.Plot(RCircos.Scatter.Data, data.col=5,track.num=1, side='in', by.fold=1)

RCircos.Histogram.Plot(RCircos.Histogram.Data,data.col = 4,track.num = 3,side = 'in')

RCircos.Link.Plot(RCircos.Link.Data,track.num=4,TRUE)
RCircos.Ribbon.Plot(ribbon.data=RCircos.Ribbon.Data, track.num=5, by.chromosome=FALSE, twist=FALSE)
dev.off()

图形：