如何在现有R语言绘图代码中正确添加指定垂直线？

你现有R代码可生成包含两个直方图和彩色矩形的绘图，但添加指定的四条垂直线时无法显示，核心原因是变量作用域冲突——x_right、mean_x等变量在base_plot_fun函数内部定义，外部代码无法直接访问；同时后续叠加直方图时重复指定xlim，可能导致绘图范围不一致，进一步影响垂直线显示。

以下提供两种不破坏原有功能的解决方案：

方案一：函数内部添加垂直线（垂直线位于矩形之上、直方图之下） #

修改base_plot_fun函数，在绘制完彩色矩形后直接添加垂直线，利用函数内部的变量作用域直接访问所需参数，同时统一绘图的x范围避免冲突：

V_max=10
V_min=0
x_lim_min=-1
x_lim_max=13

base_plot_fun <- function(x, x1, y, y1) {
  mean_x <- mean(x)
  mean_x1 <- mean(x1)
  
  # 初始化绘图，统一使用外部定义的x范围
  plot(
    0,
    0,
    ylim = range(c(y$counts, y1$counts)),
    xlim = c(x_lim_min, x_lim_max),
    xlab = "x",
    ylab = "counts",
    type = "n"
  )
  
  corners <- par('usr') # 获取绘图区域边界
  
  x_left <- corners[1]
  x_right <- quantile(x, 0.95) # 95分位数，可按需修改
  x1_left <- x_right
  x1_right <- quantile(x1, 0.05) # 5分位数，可按需修改
  
  # 绘制彩色矩形
  rect(
    xleft = x_left,
    xright = x_right,
    ybottom = corners[3],
    ytop = corners[4],
    col = 'lightblue',
    density = 100
  )
  rect(
    xleft = x1_left,
    xright = x1_right,
    ybottom = corners[3],
    ytop = corners[4],
    col = 'pink',
    density = 100
  )
  
  # 直接在函数内添加垂直线
  abline(v=x_right,col='steelblue',lwd=3)
  abline(v=x1_right,col='firebrick',lwd=3)
  abline(v=mean_x,col='forestgreen',lwd=3)
  abline(v=mean_x1,col='forestgreen',lwd=3)
}

x <- rnorm(100,V_min,1)
x1 <- rnorm(100,V_max,1)
y <- hist(x,plot=FALSE,breaks = 20)
y1 <- hist(x1, plot=FALSE,breaks=20)

base_plot_fun(x = x, x1 = x1, y = y, y1 = y1)
plot(y, col='gray48',add=T)
plot(y1, col='gray48',add=T)

方案二：函数返回变量，外部添加垂直线（垂直线位于直方图上层） #

如果需要垂直线显示在直方图上方，可让base_plot_fun返回所需变量，在叠加完直方图后再调用abline：

V_max=10
V_min=0
x_lim_min=-1
x_lim_max=13

base_plot_fun <- function(x, x1, y, y1) {
  mean_x <- mean(x)
  mean_x1 <- mean(x1)
  
  # 初始化绘图，统一使用外部定义的x范围
  plot(
    0,
    0,
    ylim = range(c(y$counts, y1$counts)),
    xlim = c(x_lim_min, x_lim_max),
    xlab = "x",
    ylab = "counts",
    type = "n"
  )
  
  corners <- par('usr') # 获取绘图区域边界
  
  x_left <- corners[1]
  x_right <- quantile(x, 0.95) # 95分位数，可按需修改
  x1_left <- x_right
  x1_right <- quantile(x1, 0.05) # 5分位数，可按需修改
  
  # 绘制彩色矩形
  rect(
    xleft = x_left,
    xright = x_right,
    ybottom = corners[3],
    ytop = corners[4],
    col = 'lightblue',
    density = 100
  )
  rect(
    xleft = x1_left,
    xright = x1_right,
    ybottom = corners[3],
    ytop = corners[4],
    col = 'pink',
    density = 100
  )
  
  # 返回需要的垂直线参数
  return(list(
    x_right = x_right,
    x1_right = x1_right,
    mean_x = mean_x,
    mean_x1 = mean_x1
  ))
}

x <- rnorm(100,V_min,1)
x1 <- rnorm(100,V_max,1)
y <- hist(x,plot=FALSE,breaks = 20)
y1 <- hist(x1, plot=FALSE,breaks=20)

# 调用函数并获取垂直线参数
line_params <- base_plot_fun(x = x, x1 = x1, y = y, y1 = y1)
plot(y, col='gray48',add=T)
plot(y1, col='gray48',add=T)

# 在直方图上层添加垂直线
abline(v=line_params$x_right,col='steelblue',lwd=3)
abline(v=line_params$x1_right,col='firebrick',lwd=3)
abline(v=line_params$mean_x,col='forestgreen',lwd=3)
abline(v=line_params$mean_x1,col='forestgreen',lwd=3)

关键修改说明 #

1. 统一x范围：在初始plot中直接使用外部定义的x_lim_min和x_lim_max，避免后续叠加直方图时重复指定xlim导致绘图范围偏移，确保垂直线处于可视区域内。
2. 变量作用域处理：方案一利用函数内部作用域直接访问参数；方案二通过返回变量实现外部访问，可灵活控制垂直线的绘制层级。

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