matlab - Julia中的峰度函数

Question

所以我一直在玩 Julia，我发现计算概率分布峰度的函数在 Julia 和 MATLAB 之间的实现方式不同。

在 Julia 中，执行以下操作：

using Distributions
dist = Beta(3, 5)
x = rand(dist, 10000)
kurtosis(x) #gives a value approximately around -0.42

在 MATLAB 中执行：

x = betarnd(3, 5, [1, 10000]);
kurtosis(x) %gives something approximately around 2.60

这里发生了什么事？为什么两种语言之间的峰度不同？

Answer 1

如此处所述：http ://www.itl.nist.gov/div898/handbook/eda/section3/eda35b.htm

我们经常使用过度峰度 (Kurtosis - 3) 以使正态分布的 (Excess) Kurtosis 变为零。如Julia中使用的distributions.jl 文档所示。kurtosis(x)

Matlab 不使用过度测量（文档中甚至有一条注释提到了这个潜在问题）。