python - scipy.stats.spearmanr : Warning: divide by zero encountered in divide

Question

我使用 scipy.stats.spearmanr(a,b) 我得到:

Warning: divide by zero encountered in divide

操作正确结束，但显示警告。 a 和 b 都是“正常”数据(没有“仅零”向量等)。知道原因是什么，或者如何暂停警告？

编辑:

这是 spearmanr 中有问题的行:

/usr/lib/python2.7/dist-packages/scipy/stats/stats.pyc in spearmanr(a, b, axis)
   2226     rs = np.corrcoef(ar,br,rowvar=axisout)
   2227 
-> 2228     t = rs * np.sqrt((n-2) / ((rs+1.0)*(1.0-rs)))
   2229     prob = distributions.t.sf(np.abs(t),n-2)*2
   2230 

Answer 1

除以零是设计使然，对角线上的 rs=1。对于任何值都会发生这种情况。

但是，在 scipy 0.9 中，此错误已在 spearmanr 函数中本地消除。相应的sourceline包含

np.seterr(divide='ignore') # rs 可以有等于 1 的元素