使用“gridExtra”和多个方面并排绘制图形

Question

我有两个不同大小的数据库，dt 和 dt1。我想使用 gridExtra 包中的命令 grid.arrange 并排显示 g1 和 g2。如果可能的话，我还希望看到 g1 和 g2 使用 facet_grid 或 facet_wrap 命令或使用 gridExtra 但使用 facet_grid\facet_wrap 视觉。我在互联网上进行了长时间的搜索，但无法使用下面的代码获取这些图形。

set.seed(000)
m <- matrix(rnorm(1000,0,1),1000,1)
dt <- data.frame(m)
names(dt) <- c("X")

library(ggplot2)

g1 <- ggplot(dt, aes(x=X)) 
g1 <- g1+geom_histogram(aes(y=..density..),      # Histogram with density instead of count on y-axis
                        binwidth=.5,
                        colour="black", fill="white",breaks=seq(-2, 2, by = 0.1)) 
g1 <- g1 + stat_function(fun=dnorm,
                         color="black",geom="area", fill="gray", alpha=0.1,
                         args=list(mean=mean(dt$X), 
                                   sd=sd(dt$X)))
g1 <- g1+  geom_vline(aes(xintercept=0, linetype="Valor Verdadeiro"),show.legend =TRUE)
g1 <- g1+  geom_vline(aes(xintercept=mean(dt$X, na.rm=T),    linetype="Valor Estimado"),show.legend =TRUE)
g1 <- g1+  scale_linetype_manual(values=c("dotdash","solid")) # Overlay with transparent density plot
g1 <- g1+  xlab(expression(paste(gamma[1])))+ylab("Densidade")
g1 <- g1+  theme(plot.margin=unit(c(0.5, 0.5, 0.5, 0.5), units="line"),
                 legend.position = "top",
                 legend.justification = c(0,0),
                 legend.box.just = "top",
                 legend.margin = margin(0,0,-10,-5),
                 legend.title=element_blank(),
                 legend.direction = "horizontal",
                 legend.background = element_rect(fill="transparent", size=.5, linetype="dotted"))
g1 <- g1+ guides(linetype = guide_legend(override.aes = list(size = 1)))


# Adjust key height and width
g1 = g1 + theme(
  legend.key.height = unit(.6, "cm"),
  legend.key.width = unit(1, "cm"))

# Get the ggplot Grob
gt = ggplotGrob(g1)

# grid.ls(grid.force(gt))  # To get a list of editable grobs

# Edit the relevant keys
library(grid)
gt <- editGrob(grid.force(gt), gPath("key-1-[3,7]-[1,2]"), 
               grep = TRUE, global = TRUE,
               x0 = unit(0, "npc"), y0 = unit(0.5, "npc"), 
               x1 = unit(1, "npc"), y1 = unit(0.5, "npc")) 

# Draw it
grid.newpage()
g1 <- grid.draw(gt)

m1 <- matrix(rnorm(2000,0,1),2000,1)
dt1 <- data.frame(m1)
names(dt1) <- c("Z")
library(ggplot2)

g2 <- ggplot(dt1, aes(x=Z)) 
g2 <- g2+geom_histogram(aes(y=..density..),      # Histogram with density instead of count on y-axis
                        binwidth=.5,
                        colour="black", fill="white",breaks=seq(-2, 2, by = 0.1)) 
g2 <- g2 + stat_function(fun=dnorm,
                         color="black",geom="area", fill="gray", alpha=0.1,
                         args=list(mean=mean(dt1$Z), 
                                   sd=sd(dt1$Z)))
g2 <- g2+  geom_vline(aes(xintercept=0, linetype="Valor Verdadeiro"),show.legend =TRUE)
g2 <- g2+  geom_vline(aes(xintercept=mean(dt1$Z, na.rm=T),    linetype="Valor Estimado"),show.legend =TRUE)
g2 <- g2+  scale_linetype_manual(values=c("dotdash","solid")) # Overlay with transparent density plot
g2 <- g2+  xlab(expression(paste(gamma[1])))+ylab("Densidade")
g2 <- g2+  theme(plot.margin=unit(c(0.5, 0.5, 0.5, 0.5), units="line"),
                 legend.position = "top",
                 legend.justification = c(0,0),
                 legend.box.just = "top",
                 legend.margin = margin(0,0,-10,-5),
                 legend.title=element_blank(),
                 legend.direction = "horizontal",
                 legend.background = element_rect(fill="transparent", size=.5, linetype="dotted"))
g2 <- g2+ guides(linetype = guide_legend(override.aes = list(size = 1)))


# Adjust key height and width
g2 = g2 + theme(
  legend.key.height = unit(.6, "cm"),
  legend.key.width = unit(1, "cm"))

# Get the ggplot Grob
gt2 = ggplotGrob(g2)

# grid.ls(grid.force(gt))  # To get a list of editable grobs

# Edit the relevant keys
library(grid)
gt2 <- editGrob(grid.force(gt2), gPath("key-1-[3,7]-[1,2]"), 
                grep = TRUE, global = TRUE,
                x0 = unit(0, "npc"), y0 = unit(0.5, "npc"), 
                x1 = unit(1, "npc"), y1 = unit(0.5, "npc")) 

# Draw it
grid.newpage()
g2 <- grid.draw(gt2)

#library(gridExtra)
#grid.arrange(grid.draw(gt),grid.draw(gt2))

Answer 1

这是否符合您的要求grid.arrange(gt, gt2, ncol = 2)？

(g1 和 g2 在你的代码中都是 NULL 因为它们是通过调用 grid.draw 创建的，这不会 return 任何东西)

要使用 facet_wrap，您需要将所有数据放入一个长格式的数据框中：

library(tidyr)
df <- cbind.data.frame(dt, dt1)
df <- gather(df, key = "db", value = "value")

然后剧情：

p <- ggplot(df, aes(x = value)) + 
geom_histogram(aes(y = ..density..), 
               binwidth = .5,
               breaks = seq(-2, 2, by = .1)) +
facet_wrap(~ db)

使用“gridExtra”和多个方面并排绘制图形

Graphs side by side using `gridExtra` and multiple facets

r

dataframe

ggplot2

facet

gridextra