R中箱形图的平均钻石图

Question

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

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

如何在R？

中叠加这两个图

Answer 1

以下功能易于编写。 输入x和y作为钻石的平均值，上限，下限和宽度。您也可以将其他参数（如行类型和col颜色等lty）传递给段。

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)