ValueError: qk and pk must have same shape - scipy.spatial.distance.jensenshannon

Question

我正在调用下面的 jensen_shannon(query,matrix) 函数来查找文档矩阵中文档查询的最相似文档

def jensen_shannon(query, matrix):
"""
This function implements a Jensen-Shannon similarity
between the input query (an LDA topic distribution for a document)
and the entire corpus of topic distributions.
It returns an array of length M where M is the number of documents in the corpus
"""
# lets keep with the p,q notation above
p = query[None,:].T # take transpose
q = matrix.T # transpose matrix
m = 0.5*(p + q)
return np.sqrt(0.5*(entropy(p,m) + entropy(q,m)))

查询形状：(100,)

矩阵的形状：(10804, 100)

错误回溯：

ValueError                                Traceback (most recent call last)
<ipython-input-103-86cb68dd862d> in <module>
      1 # this is surprisingly fast
----> 2 most_sim_ids = get_most_similar_documents(new_doc_distribution,doc_topic_dist)

<ipython-input-102-c0fb95224e87> in get_most_similar_documents(query, matrix, k)
      6     print(query.shape)
      7     print(matrix.shape)
----> 8     sims = jensen_shannon(query,matrix) # list of jensen shannon distances
      9     return sims.argsort()[:k] # the top k positional index of the smallest Jensen Shannon distances

<ipython-input-74-6ffb0ec54e9a> in jensen_shannon(query, matrix)
     10     q = matrix.T # transpose matrix
     11     m = 0.5*(p + q)
---> 12     return np.sqrt(0.5*(entropy(p,m) + entropy(q,m)))

~/venv/lib/python3.6/site-packages/scipy/stats/_distn_infrastructure.py in entropy(pk, qk, base, axis)
   2668         qk = asarray(qk)
   2669         if qk.shape != pk.shape:
-> 2670             raise ValueError("qk and pk must have same shape.")
   2671         qk = 1.0*qk / np.sum(qk, axis=axis, keepdims=True)
   2672         vec = rel_entr(pk, qk)

ValueError: qk and pk must have same shape.

Add axis parameter for scipy.spatial.distance.jensenshannon 但它不接受函数中的轴参数。

有人知道我错过了什么吗？任何领导高度赞赏。谢谢

仅供参考：我正在尝试这个 kaggle 代码 https://www.kaggle.com/ktattan/lda-and-document-similarity/data

Answer 1

运行我自己进去了。 scipy 版本 1.3.0 的行为仍然符合您使用的 Jensen-Shannon 公式的预期。

Answer 2

试试这个：

p = query[None,:].T + np.zeros([100, 10804])

100 = 主题数 10804 = 文件数量

ValueError: qk and pk must have same shape - scipy.spatial.distance.jensenshannon

ValueError: qk and pk must have same shape - scipy.spatial.distance.jensenshannon

python

scipy

topic-modeling