Как мне перекрыть оси фигуры и расширить заголовок одной оси на две оси в ggplot?

Question

У меня есть фигура, которая включает четыре сюжета. По всей оси Y отложена температура, а по двум осям X - Time1 и Time2. Я складываю и выравниваю все четыре графика в фигуру 2x2, используяpatchwork пакет и plot_layout()функция в R. Поскольку все оси Y являются одной и той же переменной, я хочу показать только один заголовок оси Y в левой части графика, но я не могу понять, как это сделать. Я также хочу скрыть оси Y на двух графиках в правом столбце, но когда я перемещаю фигуры, правый столбец перекрывает левый столбец (я хочу противоположное).

Вот пример того, что я сделал до сих пор.

library(ggplot2)
library(patchwork)

df1 <- data.frame(x1 = rnorm(20), y1 = rnorm(20))
df2 <- data.frame(x2 = rnorm(20), y2 = rnorm(20))
df3 <- data.frame(x3 = rnorm(20), y3 = rnorm(20))
df4 <- data.frame(x4 = rnorm(20), y4 = rnorm(20))

p1 <- ggplot(data = df1, aes(x = x1, y = y1)) +
  geom_point() +
  ylab("Temperature") +
  theme_bw() +
  theme(panel.grid = element_blank(),
           plot.margin = unit(c(0,0,0,0),"cm"))

p2 <- ggplot(data = df2, aes(x = x2, y = y2)) +
  geom_point() +
  ylab("Temperature") +
  theme_bw() +
  theme(panel.grid = element_blank(),
        plot.margin = unit(c(0,0,0,0),"cm"))

p3 <- ggplot(data = df3, aes(x = x3, y = y3)) +
  geom_point() +
  xlab("Time1") +
  ylab("Temperature") +
  theme_bw() +
  theme(panel.grid = element_blank(),
        plot.margin = unit(c(-0.5,0,0,0),"cm"))

p4 <- ggplot(data = df4, aes(x = x4, y = y4)) +
  geom_point() +
  xlab("Time2") +
  ylab("Temperature") +
  theme_bw() +
  theme(panel.grid = element_blank(),
        plot.margin = unit(c(-0.5,0,0,0),"cm"))

p5 <- p1/p2/p3/p4 + plot_layout(ncol = 2, heights = c(1,1))
p5

Что дает эту цифру:

Однако мне бы хотелось, чтобы фигура выглядела примерно так:

Answer 1

Как насчет использования facet_grid? Но сначала нам нужно очистить данные:

library(dplyr)
library(tidyr)
libary(ggplot2)
df1 <- data.frame(Sample = 1, Time1 = rnorm(20), y = rnorm(20))
df2 <- data.frame(Sample = 1, Time2 = rnorm(20), y = rnorm(20))
df3 <- data.frame(Sample = 2, Time1 = rnorm(20), y = rnorm(20))
df4 <- data.frame(Sample = 2, Time2 = rnorm(20), y = rnorm(20)) 

data <- bind_rows(mget(ls(pattern = "df[0-9]"))) %>%
  pivot_longer(cols = c(-y,-Sample), 
               names_to = "Time",
               values_to = "x") %>%
  dplyr::filter(complete.cases(.))
data
## A tibble: 80 x 4
#   Sample       y Time       x
#    <dbl>   <dbl> <chr>  <dbl>
# 1      1 -0.449  Time1  0.768
# 2      1 -1.10   Time1  0.410
# 3      1 -0.529  Time1 -1.98 
# 4      1 -0.485  Time1 -1.09 
# 5      1  0.128  Time1  1.06 
## … with 75 more rows

Теперь, когда у нас есть данные в правильной форме, мы можем facet_wrap по времени:

ggplot(data = data, aes(x = x, y = y)) +
  geom_point() +
  ylab("Temperature") +
  xlab(element_blank()) + 
  facet_grid(cols = vars(Time), rows = vars(Sample), scales = "free",switch = "both") + 
  theme_bw() +
  theme(panel.grid = element_blank(),
           plot.margin = unit(c(0,0,0,0),"cm"),
        strip.background = element_blank(),
        strip.placement = "outside", strip.text.y = element_blank())

https://i.sta ck.imgur.com/xT4NE.png

Обратите внимание, что нет возможности иметь свободную ось Y в одной строке или свободную ось X в одном столбце.