如何執行 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)

Answer 2

from scipy import stats
stats.kstest(myarray, 'gamma', (15.5, 7))