matlab - Juliaの尖度関数

Question

Julia をいじってみたところ、確率分布の尖度を計算する関数が Julia と MATLAB で異なる方法で実装されていることがわかりました。

ジュリアで、次のようにします。

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

ここで何が起こっているのですか？2 つの言語で尖度が異なるのはなぜですか?

Answer 1

ここで説明されているように: http://www.itl.nist.gov/div898/handbook/eda/section3/eda35b.htm

正規分布の（超過）尖度がゼロになるように、過剰尖度（Kurtosis - 3）をよく使用します。Juliaで使用されているdistributions.jl docsに示されているように。kurtosis(x)

Matlab は過剰測定を使用しません (ドキュメントには、この潜在的な問題について言及しているメモさえあり​​ます)。