我无法理解为什么scipy.stats.hypergeom.pmf()
根据它own documentation nan
引发import scipy.stats as stats
(M, n, N) = (1, 0, 1)
pmf = stats.hypergeom.pmf(np.array([0, 1]), M=M, n=n, N=N)
print pmf
时似乎不应该这样做。
考虑这个代码块:
[ nan nan]
输出为pmf(k, M, n, N) = choose(n, k) * choose(M - n, N - k) / choose(M, N)
。
然而,根据文件:
概率质量函数定义为:
scipy
当我使用nan
运行该代码时,我得到了正确的非from scipy.misc import comb
(M, n, N) = (1, 0, 1)
k=np.array([0, 1])
pmfGood = comb(n, k) * comb(M - n, N - k) / comb(M, N)
print pmfGood
[ 1. 0.]
答案:
scipy.stats.hypergeom.pmf()
以下是我使用的软件包的一些版本信息:
我使用Scipy版本0.15.1; Numpy版本1.8.2;和Python 2.7.9 [Anaconda 2.1.0(x86_64)]
这是一个错误吗?我不理解{{1}} ??