我需要在scipy中使用normaltest来测试数据集是否为正态分布。但是我似乎找不到如何使用的好例子scipy.stats.normaltest。
scipy.stats.normaltest
我的数据集有100多个值。
In [12]: import scipy.stats as stats In [13]: x = stats.norm.rvs(size = 100) In [14]: stats.normaltest(x) Out[14]: (1.627533590094232, 0.44318552909231262)
normaltest返回卡方统计量的2元组以及关联的p值。给定x来自正态分布的零假设,p值表示看到较大(或更大)的卡方统计量的概率。
normaltest
x
如果p值很小,则意味着数据不太可能来自正态分布。例如:
In [15]: y = stats.uniform.rvs(size = 100) In [16]: stats.normaltest(y) Out[16]: (31.487039026711866, 1.4543748291516241e-07)