julia-lang 中的直方图计算

Question

引用 julia-lang 文档:

hist(v[, n]) → e, counts

计算 v 的直方图，可选择使用大约 n 个 bin。返回值是一个范围 e，它对应于 bin 的边缘，以及包含每个 bin 中 v 元素数量的计数。注意:Julia 在计算中不会忽略 NaN 值。

我选择一个样本数据范围

testdata=0:1:10;

然后使用 hist 函数计算 1 到 5 个 bin 的直方图

hist(testdata,1) # => (-10.0:10.0:10.0,[1,10])
hist(testdata,2) # => (-5.0:5.0:10.0,[1,5,5])
hist(testdata,3) # => (-5.0:5.0:10.0,[1,5,5])
hist(testdata,4) # => (-5.0:5.0:10.0,[1,5,5])
hist(testdata,5) # => (-2.0:2.0:10.0,[1,2,2,2,2,2])

如您所见，当我想要 1 个 bin 时，它会计算 2 个 bin，而当我想要 2 个 bin 时，它会计算 3 个。

为什么会这样？

Answer 1

作为编写底层函数的人:目标是根据以 10 为基数的计数系统(即 10k、2×10k、5×10k)获得“不错”的 bin 宽度。如果您想要更多控制，您还可以指定确切的 bin 边缘。