R/RCircos.R

   1 datN=read.delim('Normal.cg',header=F)
   2 datB=read.delim('brother2.cg',header=F)
   3
   4 library(RCircos)
   5 data(UCSC.HG19.Human.CytoBandIdeogram);
   6 cyto.info <- UCSC.HG19.Human.CytoBandIdeogram;
   7 chr.exclude <- NULL;
   8 num.inside <- 5;
   9 num.outside <- 0;
  10 #chr.exclude = paste("chr", seq(5,22), sep="")
  11
  12 RCircos.Set.Core.Components(cyto.info,chr.exclude, num.inside, num.outside);
  13 RCircos.Set.Plot.Area();
  14 RCircos.Chromosome.Ideogram.Plot();
  15
  16 RCircos.Histogram.Plot(hist.data=datB,data.col=6, track.num=1, side="in");
  17 RCircos.Histogram.Plot(hist.data=datN,data.col=6, track.num=2, side="in");
  18
  19 datN['D'] <- datB$V6-datN$V6
  20 datN['X'] <- 20 * datN['D']
  21 RCircos.Histogram.Plot(hist.data=datN,data.col=8, track.num=3, side="in");
  22 RCircos.Histogram.Plot(hist.data=datN,data.col=7, track.num=5, side="in");
  23
  24 # https://github.com/cran/RCircos/blob/master/demo/RCircos.Demo.Human.R
  25 data(RCircos.Histogram.Data);
  26 RCircos.Histogram.Plot(hist.data=RCircos.Histogram.Data,data.col=4, track.num=6, side="in");