【发布时间】:2016-04-24 16:26:56
【问题描述】:
我编写了以下代码来演示平坦的统计分布:
m <- trunc(runif(100000,1,21))
table(m)
hist(m)
表 (m) 中的值与预期一致:
5010 5053 4951 5030 5034 5131 4895 5032 4974 5073 4947 5030 4939 5011 5069
4910 4992 4990 4946 4983
但是 hist (m) 在第一个 bin 中显示的计数是预期的两倍:
我尝试使用 20、21 或其他值设置中断:hist (m, breaks = 21),但问题仍然存在。
我认为这可能与“漂亮”运算符有关,但花了一个小时研究这个,我仍然看不出如何获得合理的直方图。
【问题讨论】:
-
有关如何计算休息时间的详细说明,您可以查看this post。中断计算过程从
ceiling(log2(length(x)) + 1)开始,随后传递给pretty,如链接文章中所述。