jarque bera測試結果

Question

我不了解Jarque Bera測試的結果。

from statsmodels.stats.stattools import jarque_bera    
np.random.seed(123)
jarque_bera(np.random.normal(-5, 1, 1000))

結果：

(0.1675179797931011,
 0.9196528750223983,
 -0.029040113501245704,
 2.9745614712223074)

第三值看起來像P值。 我認為其他的是峰度和偏斜，而我不知道的是第四。

所以我測試了我的理論，但是按照下面的代碼是錯誤的：

import scipy.stats as stats

print(stats.skew(np.random.normal(-5, 1, 1000)))
print(stats.kurtosis(np.random.normal(-5, 1, 1000)))

-0.19743173433793879
-0.11038007419823126

Answer 1

您需要n> 2000才能使Jarque Bera測試有效

輸出給您； 測試stat，p值，偏斜，峰度按該順序排列。 不知道為什么這不在文檔中嗎？

同樣，已實施的Jarque Bera測試使用的是Pearson對峰度的定義，而不是Fisher的定義 ，因此...

from statsmodels.stats.stattools import jarque_bera  
import scipy.stats as stats
import numpy as np

np.random.seed(123)
samples = np.random.normal(-5, 1, 3000)

print(jarque_bera(samples))
print(stats.skew(samples))
print(stats.kurtosis(samples, fisher=False))

輸出..

(3.9600892567754835, 0.13806307564092868, -0.08899286958111645, 3.0013381737844793)
-0.08899286958111645
3.0013381737844793