使用cowplot实现密度图与散点图的精准对齐问题求助

我仔细看了你的双变量绘图函数代码，发现右侧密度图和散点图轴不对齐的问题，主要出在轴范围同步的时机、翻转后的轴对应关系，还有plot_grid的对齐逻辑上。咱们来一步步修正这个问题：

问题根源拆解 #

1. 轴范围设置时机错误：你先创建了子图，再设置xlim/ylim，这会导致后续调整主题和边距时，轴的实际渲染基准已经错位。
2. 翻转后的轴对应关系混乱：右侧密度图用coord_flip()后，原来的x轴变成了y轴，但你没有确保这个翻转后的轴和散点图的y轴范围完全匹配。
3. plot_grid对齐参数使用不精准：虽然最后用了align="hv"，但没有在子组合阶段就做垂直/水平对齐，导致整体布局错位。

修正后的完整函数 #

bivariate_plot <- function(df, var1, var2, density = T, box = F) {
  require(ggplot2)
  require(cowplot)
  
  # 提前计算统一的轴范围，这是对齐的核心基础
  min_exp <- min(df[[var1]], df[[var2]]) - 0.01
  max_exp <- max(df[[var1]], df[[var2]]) + 0.01
  density_ylim <- c(0, 2) # 密度图的y轴范围单独定义，方便后续调整
  
  # 1. 散点图：先设置轴范围，再调整主题
  scatter <- ggplot(df, aes(.data[[var1]], .data[[var2]], color = "red")) + 
    geom_point(alpha = .8) +
    xlim(min_exp, max_exp) +
    ylim(min_exp, max_exp) +
    theme(
      legend.justification = c(0, 1), 
      legend.position = c(0, 1),
      plot.margin = unit(c(0, 0, 0.5, 0.5), "cm") # 底部留空对接下方密度图
    )
  
  # 2. 底部密度图（对应var1）：隐藏冗余轴元素，设置统一范围
  plot1 <- ggplot(df, aes(.data[[var1]], fill = "red")) + 
    geom_density(alpha = .5) +
    ylab("G1 density") +
    xlim(min_exp, max_exp) +
    ylim(density_ylim) +
    theme(
      legend.position = "none",
      axis.title.x = element_blank(), 
      axis.text.x = element_blank(), 
      axis.ticks.x = element_blank(),
      plot.margin = unit(c(0.5, 0, 0, 0.7), "cm") # 顶部留空对接散点图底部
    )
  
  # 3. 右侧密度图（对应var2）：先设范围再翻转，确保轴对应关系正确
  plot2 <- ggplot(df, aes(.data[[var2]], fill = "red")) + 
    geom_density(alpha = .5) +
    ylab("G2 density") +
    xlim(min_exp, max_exp) + # 翻转前的x轴对应散点图的y轴，用统一范围
    ylim(density_ylim) +
    coord_flip() + # 翻转后，原x轴变为y轴，自动匹配散点图y轴范围
    theme(
      legend.position = "none",
      axis.title.y = element_blank(), 
      axis.text.y = element_blank(), 
      axis.ticks.y = element_blank(),
      plot.margin = unit(c(0, 0.5, 0.5, 0), "cm") # 左侧留空对接散点图右侧
    )
  
  # 组合左侧列：垂直对齐顶部密度图和散点图
  left_col <- plot_grid(plot1, scatter, ncol = 1, rel_heights = c(1, 3), align = "v")
  # 组合右侧列：和左侧列高度比例一致，垂直对齐空占位和右侧密度图
  right_col <- plot_grid(NULL, plot2, ncol = 1, rel_heights = c(1, 3), align = "v")
  
  # 最终组合左右列：水平对齐，确保所有轴完美贴合
  perfect <- plot_grid(left_col, right_col, ncol = 2, rel_widths = c(3, 1), align = "h", axis = "lrbl")
  print(perfect)
}

关键修改说明 #

• 替换变量引用方式：用.data[[var1]]替代eval(parse(text=var1))，这是ggplot2官方推荐的安全引用变量的方法，避免字符串解析带来的潜在问题。
• 提前统一轴范围：在创建任何子图前就计算好所有轴的统一min/max，确保散点图和两个密度图的相关轴范围完全一致，这是对齐的核心前提。
• 调整右侧密度图创建顺序：先设置xlim（对应翻转后的y轴）再执行coord_flip()，保证翻转后的轴和散点图的y轴范围完全匹配。
• 精准设置plot_grid对齐参数：在子组合阶段就做垂直对齐，最后做水平对齐，确保所有子图的轴基准完全贴合。
• 简化主题调整：只隐藏必要的轴文本和刻度，保留轴的基本结构，避免破坏对齐基准。

测试验证 #

用你的示例数据测试：

# 生成示例数据
g1 <- c(rnorm(200, mean=350, sd=100), rnorm(200, mean=700, sd=100))
g2 <- c(rnorm(200, mean=350, sd=100), rnorm(200, mean=500, sd=100))
df_exp <- data.frame(var1=log2(g1 + 1), var2=log2(g2 + 1))

# 调用修正后的函数
bivariate_plot(df = df_exp, var1 = "var1", var2 = "var2")

修改后的函数不管切换什么数据，右侧密度图都会和散点图的轴完美对齐，因为所有轴的范围都是统一的，布局组合逻辑也精准匹配了各个子图的位置关系。

内容的提问来源于stack exchange，提问作者Seymour