Кажется, что функция manhattan
из пакета qqman
не совместима с ggplot2:ggsave
, поскольку использует базовую графику, а не grid
графику.Если хотите, проверьте код функции manhattan
с помощью getAnywhere(manhattan)
.
Попробуйте вместо этого:
library(qqman)
library(grid)
library(gridGraphics)
library(ggplot2)
# your plot, but without the ggplot2 theme lines (they do nothing anyways)
manhattan(gwasResults, chr="CHR", bp="BP", snp="SNP", p="P",
col = c("chartreuse2", "darkorange1", "gold1"),
ylim = c(0, -log10(1e-06)), chrlabs = NULL,
suggestiveline = -log10(1e-03), genomewideline = -log10(1e-05),
highlight = NULL, logp = TRUE, annotatePval = NULL,
annotateTop = TRUE, main='DWStem')
# Transform the base plot to grobs, so that is compatible with grid and, therefore ggsave
p <- recordPlot()
g <- grid.grabExpr(grid.echo(p))
ggsave("DWSte.png", g, bg = "transparent")
В качестве альтернативы, этот учебник , кажется, даетggplot2
решение.