在箱线图中部分更改大纲样式

Question

假设我有以下数据集

data <- c(
9.5, 27.9, 7.8, 17.8, 31.4, 25.9, 27.4, 
25.2, 31.1, 34.7, 42, 29.1, 32.5, 30.3, 33, 33.8, 41.1, 34.5, 62)

当我在r中绘制箱形图时

boxplot(data)

我得到三个离群值7.8、9.5和62，在图中用三个小圆圈表示。

在这里，我想将最大异常值（即62）的pch更改为实心圆，但不要将其他两个较小的轮廓线更改为pch。

以下是我尝试过的方法，但是不起作用：

boxplot(data, outpch = ifelse(data >= 60, 16, 1))

有没有办法做到这一点？ 谢谢

Answer 1

我认为您无法直接在boxplot函数中执行此操作，因为outpch中的boxplot参数不需要向量，但是我们可以使用points函数以不同的方式显示离群值。

bp <- boxplot(data, outpch = NA) 
with(bp, points(group, out, pch = ifelse(out >=60, 16, 1)))