ggplot2中的双Y轴用于多面板图形

Question

我们不能轻易地使用ggplot2创建双y轴。

• kohske提供了一种解决方案，用于为单面板图添加双y轴（http://rpubs.com/kohske/dual_axis_in_ggplot2）
• @Andreas展示了一个示例，使用ggplot2对齐两个绘图（仍然是单面板图）（Plot with 2 y axes, one y axis on the left, and another y axis on the right）

如果图形具有多个面板（由geom_wrap或geom_grid创建），则这些技巧无法使用，例如下面的4面板图。我想将p1的轴添加到合并图的面板右侧。

编辑：我在此处删除了我的代码，并将自己的解决方案发布为答案。以防有人需要它。

- Bangyou

---

我尝试过这个，但对我没有用。 - Bangyou

要重新打开这个问题，您应该提供一个解释，为什么它对您没有起作用。 - mhlester

抱歉，我在手机上点击了重新打开按钮，但是失去了连接。我会尽快添加我的解释。 - Bangyou

我已经对我的问题进行了解释，并且取得了一些获取轴grobs的进展。然而，当调用gtable_add_grob时，我遇到了一个错误。 - Bangyou

@Bangyou 只是为了明确起见，你可以回答自己的问题，这也是推荐的做法。 - MattLBeck

Answer 1

我发布了自己的解决方案作为答案，以防有人需要它。

library(ggplot2)
library(gtable)
library(grid)
df1 <- expand.grid(list(x = 1:10, z = 1:4, col = 1:2))
df2 <- df1
df1$y <- df1$x * df1$x + df1$col * 10
df2$y <- 4 * df2$x + df1$col * 10

p1 <- ggplot(df1)
p1 <- p1 + geom_point(aes(x, y, colour = factor(col)))
p1 <- p1 + facet_wrap(~z)
p1 <- p1 + ylab('Y label for figure 1')
p1 <- p1 + theme(legend.position = 'bottom',
    plot.margin = unit(c(1, 2, 0.5, 0.5), 'cm'))

p2 <- ggplot(df2)
p2 <- p2 + geom_line(aes(x, y, colour = factor(col)))
p2 <- p2 + facet_wrap(~z)
p2 <- p2 + ylab('Y label for figure 2')
p2 <- p2 + theme(legend.position = 'bottom',
    plot.margin = unit(c(1, 2, 0.5, 0.5), 'cm'))

g1 <- ggplot_gtable(ggplot_build(p1))
g2 <- ggplot_gtable(ggplot_build(p2))

combo_grob <- g2
pos <- length(combo_grob) - 1
combo_grob$grobs[[pos]] <- cbind(g1$grobs[[pos]],
    g2$grobs[[pos]], size = 'first')

panel_num <- length(unique(df1$z))
for (i in seq(panel_num))
{
    # grid.ls(g1$grobs[[i + 1]])
    panel_grob <- getGrob(g1$grobs[[i + 1]], 'geom_point.points',
        grep = TRUE, global = TRUE)
    combo_grob$grobs[[i + 1]] <- addGrob(combo_grob$grobs[[i + 1]], 
        panel_grob)
}       


pos_a <- grep('axis_l', names(g1$grobs))
axis <- g1$grobs[pos_a]
for (i in seq(along = axis))
{
    if (i %in% c(2, 4))
    {
        pp <- c(subset(g1$layout, name == paste0('panel-', i), se = t:r))

        ax <- axis[[1]]$children[[2]]
        ax$widths <- rev(ax$widths)
        ax$grobs <- rev(ax$grobs)
        ax$grobs[[1]]$x <- ax$grobs[[1]]$x - unit(1, "npc") + unit(0.5, "cm")
        ax$grobs[[2]]$x <- ax$grobs[[2]]$x - unit(1, "npc") + unit(0.8, "cm")
        combo_grob <- gtable_add_cols(combo_grob, g2$widths[g2$layout[pos_a[i],]$l], length(combo_grob$widths) - 1)
        combo_grob <- gtable_add_grob(combo_grob, ax,  pp$t, length(combo_grob$widths) - 1, pp$b)
    }
}

pp <- c(subset(g1$layout, name == 'ylab', se = t:r))

ia <- which(g1$layout$name == "ylab")
ga <- g1$grobs[[ia]]
ga$rot <- 270
ga$x <- ga$x - unit(1, "npc") + unit(1.5, "cm")

combo_grob <- gtable_add_cols(combo_grob, g2$widths[g2$layout[ia,]$l], length(combo_grob$widths) - 1)
combo_grob <- gtable_add_grob(combo_grob, ga, pp$t, length(combo_grob$widths) - 1, pp$b)
combo_grob$layout$clip <- "off"

grid.draw(combo_grob)

P1：

P2：

合并后：

- Bangyou

---

不确定在 ggplot2 v2.0.0 或 grid 3.1.2 或 gtable 0.1.2 中是否有所改变，但我收到一个错误，说“Error in addGrob(combo_grob$grobs[[i + 1]], panel_grob) : it is only valid to add a 'grob' to a "gTree"”，你有什么想法吗？我正在使用 geom_bar 来绘制我的 p1。 - Dan