Cowplot plot_grid (): Titre du plot

Aug 26 2020

J'essaye d'inclure un titre dans un graphique assemblé construit en utilisant plot_grid(). J'ai essayé à la fois la solution proposée dans le cowplotguide de référence ici par @Claus Wilke, et la solution de contournement suggérée ici dans le cadre des questions elles-mêmes (extraire le titre d'un ggplotobjet vide ).

Dans les deux cas, je suis capable de construire un objet avec le titre seul. Cependant, lorsque j'essaie de l'inclure dans un graphique à une colonne avec deux autres objets, quelque chose d'étrange se produit. Si j'essaie de l'inclure en haut (par exemple, en spécifiant un plot_grid()paramètre rel_heights(0.1,1,0.5), comme le suggère le premier lien que j'ai publié ci-dessus, la sortie n'inclut qu'un espace vide. Inversement, si j'essaie d'inverser l'ordre des graphiques (c'est-à-dire en incluant le titre ci-dessous) les graphiques principaux), le résultat est exactement comme il se doit et le titre est tracé sous les graphiques.

Une idée pour pourquoi cela se produit-il? ma configuration n'est pas remarquablement différente des exemples faits dans les références que j'ai incluses, sauf pour l'orientation du graphique principal (qui est quadrillé comme une seule colonne plutôt que sur une seule ligne)

MON CODE


# 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))

Une idée pour pourquoi cela se produit-il? J'ai essayé de construire le graphique en construisant le graphique en colonnes avant et en ajoutant le titre dans une deuxième étape, mais sans succès. Ci-dessous, je fournis un sous-échantillon des données.

Cowplot plot_grid (): Titre du plot

Réponses