运行主函数时出现断言错误

Question

我用python实现kmeans算法，代码如下。我用一些简单的数据测试代码。如下所示，哪个存储在一个名为数据.txt
25
37
-1-2
-3-3
54
4-4
3-7
3.5-9节

我的问题是，在迭代过程中，有些簇看起来是空的，即（簇数）<；k，经过我的分析，似乎会发生这种情况，但是在搜索web之后，我发现kmeans算法中没有人处理这个问题。在

所以我不知道错在哪里？是因为我的测试数据很简单吗

import sys
import numpy as np
from math import sqrt

"""
useage: python mykmeans.py mydata.txt k

"""

GAP = 2
MIN_VAL = 1000000

def get_distance(point1, point2):
    dis = sqrt(pow(point1[0] - point2[0], 2) + pow(point1[1] - point2[1], 2))

    return dis


def cluster_dis(centroid, cluster):
    dis = 0.0
    for point in cluster:
        dis += get_distance(centroid, point)

    return dis

def update_centroids(centroids, cluster_id, cluster):
    x, y = 0.0, 0.0
    length = len(cluster)
    if length == 0:  # TODO： this is my question? do we need to examine this?
        return

    for item in cluster:
        x += item[0]
        y += item[1]
    centroids[cluster_id] = (x / length, y / length)


def kmeans(data, k):
    assert k <= len(data)

    seed_ids = np.random.randint(0, len(data), k)
    centroids = [data[idx] for idx in seed_ids]
    clusters = [[] for _ in xrange(k)]
    cluster_idx = [-1] * len(data)

    pre_dis = 0
    while True:
        for point_id, point in enumerate(data):
            min_distance, tmp_id = MIN_VAL, -1
            for seed_id, seed in enumerate(centroids):
                distance = get_distance(seed, point)
                if distance < min_distance:
                    min_distance = distance
                    tmp_id = seed_id
            if cluster_idx[point_id] != -1:
                dex = clusters[cluster_idx[point_id]].index(point)
                del clusters[cluster_idx[point_id]][dex]
            clusters[tmp_id].append(point)
            cluster_idx[point_id] = tmp_id

        now_dis = 0.0
        for cluster_id, cluster in enumerate(clusters):
            now_dis += cluster_dis(centroids[cluster_id], cluster)
            update_centroids(centroids, cluster_id, cluster)

        delta_dis = now_dis - pre_dis
        pre_dis = now_dis

        if delta_dis < GAP:
            break

    print(centroids)
    print(clusters)

    return centroids, clusters

def get_data(file_name):
    try:
        fr = open(file_name)
        lines = fr.read().splitlines()
    except IOError, e:
        pass
    finally:
        fr.close()

    data = []
    for line in lines:
        tmp = line.split()
        x, y = float(tmp[0]), float(tmp[1])
        data.append([x, y])

    return data

def main():
    args = sys.argv[1:]
    assert len(args) > 1
    file_name, k = args[0], int(args[1])

    data = get_data(file_name)
    kmeans(data, k)


if __name__ == '__main__':
    main()

Answer 1

k-means有可能诱导一个空簇。这是one example如图所示。我还复制了下面的数字，以防链接某天到期。在

下面的第一张图显示了7个点的分布情况。最初选择3、5和6作为群集中心。在

下面的“+”表示第一次迭代后簇中心的变化，同一颜色表示对应点在同一簇中。在

从下图中可以看到，经过2次迭代后，蓝色的集群变为空，实际上有2个集群，而不是初始化值3。在

因此，空簇可能是由于初始化和“不正确”的簇号造成的。您可以在代码中尝试不同的k，并多次运行该程序以观察聚类结果，从而使其更加健壮。在