1

R で箱ひげ図を作成したいデータ セットがあります。R コンソールで次のコマンドを使用して、簡単に行うことができました (ここではサンプル データのみを使用しました)。

x <- rnorm(100);

y <- rnorm(100);

boxplot(x, y, col=rainbow(2))

ここに画像の説明を入力

今、私が本当に望んでいるのは、上記のボックスプロットの上にデータセットの (平均 +/- (2 または 3)*シグマ) マークを表示することです。私は標準偏差と平均が堅牢な統計ではないことを知っていますが、私のデータセットではその情報が役立つ場合があります。つまり、意図した図は次のようになります。

ここに画像の説明を入力 (異なる色の) 点線は、対応するプロットの 95% または 99.7% 信頼区間、または (平均 +/- (2 または 3)*シグマ) を表します。

注意: 残念ながら、私は十分な評判を持っていないため、ここに数値を追加することはできません。誰かが上記のリンクからここにアップロードしていただければ幸いです.

4

1 に答える 1

2

非常にシンプルで簡単にできることは、関数でnotch=TRUE引数を使用することです ( ?boxplotを参照)。これは、平均の 95% 信頼区間とまったく同じではありませんが、代わりに中央値の 95% CI に類似しています (これは、堅牢性に関する懸念に対処します)。ここで説明します： boxplot()

ノッチ (要求された場合) は、+/-1.58 IQR/sqrt(n) まで拡張されます。これは、Chambers らの 1.57 の式と同じ計算に基づいているようです。(1983, p. 62)、McGill et al. (1978 年、16 ページ)。それらは、中央値の漸近的正規性と、比較される 2 つの中央値のサンプル サイズがほぼ等しいことに基づいており、サンプルの基になる分布にかなり鈍感であると言われています。アイデアは、2 つの中央値の差におよそ 95% の信頼区間を与えることであると思われます。

于 2013-06-17T19:23:44.723 に答える