首先，boxplot() 其实会回传一个物件，我叫 z 物件好了。
这个 z 物件本身是一个 list，里头就是 boxplot 的各种资讯。
像是 z$stat 就是各变量的 quantiles，z$name 是各变量的名称。
你可以观察看看 bx.p <- boxplot(split(rt(100, 4), gl(5, 20))) 的
bx.p 实际上是回传了什么给你。
这个 z 物件可以自己建，就像建一个 list 一样。
然后再以 bxp(z物件) 画出一个 boxplot。
请看以下按你的需求设计的例子。
# 二组虚假资料
a <- c(1:100)
b <- c(50:1000)
# 建立 z 物件（本身是一个 list）
bxpObj <- list()
# 一定要有 bxpObj$stats 这项，里头是一个 matrix，
# 其中栏是变量，列是各 quantile，要有5列。
# 各 quantile 分别会被画成下底线、下盒线、中位线、上盒线和上底线。
bxpObj$stats <- rbind(
cbind(
quantile(a, c(0.1, 0.25, 0.5, 0.75, 0.9)),
quantile(b, c(0.1, 0.25, 0.5, 0.75, 0.9))
)
)
# bxpObj$name 这项就是变量名称
bxpObj$name <- c("a", "b")
# 重点来了：用 bxp() 画
bxp(bxpObj)
# 直接用 points() 补画你要求的 20% 和 80% quantile。
points(1, quantile(a, 0.2), pch=1) # a的20% quantile
points(1, quantile(a, 0.8), pch=2) # a的80% quantile
points(2, quantile(b, 0.2), pch=3) # b的20% quantile
points(2, quantile(b, 0.8), pch=4) # b的80% quantile
最后，如果你想自行改写成一个新的 function，也要视你的资料而定。
例如是 myboxplot(变量，组别) 还是 myboxplot(matrix) 的写法会不太一样。
但重点就是我前面写的：建一个 z 物件再喂 bxp() 吃下去就对了。
※ 引述《ming790925 (MING)》之铭言：
: [问题类型]:
: 程式咨询
: [软件熟悉度]:
: 使用者(已经有用R 做过不少作品)
: [问题叙述]:
: 我想画一个 boxplot 但不要内建那种，
: 因为我在意的是某几格特定的百分位数，
: 还有我不希望画出outliers
: 因为找了满久的 都没有找到能达到这样的套件
: 想请问各位是否有印象有什么套件已经能完成?
: 举例来说:
: 除了中间的盒子外，我想要延伸出去到90% 10% quantile的值
: 并在80% 20%也加上标记
: 感激!!!!!
: [关键字]:
: boxplot quntile