为什么 Jarque-Bera Test 正确值只返回 0？

Question

我测试 jarque_bera 的代码如下：

for x in data.columns.values:
    print(x, "-", stats.jarque_bera(data[x]))

结果是

这是否意味着我的数据帧的所有列都没有标准化？

Answer 1

jarque_bera测试正在测试零假设，即您的值以与正态分布相同的偏度和峰度分布。

发生的情况是，对于每个测试，您的测试统计量足够大，以至于 p 值非常小，因此 python 仅显示 0。Scipy 使用cdf计算 p 值，该值仅下降到e-16 。 使用sf可以给你小数，直到它下溢。

from scipy import stats

stats.chi2.cdf(1154.1959249402914, 2)  
#1.0

# scipy reported p-value
1 - stats.chi2.cdf(1154.1959249402914, 2)
#0.0

# This function gives us decimals
stats.chi2.sf(1154.1959249402914, 2)
#2.3417438136962163e-251