R做圈图

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 library (GenomicRanges)   set.seed (1) N <- 100   gr <-  GRanges (seqnames =  sample ( c ( "chr1" ,  "chr2" ,  "chr3" ), size = N, replace =  T

>>阅读原文<<