如何用R的par与pdf绘制多组柱状图，避免截断并让标题居中于柱子

针对你遇到的两个绘图问题，我整理了具体的原因分析和可落地的解决方案，还有修改后的示例代码：

问题1：让标题仅在柱子上方居中 #

默认的main参数是相对于整个绘图面板居中，而非柱子所在的区域——这是因为你手动固定了xlim=c(0,1)，导致绘图区的x范围远大于柱子实际占据的空间。想要让标题精准对齐柱子区域，最优方案是放弃main参数，手动计算柱子的中心位置，用text()函数添加标题：

• 调用barplot()时，保存返回的柱子中心x坐标；
• 根据柱子宽度，计算整个柱子组的左右边界（最小x减去半个柱子宽度，最大x加上半个柱子宽度），取中间值作为标题的x坐标；
• 用text()在该位置添加标题，y坐标设置为柱子最大高度加上一点偏移量，避免和柱子重叠。

问题2：避免柱子截断且保持宽度一致 #

柱子截断的核心原因是固定了xlim=c(0,1)，当柱子数量增多时，柱子总宽度（柱子宽度总和+柱子间隙）会超过xlim的上限。要同时满足柱子宽度统一和不截断的要求，需要动态调整x范围：

• 提前固定width参数（比如你设置的0.2）和柱子间隙space参数（比如设为0.1，确保不同图表的间隙一致）；
• 针对每个图表，根据柱子数量计算总宽度：总宽度 = 柱子数量*width + (柱子数量-1)*space；
• 把xlim设置为c(0, 总宽度)，这样每个面板的x范围刚好容纳所有柱子，既不会截断，又能保持柱子宽度统一。

修改后的示例代码 #

# 创建数据
a <- c(1:100)
b <- c(1:200)
c <- c(1:300)
d <- c(1:400)
e <- c(1:500)
f <- c(1:600)

# 构建数据框
test <- as.data.frame(cbind(a,b))
test1 <- as.data.frame(cbind(a,b,c))
test2 <- as.data.frame(cbind(a,b,c,d))
test3 <- as.data.frame(cbind(a,b,c,d,e))
test4 <- as.data.frame(cbind(a,b,c,d,e,f))

# 计算列均值
a0 <- mean(a)  # 修正原代码中mean(a)的调用，确保是单个数值
a1 <- colMeans(test)
a2 <- colMeans(test1)
a3 <- colMeans(test2)
a4 <- colMeans(test3)
a5 <- colMeans(test4)

# 统一设置柱子宽度和间隙
bar_width <- 0.2
bar_space <- 0.1
# 统一y轴范围，让图表对比更直观
y_max <- max(a0,a1,a2,a3,a4,a5) * 1.1

# 绘图
pdf(file= "/Users/Highf_000/Desktop/prac1.pdf")
par(mfrow = c(2, 3), 
    oma = c(2, 2, 0, 0), 
    mar = c(5, 5, 2, 1)+0.1, 
    mgp = c(2, 1, 0), 
    xpd = NA)

# 第一个图：1根柱子
n_bars <- length(a0)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a0, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

# 第二个图：2根柱子
n_bars <- length(a1)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a1, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

# 第三个图：3根柱子
n_bars <- length(a2)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a2, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

# 第四个图：4根柱子
n_bars <- length(a3)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a3, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

# 第五个图：5根柱子
n_bars <- length(a4)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a4, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

# 第六个图：6根柱子
n_bars <- length(a5)
total_width <- n_bars*bar_width + (n_bars-1)*bar_space
x_coords <- barplot(a5, width = bar_width, space = bar_space, 
                    xlim = c(0, total_width), ylim = c(0, y_max))
title_x <- mean(range(x_coords) + c(-bar_width/2, bar_width/2))
text(title_x, y_max*0.95, "Learning", pos=3, font=2)

dev.off()

额外优化建议 #

如果觉得重复代码太多，可以封装一个小函数来简化绘图流程：

plot_custom_bar <- function(data, title, width=0.2, space=0.1, y_max) {
  n_bars <- length(data)
  total_width <- n_bars*width + (n_bars-1)*space
  x_coords <- barplot(data, width=width, space=space, 
                      xlim=c(0, total_width), ylim=c(0, y_max))
  title_x <- mean(range(x_coords) + c(-width/2, width/2))
  text(title_x, y_max*0.95, title, pos=3, font=2)
}

之后调用这个函数绘图会简洁很多：

# 调用示例
plot_custom_bar(a0, "Learning", bar_width, bar_space, y_max)
plot_custom_bar(a1, "Learning", bar_width, bar_space, y_max)
# ... 以此类推

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