Scipy Normaltest如何使用？

Question

我需要在scipy中使用normaltest来测试数据集是否为正态分布。但是我似乎找不到如何使用的好例子scipy.stats.normaltest。

我的数据集有100多个值。

Answer 1

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值表示看到较大（或更大）的卡方统计量的概率。

如果p值很小，则意味着数据不太可能来自正态分布。例如：

In [15]: y = stats.uniform.rvs(size = 100)

In [16]: stats.normaltest(y)
Out[16]: (31.487039026711866, 1.4543748291516241e-07)