R 代码用于绘制多个，将图紧密对齐？

Question

我正在尝试将 plot 多个条形图和散点图紧密地结合在一起（即没有间隙将图分开）。 我正在绘制两条不同河流的流量和盐度值，尽管其中一条有一年前/一年后的数据，但重要的是要使年份保持一致。 我希望没有间隙、轴刻度标签、轴标签、绘图之间的空格，尽管每个 plot 的刻度线将有助于使其可读。 我还希望整个过程中有非常浅灰色的垂直线，以便于比较给定年份两条河流的不同值。

我已经尝试在 R 基本计划和 ggplot2 中执行此操作，但无法成功。 我找不到一种简单的方法来对齐 x 轴（它们不会从点到条形图以及从河流到河流对齐），也无法让这些图垂直相邻。

我想将 output 导出为 multi-panel.png 图像。

这是示例数据以及到目前为止我在 ggplot 中的尝试方式：

river1 = data.frame(year     =  1991:2017, 
                    flow     =  rnorm(n=27, mean = 900, sd = 350),
                    salinity =  rnorm(n=27, mean = 3000, sd = 800))
river2 = data.frame(year     =  1992:2019, 
                    flow     =  rnorm(n=28, mean = 20000, sd = 2400),
                    salinity =  rnorm(n=28, mean = 700, sd = 1500))

library(ggplot2)
(r1sal = ggplot(river1, aes(x = year, y = salinity)) + 
    geom_point(stat = "identity", color="darkred", fill="darkred", size = 2) + 
    labs(title = "River 1 flow and flow-weighted salinity", x = "", y = "Streamflow (ML)") +
    scale_x_continuous(breaks= seq(from=1991, to= 2019, by=2)) + theme_classic())

(r1flow = ggplot(river1, aes(x = year, y = flow)) + 
    geom_bar(stat = "identity", color="blue", fill="lightblue") + 
    labs(x = "Water year", y = "Streamflow (ML)") +
    scale_x_continuous(breaks= seq(from=1991, to= 2019, by=2)) + theme_classic())

(r2sal = ggplot(river2, aes(x = year, y = salinity)) + 
    geom_point(stat = "identity", color="darkred", fill="darkred", size = 2) + 
    labs(title = "River 2 flow and flow-weighted salinity", x = "Water year", y = "Streamflow (ML)") +
    scale_x_continuous(breaks= seq(from=1991, to= 2019, by=2)) + theme_classic())

(r2flow = ggplot(river2, aes(x=year, y=flow)) + 
       geom_bar(stat = "identity", color="blue", fill="lightblue") + 
       labs( x = "Water year", y = "Streamflow (ML)") +
       scale_x_continuous(breaks= seq(from=1991, to= 2019, by=2)) + theme_classic())

library(ggpubr)
ggarrange(r1sal, r1flow, r2sal, r2flow,
            labels = c("A", "B", "C","D"),
            ncol = 1, nrow = 4)

我知道必须有一个简单的方法来做到这一点。 非常感谢您的帮助！

Answer 1

我会考虑使用方面。

river1 = data.frame(year     =  1991:2017, 
                    flow     =  rnorm(n=27, mean = 900, sd = 350),
                    salinity =  rnorm(n=27, mean = 3000, sd = 800))
river2 = data.frame(year     =  1992:2019, 
                    flow     =  rnorm(n=28, mean = 20000, sd = 2400),
                    salinity =  rnorm(n=28, mean = 700, sd = 1500))


library(tidyverse)

rivers <- bind_rows(river1 = river1, river2 = river2, .id = "River") %>% 
  pivot_longer(cols = c("flow", "salinity"), names_to = "variable", values_to = "value") 

ggplot(data = filter(rivers, variable == "salinity" ), mapping = aes(x = year, y = value)) +
  geom_col(data = filter(rivers, variable == "flow" ), fill = "blue") +
  geom_point(color="darkred") +
  geom_line(color="darkred") +
  facet_grid(River + variable ~ ., scales = "free_y")

^{由reprex package (v0.3.0) 于 2020 年 5 月 11 日创建}

另一种方法是使用主题将 axis.text 设置为 element_blank() 然后组合数字（我会用patchwork包来做到这一点）