Get axis limits from ggplot object [duplicate]

Question

Sometimes I will make two plots that are similair, but have different data. As such, the X- and Y-axis will have different ranges.

library(ggplot2)
library(ggpubr)
#> Loading required package: magrittr
df1 <- data.frame(x=runif(10)*2,y = runif(10)*2)
df2 <- data.frame(x=runif(10)*3,y = runif(10)*1)
p1 <- qplot(x = x, y = y, data = df1, geom = "line")
p2 <- qplot(x = x, y = y, data = df2, geom = "line")

ggarrange(p1,p2)

^{Created on 2020-07-09 by the reprex package (v0.3.0)}

This can manually be overcome by explicitly stating a range with xlim and ylim, but this is both tiresome and could lead to some data being outside of the specified range if one is not careful.

An ideal solution would be to dynamically obtain the limits from the p1, and if these are larger than those of p2, use those instead for p2. For example p2 + xlim(getLimits(p1)).

Is something like this supported?

EDIT: This question was suggested, but the answer applies to an older version of ggplot2. Additionally, the supported p1$coordinates$limits only return the manually specified limits, which defeats the purpose.

Answer 1

ggplot now has layer_*** convenience functions for extracting information from ggplots. In this case you can use the layer_scales function:

layer_scales(p1)$y$get_limits()

[1] 0.1499588 1.9527970

So you could do something like:

library(tidyverse)
library(patchwork)
library(ggpubr)
theme_set(theme_bw())

set.seed(2)
df1 <- data.frame(x=runif(10)*2,y = runif(10)*2)
df2 <- data.frame(x=runif(10)*3,y = runif(10)*1)
p1 <- qplot(x = x, y = y, data = df1, geom = "line")
p2 <- qplot(x = x, y = y, data = df2, geom = "line")

fnc = function(...) {
  
  p = list(...)
  
  yr = map(p, ~layer_scales(.x)$y$get_limits()) %>% 
    unlist %>% range
  
  xr = map(p, ~layer_scales(.x)$x$get_limits()) %>% 
    unlist %>% range
  
  p %>% map(~.x + xlim(xr) + ylim(yr))
}

wrap_plots(fnc(p1, p2))