我正在用geom_boxplot在ggplot2中制作一个盒子图。然而,我发现胡须的长度是不正确的,我不知道为什么。这是我的数据:
value = c(1.3739117,0.8709891,3.4510461,0.8470309,1.4838725,0.6942611,1.3095816,3.0444649,19.2785424,1.0866242,0.9376845,2.2343836, 20.7975509, 20.3102489, 18.0046679,1.4197519)
data = data.frame(value)
ggplot(data, aes(y = value)) +
stat_boxplot(geom = "errorbar", width = 0.3) +
geom_boxplot(width = 0.5)
我看到了这样的情节:
第3分位数与上晶须重叠。我手动进行了计算,结果如下:
summary(data)
Min. : 0.6943
1st Qu.: 1.0494
Median : 1.4518
Mean : 6.0715
3rd Qu.: 7.0895
Max. :20.7976
根据geom_boxplot的解释:上晶须从铰链延伸到最大值不超过1.5 * IQR (其中IQR是四分位数范围,或第一和第三四分位数之间的距离)。较低的晶须从铰链延伸到铰链最多1.5 * IQR的最小值。
在我的例子中,IQR是: 7.0895-1.0494 = 6.0401。
较低的晶须应为: 0.6943 - 1.5*6.0401 = -8.36585。
上晶须应为: 7.0895 + 1.5*6.0401 = 16.14965。
我知道负的低晶须是没有意义的,所以这里用最小值代替。但是为什么上面的胡须没有显示出来呢?我很困惑,在网上找不到解决这个问题的例子。一些我误解了的格格图设置?我真的很感谢你的帮助和建议!
发布于 2019-12-27 07:09:34
引用的部分如下:
上晶须从铰链延伸到最大值不超过1.5 * IQR从铰链(其中IQR是四分位数的范围,或第一和第三个四分位数之间的距离)。
所谓“值”,是指原始数据点中的值。如果绘制数据,在7.09和16.15 (+1.5*IQR)的顶部铰链之间没有值。如果这些四分位数是从数据中产生的,其中一个值位于这个范围内,那么上面的晶须就会出现在那里。
ggplot(data, aes(y = value)) +
geom_jitter(aes(x = 0.5), width = 0.05) +
stat_boxplot(geom = "errorbar", width = 0.3,
color = "red", size = 1.5) +
geom_boxplot(width = 0.5, alpha = 0.5) +
geom_hline(yintercept = c(7.09, 16.15), lty = "dashed")
https://stackoverflow.com/questions/59497080
复制相似问题