python - 如何使用 scipy(python) 计算伽马分布的熵？

Question

我使用 scipy 尝试了Saullo Castro 的示例代码来获取 gamma 分布的参数。

import scipy
import scipy.stats

size = 30000
x = scipy.arange(size)
y = scipy.int_(scipy.round_(scipy.stats.vonmises.rvs(5,size=size)*47))
dist = getattr(scipy.stats, 'gamma')
param = dist.fit(y)
dist.entropy(param)

给我熵结果，这是没有意义的。我在这里错在哪里？

array([ 5.54183717,         nan, -0.56062118])

在 scipy 文档中，它指出

entropy(a, loc=0, scale=1)  (Differential) entropy of the RV.

在哪里

a is the array of parameters of the gamma distribution