我正在尝试使用以下代码在R中绘制Cauchy分布的直方图:
X = rcauchy(10^5)
hist(X)
无论我在hist()
函数中尝试了哪些选项,我的直方图上都看不到两个以上的条形图(基本上一个用于负值,一个用于正值)。
但是,当我使用正态分布(或其他)时,它工作正常。
答案 0 :(得分:5)
这是由分布的属性产生的。
大多数值相对接近于零,但非常大的绝对值比正态分布更可能。大约1%的值的绝对值大于50,大于500的0.1%。
尝试仅绘制部分值:
hist(X[abs(X)<1])
hist(X[abs(X)<5])
hist(X[abs(X)<50])
hist(X)
您还可以查看累积分布函数:
plot(ecdf(X))
并查看boxplot:
boxplot(X)