使用希尔伯特曲线(Hilbert Curve)可视化基因组数据。
基本信息
| 软件包 | 链接 |
|---|---|
| 编程语言 | R |
| Bioconductor | https://bioconductor.org/packages/HilbertCurve/ |
| GitHub | https://github.com/jokergoo/HilbertCurve |
| 文档 | https://jokergoo.github.io/HilbertCurve/ |
| 论文 | Zuguang Gu, et al., HilbertCurve: an R/Bioconductor package for high-resolution visualization of genomic data. Bioinformatics 2016. |
例子
人类染色体1上的基因分布。
library(GenomicRanges)
library(HilbertCurve)
library(circlize)
load(system.file("extdata", "refseq_chr1.RData", package = "HilbertCurve"))
hc = GenomicHilbertCurve(chr = "chr1", level = 5, reference = TRUE,
reference_gp = gpar(lty = 1, col = "grey"), arrow = FALSE)
hc_segments(hc, g, gp = gpar(lwd = 6, col = rand_color(length(g))))