kmeans с расстоянием L1 в python

Question

kmeans с расстоянием L1 в python

Заданы векторы признаков NxM в виде матрицы numpy. Есть ли процедура, которая может кластерного алгоритма Kmeans с использованием дистанционных Л1 (расстояние Манхэттена)?

758 3

python numpy k-means

3 ответов:

Comments

Ничего не найдено.

JoshAdel · Accepted Answer · 2011-06-06 17:48:50

Я не думаю, что это предлагается явно в scipy, но вы должны взглянуть на следующее:

Http://projects.scipy.org/scipy/ticket/612

junwangbuaa · Accepted Answer · 2012-04-19 05:17:05

Вот один алгоритм Kmeans, использующий расстояние L1 (Manhattan distance). Для обобщения вектор признаков представлен в виде списка,который легко преобразовать в матрицу numpy.

    import random
    #Manhattan Distance
    def L1(v1,v2):
      if(len(v1)!=len(v2):
        print “error”
        return -1
      return sum([abs(v1[i]-v2[i]) for i in range(len(v1))])

    # kmeans with L1 distance. 
    # rows refers to the NxM feature vectors
    def kcluster(rows,distance=L1,k=4):# Cited from Programming Collective Intelligence 
        # Determine the minimum and maximum values for each point
        ranges=[(min([row[i] for row in rows]),max([row[i] for row in rows])) for i in range(len(rows[0]))]

        # Create k randomly placed centroids
        clusters=[[random.random( )*(ranges[i][1]-ranges[i][0])+ranges[i][0] for i in range(len(rows[0]))] for j in range(k)]

        lastmatches=None
        for t in range(100):
            print 'Iteration %d' % t
            bestmatches=[[] for i in range(k)]
            # Find which centroid is the closest for each row
            for j in range(len(rows)):
                row=rows[j]
                bestmatch=0
                for i in range(k):
                    d=distance(clusters[i],row)
                    if d<distance(clusters[bestmatch],row): 
                        bestmatch=i
                bestmatches[bestmatch].append(j)
            ## If the results are the same as last time, this is complete
            if bestmatches==lastmatches:
                break
            lastmatches=bestmatches

            # Move the centroids to the average of their members
            for i in range(k):
                avgs=[0.0]*len(rows[0])
                if len(bestmatches[i])>0:
                    for rowid in bestmatches[i]:
                        for m in range(len(rows[rowid])):
                            avgs[m]+=rows[rowid][m]
                    for j in range(len(avgs)):
                        avgs[j]/=len(bestmatches[i])
                    clusters[i]=avgs
        return bestmatches

denis · Accepted Answer · 2017-05-23 15:09:49

Есть код под -Это-можно-укажите свое расстояние-функция-через-scikits-учиться-методом k-средних, который использует любой из 20 с лишним показателей в составляющей.пространственный.расстояние. Смотреть также L1-или-L. 5-метрики для кластеризации ; не могли бы вы прокомментировать ваши результаты с L1 против L2 ?