稀疏输入的非负矩阵分解

时间:2013-09-02 14:46:15

标签: matrix factorization

任何人都可以推荐一套工具来对稀疏输入数据[大小为50kx50k的矩阵]执行标准NMF应用程序,谢谢!

1 个答案:

答案 0 :(得分:5)

scikit-learn有NMF for sparse matrices的实现。但是,您将需要GitHub的最新版本,因为所有已发布的版本(最多包括0.14)都存在可伸缩性问题。接下来是一个演示。

加载一些数据:二十个新闻组文本语料库。

>>> from sklearn.datasets import fetch_20newsgroups
>>> from sklearn.decomposition import NMF
>>> from sklearn.feature_extraction.text import CountVectorizer
>>> from sklearn.preprocessing import normalize
>>> data = fetch_20newsgroups().data
>>> X = CountVectorizer(dtype=float).fit_transform(data)
>>> X = normalize(X)
>>> X
<11314x130107 sparse matrix of type '<type 'numpy.float64'>'
    with 1787565 stored elements in Compressed Sparse Column format>

现在适用于包含10个组件的NMF模型。

>>> nmf = NMF(n_components=10, tol=.01)
>>> Xnmf = nmf.fit_transform(X)

我调整了公差选项,以便在几秒钟内完成收敛。使用默认容差,需要更长的时间。此问题的内存使用量约为360MB。

免责声明:我是此实施的贡献者,因此这不是公正的建议。