为什么JarqueBera测试正确值只返回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值，该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