r - fiddle 图 : How is the adjacent value range determined, 为什么它与箱线图不同？

Question

理论上 vioplot 包的 violinplot 是一个 boxplot + 密度函数。

在“箱线图部分”中，

黑框对应于 IQR(确实，请参见下文)和

中线应该对应于相同的范围(相邻值，默认为 1.5 IQR)，但事实并非如此(见下文)。任何人都可以解释为什么它们不同？

require("vioplot")
a = rnorm(100)
range (a)
a = c(a,2,8,2.9,3,4, -3, -5) # add some outliers

par ( mfrow = c(1,2))
boxplot(a, range=1.5)
vioplot(a, range=1.5 )

受到以上好评:



Hintze, J. L. and R. D. Nelson (1998). Violin plots: a box plot-density trace synergism. The American Statistician, 52(2):181-4.

Answer 1

让我用一个简单的例子来说明这一点:

b <- c(1:10, 20)

par(mfrow = c(1,2))
boxplot(b, range=1.5)
vioplot(b, range=1.5 )

R 的箱线图的定义是(借用主题上的 ggplot's help):

The upper whisker extends from the hinge to the highest value that is within 1.5 * IQR of the hinge, where IQR is the inter-quartile range, or distance between the first and third quartiles.

浏览 source code vioplot，我们看到 upper[i] <- min(q3[i] + range*iqd, data.max) .

因此，让我们尝试重现上部晶须值:

# vioplot draws
quantile(b, 0.75) + 1.5 * IQR(b)
# 16

# boxplot draws
max(b[b <= quantile(b, 0.75) + 1.5 * IQR(b)])
# 10