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()