关于HIST函数的问题

clanyibin

我使用hist函数时，breaks=10时，最后作图的分组却是11组，这是为何？请大家指教！

jcwang

x=rnorm(1000)

#使用hist函数时，如果希望分成10组，则

hist(x,breaks=seq(-3,3,length.out=11))

#or

hist(x,breaks=seq(min(x),max(x),length.out=11))

clanyibin

回复第2楼的 jcwang：谢谢你，现在可以了，但是不知道这是啥原理啊，breaks不就是限定组数的吗？为啥直接输入就不行呢？

yufree

如果你看一下hist的原函数

getS3method('hist', 'default')

</p>

会发现breaks的参数是通过pretty来给出的

<br />
breaks <- pretty(range(x), n = breaks, min.n = 1)<br />
nB <- length(breaks)<br />

</p>

pretty函数如其名比较文艺原函数也很文艺我没看懂但这个函数决定了breaks的飘忽不定

jcwang

4楼言之有理，另外补充一点。

如果breaks=11，那么分组数counts就是10。

这就好比是在11根电线杆中共有10个间距。

你也可以作个验证：

x=rnorm(1000)

h=hist(x,breaks=seq(min(x),max(x),length.out=11))

length(h$breaks)-1==length(h$counts)

superdesolator

回复第5楼的 jcwang：这里你那个命令图作出来貌似不太好看：

> x=rnorm(100)

> hist(x,breaks=seq(min(x),max(x),length.out=11))

你试下，我也很奇怪为什么你前面要设定-3到3。。

可以range(x)一下再调下两端么。。。

jcwang

(1)

我写的是x=rnorm(1000)而不是x=rnorm(100)

(2)

不理解"貌似不太好看"的意思。

(3)

range(x)等于c(min(x),max(x))，因为

range(x)==c(min(x),max(x))

min(x)和max(x)可以分别写成range(x)[1]和range(x)[2]

所以下面两个写法画出的图完全一样

hist(x,breaks=seq(min(x),max(x),length.out=11))

hist(x,breaks=seq(range(x)[1],range(x)[2],length.out=11))。

(4)

关于-3,3的问题，先更正一下，应改为-4,4。

hist(x,breaks=seq(-4,4,length.out=11))

其中的-4,4是用来指定图中x轴的最小值为-4，最大值为4，

因为在<-4和>4这两个区间内几乎看不出图中直条的高度了，

所有干脆略去不画，反而绘图效果更好(不妨与下句比较一下)

hist(x,breaks=seq(-5,5,length.out=11))

这样写就是为了满足1楼所提出的要求（分成10组）。

如果不要求分多少组，写起来就省事多了，hist(x)即可。