r - Rの箱ひげ図の平均ダイヤモンドプロット

Question

平均を表すドットを含む箱ひげ図があります。

これで問題ありませんが、たとえば次のように、平均値の信頼区間を表示したいと思います。

Rでこれら2つのプロットを重ねるにはどうすればよいですか?

Answer 1

次のような関数は簡単に記述できます。ダイヤモンドの平均、上限、下限、および幅の x と y を入力します。線種の 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)