python - 如何执行 kolmogorov-Smirnov 测试以针对 Gamma 函数进行测试？

Question

我有一个点数组，我想检查它们是否服从 a=15.5 和 b=7 的 Gamma 分布。我在努力

scipy.stats.kstest(myarray, gamma.pdf(15.5, 0.143))

kstest 命令的语法是什么？该文档仅显示了正态分布的基本示例。

Answer 1

要使用 Gamma 分布，将字符串 'gamma' 作为第二个参数传递，并使用 args 参数设置形状参数、位置和规模。

这是一个例子。首先，我需要一些数据来处理:

In [46]: from scipy.stats import gamma, kstest

In [47]: np.random.seed(123)

In [48]: sample = gamma.rvs(a=15.5, loc=0, scale=1./7, size=100)

应用具有形状 15.5、位置 0 和比例 1/7 的 Gamma 分布的检验。因为 args 给出了位置参数，所以我们必须包括位置:

In [49]: kstest(sample, 'gamma', args=(15.5, 0, 1./7))
Out[49]: KstestResult(statistic=0.058075855762720785, pvalue=0.88868922608343937)

如果我们尝试不同的尺度，我们会看到 p 值下降:

In [50]: kstest(sample, 'gamma', args=(15.5, 0, 0.16))
Out[50]: KstestResult(statistic=0.18210965421938488, pvalue=0.002249786773449225)

In [51]: kstest(sample, 'gamma', args=(15.5, 0, 0.175))
Out[51]: KstestResult(statistic=0.29663539668929456, pvalue=2.6970075950316641e-08)