3

次のコードを使用して、R でコーシー分布のヒストグラムをプロットしようとしています。

X = rcauchy(10^5)
hist(X)

関数でどのオプションを試してhist()も、ヒストグラムに 2 つ以上のバーが表示されることはありません (基本的に、負の値に対して 1 つ、正の値に対して 1 つ)。

ただし、正規分布（またはその他）を使用すると問題なく動作します。

4

1 に答える 1

5

これは分布の特性によるものです。

ほとんどの値は比較的ゼロに近くなりますが、絶対値が非常に大きい場合は、正規分布よりも確率が高くなります。絶対値が 50 より大きい約 1 % の値と、500 より大きい 0.1 % の値があります。

値の一部だけをプロットしてみてください:

hist(X[abs(X)<1])
hist(X[abs(X)<5])
hist(X[abs(X)<50])
hist(X)

累積分布関数も確認できます。

plot(ecdf(X))

ボックスプロットを確認します。

boxplot(X)
于 2013-01-26T10:47:36.817 に答える