使用ggbio在R中制作弦圖


分享一個制作弦圖的R包:ggbio. 以下是一個簡單的使用實例,效果圖和代碼如下。

 

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)

 

 

ggbio包還可以制作其他復雜的圖表(如下),具體使用方法可參考該包的幫助文檔:


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM