r - 在 R 中绘制柯西分布的直方图

Question

我正在尝试使用以下代码在 R 中绘制柯西分布的直方图：

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

而且无论我在hist()函数中尝试什么选项，我在直方图上都看不到两个以上的条（基本上一个用于负值，一个用于正值）。

但是，当我使用正态分布（或其他）时，它工作正常。

Answer 1

这是由于分布的特性造成的。

大多数值都相对接近于零，但非常大的绝对值比正态分布更有可能。大约有 1% 的值的绝对值大于 50，而 0.1% 的值大于 500。

尝试仅绘制部分值：

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

您还可以查看累积分布函数：

plot(ecdf(X))

并检查箱线图：

boxplot(X)