Início / coding / Perguntas / 79529355
Accepted
Mohan
Asked: 2025-03-24 01:26:36 +0800 CST

ggplot: como colocar dois gráficos lado a lado e expandir a tela para manter os layouts padrão?

  • 772

Tenho alguns ggplots que ajustei para garantir que as anotações se encaixem bem dentro dos limites do gráfico. Por exemplo:

insira a descrição da imagem aqui

Agora eu quero colocar duas dessas imagens lado a lado. Eu tentei fazer isso com grid.arrange e patchwork, mas o layout do gráfico muda, presumivelmente porque ele está sendo encaixado em um espaço menor:

insira a descrição da imagem aqui

Tudo fica mais apertado, os pontos ficam mais próximos, o que dificulta a leitura, e as caixas de anotação ficam fora do gráfico.

O que eu quero fazer é expandir a tela exatamente o suficiente para que cada gráfico pareça exatamente como era quando plotado sozinho (ou seja, quando plotado usando o tamanho de tela padrão). Presumivelmente, precisaria ser um pouco mais do que o dobro da largura, para permitir o preenchimento. Como eu conseguiria isso?

Código de exemplo:

library(ggplot2)
library(ggtext)
library(tidyverse)
library(patchwork)

N <- 10000
rho = 0.5
  
mu <- c(0, 0)
sigma <- matrix(c(1,   rho,
                  rho, 1   ), 
                nrow = 2, ncol = 2, byrow = TRUE)
            
vs = MASS::mvrnorm(N, mu, sigma)
df = tibble(rMZall = vs[,1], rDZall = vs[,2], n_twin_pairs = runif(N, 5000, 50000))

plot <- function(min_twin_pairs) {
  ggplot(data=df |> filter(n_twin_pairs > min_twin_pairs), aes(rMZall, rDZall)) + 
    ggtitle(sprintf("Studies with > %d twin pairs", min_twin_pairs)) + 
    theme(legend.position="none") +
    # Bounds
    scale_x_continuous(name = "rMZ", limits = c(-0.05, 1.1), breaks = seq(0,1,0.25)) +
    scale_y_continuous(name = "rDZ", limits = c(-0.05, 1), breaks = seq(0,1,0.25)) +
    coord_fixed() +
    # Scatterplot
    stat_density2d(aes(fill=..level.., alpha=..level..), geom='polygon', colour='black', linewidth = 0.1) + 
    scale_fill_continuous(low="green", high="red") +
    geom_point(aes(size = sqrt(n_twin_pairs))) +
    scale_size_continuous(range = c(0, 1)) + 
    # AE and AD lines
    geom_abline(intercept = 0, slope = .5) +
    annotate(geom = "richtext", x = 1, y = .5, label = "AE model", angle = 26.6) + 
    geom_abline(intercept = 0, slope = .25) +
    annotate(geom = "richtext", x = 1, y = .25, label = "AD model", angle = 14.0) 
}

plot1000 = plot(min_twin_pairs = 1000)
plot4000 = plot(min_twin_pairs = 4000)

plot1000
plot4000

plot1000 + plot4000

1 respostas

  1. Best Answer

    Você poderia usar geomtextpath , especificamente geomtextpath::geom_labelabline(). Isso lida com o posicionamento automaticamente, então você não precisa se preocupar com ângulos, dimensionamento de tela e coisas do tipo.

    library(tidyverse)
library(geomtextpath)
library(patchwork)
set.seed(13)

plot <- function(min_twin_pairs) {
  ggplot(data=df |> filter(n_twin_pairs > min_twin_pairs), aes(rMZall, rDZall)) + 
    ggtitle(sprintf("Studies with > %d twin pairs", min_twin_pairs)) + 
    theme(legend.position="none") +
    # Bounds
    scale_x_continuous(name = "rMZ", limits = c(-0.05, 1.1), breaks = seq(0,1,0.25)) +
    scale_y_continuous(name = "rDZ", limits = c(-0.05, 1), breaks = seq(0,1,0.25)) +
    coord_fixed() +
    # Scatterplot
    stat_density2d(aes(fill=..level.., alpha=..level..), geom='polygon', colour='black', linewidth = 0.1) + 
    scale_fill_continuous(low="green", high="red") +
    geom_point(aes(size = sqrt(n_twin_pairs))) +
    scale_size_continuous(range = c(0, 1)) + 
    # AE and AD lines
    geom_labelabline(intercept = 0, slope = .5, label = "AE model", hjust = 0.85) +
    geom_labelabline(intercept = 0, slope = .25, label = "AD model", hjust = 0.85)
}

plot1000 <- plot(min_twin_pairs = 1000)
plot4000 <- plot(min_twin_pairs = 4000)

plot1000 + plot4000

    • 3

relate perguntas