在Python中进行文本聚类

Question

我需要把一些文本文件分成不同的组，所以我在研究各种选择。看起来LingPipe可以直接对普通文本进行分组，而不需要先转换成其他格式（比如向量空间等），但这是我见过的唯一一个明确说可以处理字符串的工具。

有没有什么Python工具可以直接对文本进行分组？如果没有，那处理这个问题的最佳方法是什么呢？

Answer 1

有一个叫做NLTK的Python库，它可以帮助我们进行语言分析，包括把文本进行分类和聚类。

Answer 2

看起来可以通过一些简单的UNIX命令行工具，把这些文档里的文字内容提取到文本文件中，然后再用纯Python的方式来进行实际的聚类处理。

我找到了一段关于数据聚类的代码示例：

http://www.daniweb.com/code/snippet216641.html

这里有一个Python的库可以用来做这个：

http://python-cluster.sourceforge.net/

还有另一个Python库（主要用于生物信息学）：

http://bonsai.ims.u-tokyo.ac.jp/~mdehoon/software/cluster/software.htm#pycluster

Answer 3

文本聚类的质量主要取决于两个因素：

1. 你想要聚类的文档之间的相似性。例如，通过tfidf余弦距离，我们可以很容易地区分关于体育和政治的新闻文章。但如果要根据这个标准将产品评价分成“好”或“坏”，就难多了。

2. 聚类的方法本身。如果你知道会有多少个聚类，那就可以用k均值聚类。如果你不太在意准确性，只想展示一个漂亮的树状结构来导航搜索结果，那就可以用某种层次聚类的方法。

没有一种文本聚类的解决方案能在任何情况下都表现良好。因此，仅仅拿出一些聚类软件，把数据扔进去是不够的。

说到这里，这里有一些我之前用来玩文本聚类的实验代码。文档被表示为标准化的tfidf向量，相似性通过余弦距离来衡量。聚类的方法是majorclust。

import sys
from math import log, sqrt
from itertools import combinations

def cosine_distance(a, b):
    cos = 0.0
    a_tfidf = a["tfidf"]
    for token, tfidf in b["tfidf"].iteritems():
        if token in a_tfidf:
            cos += tfidf * a_tfidf[token]
    return cos

def normalize(features):
    norm = 1.0 / sqrt(sum(i**2 for i in features.itervalues()))
    for k, v in features.iteritems():
        features[k] = v * norm
    return features

def add_tfidf_to(documents):
    tokens = {}
    for id, doc in enumerate(documents):
        tf = {}
        doc["tfidf"] = {}
        doc_tokens = doc.get("tokens", [])
        for token in doc_tokens:
            tf[token] = tf.get(token, 0) + 1
        num_tokens = len(doc_tokens)
        if num_tokens > 0:
            for token, freq in tf.iteritems():
                tokens.setdefault(token, []).append((id, float(freq) / num_tokens))

    doc_count = float(len(documents))
    for token, docs in tokens.iteritems():
        idf = log(doc_count / len(docs))
        for id, tf in docs:
            tfidf = tf * idf
            if tfidf > 0:
                documents[id]["tfidf"][token] = tfidf

    for doc in documents:
        doc["tfidf"] = normalize(doc["tfidf"])

def choose_cluster(node, cluster_lookup, edges):
    new = cluster_lookup[node]
    if node in edges:
        seen, num_seen = {}, {}
        for target, weight in edges.get(node, []):
            seen[cluster_lookup[target]] = seen.get(
                cluster_lookup[target], 0.0) + weight
        for k, v in seen.iteritems():
            num_seen.setdefault(v, []).append(k)
        new = num_seen[max(num_seen)][0]
    return new

def majorclust(graph):
    cluster_lookup = dict((node, i) for i, node in enumerate(graph.nodes))

    count = 0
    movements = set()
    finished = False
    while not finished:
        finished = True
        for node in graph.nodes:
            new = choose_cluster(node, cluster_lookup, graph.edges)
            move = (node, cluster_lookup[node], new)
            if new != cluster_lookup[node] and move not in movements:
                movements.add(move)
                cluster_lookup[node] = new
                finished = False

    clusters = {}
    for k, v in cluster_lookup.iteritems():
        clusters.setdefault(v, []).append(k)

    return clusters.values()

def get_distance_graph(documents):
    class Graph(object):
        def __init__(self):
            self.edges = {}

        def add_edge(self, n1, n2, w):
            self.edges.setdefault(n1, []).append((n2, w))
            self.edges.setdefault(n2, []).append((n1, w))

    graph = Graph()
    doc_ids = range(len(documents))
    graph.nodes = set(doc_ids)
    for a, b in combinations(doc_ids, 2):
        graph.add_edge(a, b, cosine_distance(documents[a], documents[b]))
    return graph

def get_documents():
    texts = [
        "foo blub baz",
        "foo bar baz",
        "asdf bsdf csdf",
        "foo bab blub",
        "csdf hddf kjtz",
        "123 456 890",
        "321 890 456 foo",
        "123 890 uiop",
    ]
    return [{"text": text, "tokens": text.split()}
             for i, text in enumerate(texts)]

def main(args):
    documents = get_documents()
    add_tfidf_to(documents)
    dist_graph = get_distance_graph(documents)

    for cluster in majorclust(dist_graph):
        print "========="
        for doc_id in cluster:
            print documents[doc_id]["text"]

if __name__ == '__main__':
    main(sys.argv)

在实际应用中，你应该使用一个合适的分词器，使用整数而不是字符串来表示词汇，并且不要计算O(n^2)的距离矩阵……