如何用Base R的gap.boxplot绘制带Y轴间隙的箱线图？解决报错

你遇到的报错是因为gap.boxplot的间隙检查逻辑比较严苛——哪怕某个类别（比如V3）的统计量刚好靠近间隙边界，或者函数内部计算的须线范围不小心触碰到间隙区间，就会触发这个错误。不用纠结这个，直接用Base R手动实现带间隙的箱线图更灵活，或者调整间隙范围适配gap.boxplot。

方法一：Base R手动实现带间隙的箱线图 #

这种方式完全可控，能自由调整间隙位置和样式：

# 准备数据和标签
data_N3_minmax = structure(c(2.71177214885543, 4.36849175181861, 17.626304994383, 
17.7213767231357, 1116.19429181322, 1121.5558117271, NA, NA, 
5.54438655223586, 15.5774104850416, 7.34658260987168, 25.2590667564114, 
8.44869502735733, 17.5713360077048), .Dim = c(2L, 7L)) 
labels <- c("V1", "V2", "V3", "V4", "V5", "V6", "V7")

# 计算每个变量的箱线统计量（跳过全NA的列）
box_stats <- lapply(1:ncol(data_N3_minmax), function(i) {
  x <- data_N3_minmax[,i]
  if(all(is.na(x))) NULL else boxplot.stats(x)$stats
})

# 设置间隙参数
gap_low <- 30    # 间隙下限
gap_high <- 1000 # 间隙上限
bottom_y <- c(0, gap_low)       # 底部Y轴范围
top_y <- c(gap_high, max(unlist(box_stats), na.rm=TRUE) + 10) # 顶部Y轴范围

# 初始化画布，绘制底部区间的箱线
par(mar = c(5, 5, 2, 2)) # 调整边距给间隙留空间
plot(1, type="n", xlim=c(0.5, length(labels)+0.5), ylim=bottom_y, 
     xlab="", ylab="数值", xaxt="n", yaxt="n")
axis(2, at=seq(0, gap_low, 10)) # 底部Y轴刻度
axis(1, at=1:length(labels), labels=labels) # X轴标签

# 绘制除V3外的箱线
for(i in c(1,2,4,5,6,7)){
  if(!is.null(box_stats[[i]])){
    stats <- box_stats[[i]]
    rect(i-0.2, stats[2], i+0.2, stats[4], col="lightblue") # 箱体
    segments(i, stats[3], i, stats[3], lwd=2) # 中位数线
    segments(i, stats[1], i, stats[2]) # 下须线
    segments(i, stats[4], i, stats[5]) # 上须线
    segments(i-0.1, stats[1], i+0.1, stats[1]) # 下端点
    segments(i-0.1, stats[5], i+0.1, stats[5]) # 上端点
  }
}

# 绘制间隙标记（两条斜线）
segments(0.5, gap_low, length(labels)+0.5, gap_low, lty=2)
segments(0.5, gap_low - 1, length(labels)+0.5, gap_low - 3, lty=2)
segments(0.5, gap_low - 3, length(labels)+0.5, gap_low - 1, lty=2)

# 叠加绘制顶部区间的V3箱线
par(new=TRUE)
plot(1, type="n", xlim=c(0.5, length(labels)+0.5), ylim=top_y, 
     xlab="", ylab="", xaxt="n", yaxt="n")
axis(2, at=seq(gap_high, top_y[2], 50)) # 顶部Y轴刻度

# 绘制V3的箱线
i <- 3
stats <- box_stats[[i]]
rect(i-0.2, stats[2], i+0.2, stats[4], col="lightblue")
segments(i, stats[3], i, stats[3], lwd=2)
segments(i, stats[1], i, stats[2])
segments(i, stats[4], i, stats[5])
segments(i-0.1, stats[1], i+0.1, stats[1])
segments(i-0.1, stats[5], i+0.1, stats[5])

方法二：调整间隙范围适配plotrix的gap.boxplot #

报错的核心原因是你设置的间隙区间不小心覆盖了某个箱线的统计量（比如V6的上端须线是25.26，你之前设的间隙下限30其实没问题，但函数可能对边界有额外检查）。把间隙下限调到比所有低数值类别的最大须线稍高，同时把间隙上限调到V3的最小统计量稍低，就能避开冲突：

library(plotrix)
labels <- c("V1", "V2", "V3", "V4", "V5", "V6", "V7")
# 间隙设为26到1110，避开所有箱线统计量
gap.boxplot(data_N3_minmax, gap=list(top=c(26, 1110)), xlab="", ylab="数值", names=labels)

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