如何使用scipy（python）计算伽玛分布的熵？

Question

我尝试了Saullo Castro的sample code来获取伽玛分布的参数，使用scipy。

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