只需继续使用过度绘制即可bxp
:
set.seed(123)
Mydata = sample(x=100:300, size = 500, replace = T)
Mydata = c(Mydata, 1, 500)
bp <- boxplot(Mydata, range=0, plot=FALSE)
vals <- c(
min=min(Mydata),
quantile(Mydata, c(0.025, 0.25, 0.5, 0.7, 0.75, 0.975)),
max=max(Mydata)
)
bxp(bp, whisklty=0, staplelty=0)
bp$stats[2:4,] <- c(vals[2], Inf, vals[5])
bxp(bp, whisklty=0, staplelty=0, add=TRUE)
bp$stats[2:4,] <- c(vals[2], Inf, vals[7])
bxp(bp, whisklty=1, staplelty=1, add=TRUE)