2018-11-7 基因結(jié)構(gòu)圖繪制之gggenes
最重要的包參考鏈接補(bǔ)上
library(ggplot2)
# install.packages("gggenes")
# devtools::install_github("wilkox/gggenes")
library(gggenes)
p1 <- ggplot(example_genes, aes(xmin = start, xmax = end, y = molecule, fill = gene)) +
geom_gene_arrow() +
facet_wrap(~ molecule, scales = "free", ncol = 1) +
scale_fill_brewer(palette = "Set3")
ggsave("test.pdf", p, width = 10, height = 20)
image.png
注意這里繪圖需要將Rstudio的繪圖區(qū)拉到最高（盡量高，否則會(huì)報(bào)錯(cuò)）
- 為什么會(huì)知道了藕各，參考了此鏈接
Error in grid.Call(C_convert, x, as.integer(whatfrom), as.integer(whatto), :
Viewport has zero dimension(s)
image.png
p2 <- ggplot2::ggplot(example_genes, ggplot2::aes(xmin = start, xmax = end, y =molecule, fill = gene)) +
geom_gene_arrow() +
ggplot2::facet_wrap(~ molecule, scales = "free", ncol = 1) +
ggplot2::scale_fill_brewer(palette = "Set3") +
theme_genes()
dummies <- make_alignment_dummies(
example_genes,
ggplot2::aes(xmin = start, xmax = end, y = molecule, id = gene),
on = "genE"
)
image.png
geom_gene_arrow() +
ggplot2::geom_blank(data = dummies) +
ggplot2::facet_wrap(~ molecule, scales = "free", ncol = 1) +
ggplot2::scale_fill_brewer(palette = "Set3") +
theme_genes()
image.png
p4 <- ggplot2::ggplot(example_genes, ggplot2::aes(xmin = start, xmax = end, y =molecule, fill = gene, label = gene)) +
geom_gene_arrow(arrowhead_height = unit(3, "mm"), arrowhead_width = unit(1, "mm")) +
geom_gene_label(align = "left") +
ggplot2::geom_blank(data = dummies) +
ggplot2::facet_wrap(~ molecule, scales = "free", ncol = 1) +
ggplot2::scale_fill_brewer(palette = "Set3") +
theme_genes()
image.png
example_genes$direction <- ifelse(example_genes$strand == "forward", 1, -1)
p5 <- ggplot2::ggplot(subset(example_genes, molecule == "Genome1"),
ggplot2::aes(xmin = start, xmax = end, y = strand, fill = gene,
forward = direction)) +
geom_gene_arrow() +
theme_genes()
image.png
