2012-03-14 20 views
3

で箱ひげ図上にダイヤモンドプロット平均:私は平均を表すドットと箱ひげ図を持っているR

Boxplot with mean

これは大丈夫ですが、私は、例えば、平均値の周りの信頼区間を表示したいです次のように:

Diamond plot

どのように私はRでこれら二つのプロットを重ねることができますか?

+0

...でそれを試してみてくださいx軸上で? – John

+0

私は質問を更新しました。私は2つのプロットをオーバーレイしたいと思います。 – Mulone

答えて

3

以下のような機能は、書きやすいです。 ダイヤモンドの平均値、上限値、下限値、および幅の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 
    } 

あなたはCIの限界は、上記の箱ひげ図のように縦方向に位置することを意味するのですか、あなたは推定している別の変数を持っています下の図では

plot(1,1) 
diamondCI(1, 1, 1.2, 0.7, col = 'red', lwd = 3) 
関連する問題