r - R中箱线图上的平均菱形图

Question

我有一个带有代表平均值的点的箱线图：

这没关系，但我想显示平均值周围的置信区间，例如以下列方式：

如何在 R 中叠加这两个图？

Answer 1

像下面这样的函数很容易编写。输入菱形的平均值、上限、下限和宽度的 x 和 y。您也可以将其他参数（例如用于线型的 lty 和 col 颜色）传递给线段。

diamondCI <- function(x, y, ul, ll, w = ul-ll, ...){
    hw <- w/2
    segments(x-hw, y, x+hw, y, ...)  # horizontal bar
    segments(x-hw, y, x, ul, ...) # left upper diag
    segments(x, ul, x+hw, y, ...) # right upper diag
    segments(x-hw, y, x, ll, ...) # left lower diag
    segments(x, ll, x+hw, y, ...) # right lwoer diag
    }

尝试一下...

plot(1,1)
diamondCI(1, 1, 1.2, 0.7, col = 'red', lwd = 3)