python – Scipy稀疏 – 距离矩阵(Scikit或Scipy)

前端之家收集整理的这篇文章主要介绍了python – Scipy稀疏 – 距离矩阵(Scikit或Scipy)前端之家小编觉得挺不错的,现在分享给大家,也给大家做个参考。

我试图在scikit-learn的DictVectorizer返回的Scipy稀疏矩阵上计算最近邻居聚类.但是,当我尝试使用scikit-learn计算距离矩阵时,我通过pairwise.euclidean_distances和pairwise.pairwise_distances使用’euclidean’距离得到错误消息.我的印象是scikit-learn可以计算这些距离矩阵.

我的矩阵非常稀疏,形状为:< 364402x223209稀疏矩阵类型< class'numpy.float64'>
    使用压缩稀疏行格式的728804存储元素>.

我也在Scipy中尝试了诸如pdist和kdtree之类的方法,但是还收到了其他无法处理结果的错误.

任何人都可以请我指出一个有效地允许我计算距离矩阵和/或最近邻结果的解决方案吗?

一些示例代码

  1. import numpy as np
  2. from sklearn.feature_extraction import DictVectorizer
  3. from sklearn.neighbors import NearestNeighbors
  4. from sklearn.metrics import pairwise
  5. import scipy.spatial
  6. file = 'FileLocation'
  7. data = []
  8. FILE = open(file,'r')
  9. for line in FILE:
  10. templine = line.strip().split(',')
  11. data.append({'user':str(int(templine[0])),str(int(templine[1])):int(templine[2])})
  12. FILE.close()
  13. vec = DictVectorizer()
  14. X = vec.fit_transform(data)
  15. result = scipy.spatial.KDTree(X)

错误

  1. Traceback (most recent call last):
  2. File "

同样,如果我跑:

  1. scipy.spatial.distance.pdist(X,'euclidean')

我得到以下内容

  1. Traceback (most recent call last):
  2. File "

最后,在scikit-learn中运行NearestNeighbor会导致内存错误,使用:

  1. nbrs = NearestNeighbors(n_neighbors=10,algorithm='brute')
最佳答案
首先,你不能使用稀疏矩阵的KDTree和pdist,你必须将它转换为密集(你的选择是否是你的选择):

  1. >>> X
  2. <2x3 sparse matrix of type '

第二,从the docs

Efficient brute-force neighbors searches can be very competitive for small data samples. However,as the number of samples N grows,the brute-force approach quickly becomes infeasible.

您可能想尝试’ball_tree’算法并查看它是否可以处理您的数据.

猜你在找的Python相关文章