Cowplot plot_grid (): Título en la trama

Aug 26 2020

Estoy tratando de incluir un título en un gráfico ensamblado creado con plot_grid(). Probé tanto la solución propuesta en la cowplotguía de referencia aquí por @Claus Wilke, como la solución sugerida aquí como parte de las preguntas en sí (extrayendo el título de un ggplotobjeto vacío ).

En ambos casos, puedo construir un objeto solo con el título. Sin embargo, cuando trato de incluirlo en un gráfico de una columna junto con otros dos objetos, sucede algo extraño. Si trato de incluirlo en la parte superior (por ejemplo, especificando el plot_grid()parámetro rel_heights(0.1,1,0.5), de manera similar a lo que sugiere el primer enlace que publiqué arriba, la salida incluye solo un espacio vacío. Por el contrario, si trato de invertir el orden de los gráficos (es decir, incluyendo el título a continuación) las parcelas principales), la salida es exactamente como debería ser y el título se traza debajo de las parcelas.

¿Alguna idea de por qué está sucediendo esto? mi configuración no es notablemente diferente de los ejemplos hechos en las referencias que incluí, excepto por la orientación del gráfico principal (que está cuadriculado como una sola columna en lugar de en una sola fila)

MI CÓDIGO


# Building the main plots

# ----------
# Upper plot
# ----------

plot.new() ## clean up device

par(mar = c(2, 4, 2, 4), 
    mfrow=c(1,1))

plot(x = range(dato_reg$date[dato_reg$reg == 'ECA']), y = range(dato_reg$sector_any_food_l[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans']), 
     type ="n", 
     ylab = "Number of Measures",
     #main = "Western Europe and Balcans", sub = "Food", xlab = "",
     #col = "blue", 
     bty = "n",
     xaxs="i",
     yaxs="i",
     xaxt = "n",
     cex.axis = 0.9,
     cex.lab = 0.8,
     cex.main = 0.9,
     fg = 'grey')
lines(x = dato_reg$date[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], y = dato_reg$sector_any_food_l[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], col = "darkgreen")
lines(x = dato_reg$date[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], y = dato_reg$sector_any_food_r[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], col = "maroon")

# Secondary axis
par(new = TRUE, mar = c(2, 4, 2, 4), 
    mfrow=c(1,1))
plot(x = dato_reg$date[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], y = dato_reg$Tot_cov[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], 
     type = "l", 
     xaxt = "n", yaxt = "n", # Set null superposition of axis on LHS
     ylab = "", xlab = "", 
     col = "black", lty = 2, lwd = 1.5,
     xaxs="i", yaxs="i", bty = "n", fg = 'grey')

markers <- c(0,400000,800000, 1200000)

axis(side = 4, at = markers, labels=format(markers, scientific=FALSE), fg = 'grey', cex.axis = 0.9, cex.lab = 0.8)
mtext("Total Number of Covid Cases", side = 4, line = 3, srt = -180, cex = 0.8)

legend("topleft", c("# Policies - Lib", "# Policies - Res", "Total Cases"),
       col = c("darkgreen", "maroon", 'black'), 
       lty = c(1, 1, 2),
       bty = "n", 
       cex = 0.7)

up_plot <- recordPlot()

# -----------
# Bottom plot
# -----------

plot.new() ## clean up device

par(mar = c(2, 4, 1, 4))
plot(x = dato_reg$date[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], 
     y = dato_reg$delta_cases_week_sh[dato_reg$sub_reg == 'ECA: W. Eur. and Balcans'], 
     #main = 'Weekly variation of new cases',
     type = 'l',
     col = 'darkgrey',
     cex = 0.5,
     #xaxt = "n", yaxt = "n", # Set null superposition of axis on LHS
     ylab = "", xlab = "", 
     lty = 5, lwd = 1.2,
     xaxs="i", yaxs="i", bty = "n",
     cex.axis = 0.9,
     cex.lab = 0.8,
     cex.main = 0.9,
     fg = 'grey')

bot_plot <- recordPlot()

# ------------------------
# Building the title object
# ------------------------

# First method

tit <- ggdraw() + 
        draw_label(
                "Western Europe and Balcans",
                #fontface = 'bold',
                #fontfamily = 'Open Sans',
                x = 0,
                hjust = 0#,
                #vjust = 1
        ) +
    theme_delabj() #+

# Second Method
tit2 <- ggplot() + 
        labs(title = "Western Europe and Balcans") + 
        theme_delabj()
tit2 <- ggpubr::as_ggplot(get_title(tit2))

# Does not work
plot_grid(tit, up_plot, bot_plot,
          ncol = 1,
          rel_heights = c(0.1,1,0.5))

# Works
plot_grid(up_plot, bot_plot, tit
          ncol = 1,
          rel_heights = c(1,0.5,0.1))

¿Alguna idea de por qué está sucediendo esto? Intenté construir el gráfico construyendo el diagrama de columna antes y agregando el título en un segundo paso, pero sin éxito. A continuación, proporciono una submuestra de los datos.

Cowplot plot_grid (): Título en la trama

Respuestas