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 =  TRUE ), IRanges (start =  sample (1:300, size = N, replace =  TRUE ), width =  sample (70:75,size = N, replace =  TRUE )), strand =  sample ( c ( "+" ,  "-" ,  "*" ),size = N,replace =  TRUE ), value =  rnorm (N, 10, 3), score =  rnorm (N, 100, 30),sample =  sample ( c ( "Normal" ,  "Tumor" ), size = N, replace =  TRUE ), pair =  sample ( letters ,size = N, replace =  TRUE ))   library (ggbio) seqlengths (gr) <-  c (400, 500, 700) values (gr)$to.gr <- gr[ sample (1: length (gr), size =  length (gr))] idx <-  sample (1: length (gr), size = 50) gr <- gr[idx]   ggplot () +  layout_circle (gr, geom =  "ideo" , fill =  "gray70" , radius = 7, trackWidth = 3) + layout_circle (gr, geom =  "bar" , radius = 10, trackWidth = 4, aes (fill = score, y = score)) + layout_circle (gr, geom =  "point" , color =  "red" , radius = 14,trackWidth = 3, grid =  TRUE ,  aes (y = score)) + layout_circle (gr, geom =  "link" , linked.to =  "to.gr" , radius = 6, trackWidth= 1)