Question

我想通过我的数组（3d）在0轴上进行计算机伽玛拟合。我尝试过下面的代码，但有些不对劲。

import scipy.stats as st
import numpy as np
myarray = np.random.rand(20, 3, 3)
a, myloc, myscale = st.gamma.fit(myarray, axis=0)

Answer 1

fit方法不接受axis参数。

实际上，在scipy 0.15.1（可能是旧版本）中，它 “接受”它，但它忽略了它。这看起来像一个bug。当你给它任何实际上没有使用的论据时，它应该引发一个TypeError，但这不是它目前的工作方式：

In [20]: from scipy.stats import gamma

In [21]: gamma.fit([1,2,3,3,3,4,5,7])
Out[21]: (3.6899039741659925, 0.049374155400321737, 0.93515411814698823)

In [22]: gamma.fit([1,2,3,3,3,4,5,7], axis=123)
Out[22]: (3.6899039741659925, 0.049374155400321737, 0.93515411814698823)

In [23]: gamma.fit([1,2,3,3,3,4,5,7], flibbity=123)
Out[23]: (3.6899039741659925, 0.049374155400321737, 0.93515411814698823)

错误报告：https://github.com/scipy/scipy/issues/4932

Answer 2

应该是scipy.stats而不是scipy.stat。

Gamma适合轴3d阵列

2 个答案: