Як побудувати графік виведення даних кластеризації?


15

Я спробував згрупувати набір даних (набір знаків) і отримав 2 кластери. Я хотів би це графічно представити. Трохи розгублений у поданні, оскільки у мене немає координат (x, y).

Також шукає функцію MATLAB / Python для цього.

EDIT

Я думаю, що розміщення даних робить питання зрозумілішим. У мене є два кластери, які я створив, використовуючи кластеризацію kmeans в Python (не використовуючи scipy). Вони є

class 1: a=[3222403552.0, 3222493472.0, 3222491808.0, 3222489152.0, 3222413632.0, 
3222394528.0, 3222414976.0, 3222522768.0, 3222403552.0, 3222498896.0, 3222541408.0, 
3222403552.0, 3222402816.0, 3222588192.0, 3222403552.0, 3222410272.0, 3222394560.0, 
3222402704.0, 3222298192.0, 3222409264.0, 3222414688.0, 3222522512.0, 3222404096.0, 
3222486720.0, 3222403968.0, 3222486368.0, 3222376320.0, 3222522896.0, 3222403552.0, 
3222374480.0, 3222491648.0, 3222543024.0, 3222376848.0, 3222403552.0, 3222591616.0, 
3222376944.0, 3222325568.0, 3222488864.0, 3222548416.0, 3222424176.0, 3222415024.0, 
3222403552.0, 3222407504.0, 3222489584.0, 3222407872.0, 3222402736.0, 3222402032.0, 
3222410208.0, 3222414816.0, 3222523024.0, 3222552656.0, 3222487168.0, 3222403728.0, 
3222319440.0, 3222375840.0, 3222325136.0, 3222311568.0, 3222491984.0, 3222542032.0, 
3222539984.0, 3222522256.0, 3222588336.0, 3222316784.0, 3222488304.0, 3222351360.0, 
3222545536.0, 3222323728.0, 3222413824.0, 3222415120.0, 3222403552.0, 3222514624.0, 
3222408000.0, 3222413856.0, 3222408640.0, 3222377072.0, 3222324304.0, 3222524016.0, 
3222324000.0, 3222489808.0, 3222403552.0, 3223571920.0, 3222522384.0, 3222319712.0, 
3222374512.0, 3222375456.0, 3222489968.0, 3222492752.0, 3222413920.0, 3222394448.0, 
3222403552.0, 3222403552.0, 3222540576.0, 3222407408.0, 3222415072.0, 3222388272.0, 
3222549264.0, 3222325280.0, 3222548208.0, 3222298608.0, 3222413760.0, 3222409408.0, 
3222542528.0, 3222473296.0, 3222428384.0, 3222413696.0, 3222486224.0, 3222361280.0, 
3222522640.0, 3222492080.0, 3222472144.0, 3222376560.0, 3222378736.0, 3222364544.0, 
3222407776.0, 3222359872.0, 3222492928.0, 3222440496.0, 3222499408.0, 3222450272.0, 
3222351904.0, 3222352480.0, 3222413952.0, 3222556416.0, 3222410304.0, 3222399984.0, 
3222494736.0, 3222388288.0, 3222403552.0, 3222323824.0, 3222523616.0, 3222394656.0, 
3222404672.0, 3222405984.0, 3222490432.0, 3222407296.0, 3222394720.0, 3222596624.0, 
3222597520.0, 3222598048.0, 3222403552.0, 3222403552.0, 3222403552.0, 3222324448.0, 
3222408976.0, 3222448160.0, 3222366320.0, 3222489344.0, 3222403552.0, 3222494480.0, 
3222382032.0, 3222450432.0, 3222352000.0, 3222352528.0, 3222414032.0, 3222728448.0, 
3222299456.0, 3222400016.0, 3222495056.0, 3222388848.0, 3222403552.0, 3222487568.0, 
3222523744.0, 3222394624.0, 3222408112.0, 3222406496.0, 3222405616.0, 3222592160.0, 
3222549360.0, 3222438560.0, 3222597024.0, 3222597616.0, 3222598128.0, 3222403552.0, 
3222403552.0, 3222403552.0, 3222499056.0, 3222408512.0, 3222402064.0, 3222368992.0, 
3222511376.0, 3222414624.0, 3222554816.0, 3222494608.0, 3222449792.0, 3222351952.0, 
3222352272.0, 3222394736.0, 3222311856.0, 3222414288.0, 3222402448.0, 3222401056.0, 
3222413568.0, 3222298848.0, 3222297184.0, 3222488000.0, 3222490528.0, 3222394688.0, 
3222408224.0, 3222406672.0, 3222404896.0, 3222443120.0, 3222403552.0, 3222596400.0, 
3222597120.0, 3222597712.0, 3222400896.0, 3222403552.0, 3222403552.0, 3222403552.0, 
3222299200.0, 3222321296.0, 3222364176.0, 3222602208.0, 3222513040.0, 3222414656.0, 
3222564864.0, 3222407904.0, 3222449984.0, 3222352096.0, 3222352432.0, 3222452832.0, 
3222368560.0, 3222414368.0, 3222399376.0, 3222298352.0, 3222573152.0, 3222438080.0, 
3222409168.0, 3222523488.0, 3222394592.0, 3222405136.0, 3222490624.0, 3222406928.0, 
3222407104.0, 3222442464.0, 3222403552.0, 3222596512.0, 3222597216.0, 3222597968.0, 
3222438208.0, 3222403552.0, 3222403552.0, 3222403552.0]

class 2: b=[3498543128.0, 3498542920.0, 3498543252.0, 3498543752.0, 3498544872.0, 
3498544528.0, 3498543024.0, 3498542548.0, 3498542232.0]

Я хотів би побудувати це. Я спробував наступне і отримав наступний результат, коли я задумав aі b.

pylab.plot(a,'x')
pylab.plot(b,'o')
pylab.show()

введіть тут опис зображення

чи можу я отримати кращу візуалізацію кластеризації?


1
Це насправді залежить від того, що ви зробили кластеризацію :) Якщо ви покажете невеликий приклад даних, які у вас є, я впевнений, що ви отримаєте відповідь
david w

1
Використання різних кольорів та маркерів, як правило, є найпростішим / найпростішим для читання. Якщо у вас є два кластери, ви можете просто надрукувати 0/1 або O / X для різних значень.
Марцін

Скажіть, будь ласка, що ви маєте на увазі під "набором знаків". Скільки змінних у вас є, щоб охарактеризувати кластери? Також ви впевнені, що 2 - це найкраща кількість кластерів для використання? Багато разів доводиться ітеративно використовувати програми кластерного аналізу; спочатку можна отримати всього 2, але за допомогою деяких коригувань ви можете отримати більш цікаве та інформативне вище число.
rolando2

Я використовував kmeans, де я маю чітко
вказати

@ user2721, чи можете ви показати нам, як ви використовуєте kmeans?
Сигур

Відповіді:


30

Зазвичай ви плануєте початкові значення в розсипці (або матриці розсипань, якщо їх у вас багато), і використовуєте колір, щоб показати свої групи.

Ви попросили відповіді в python, і ви фактично все кластеризуєте та плануєте за допомогою scipy, numpy та matplotlib:

Почніть з внесення деяких даних

import numpy as np
from scipy import cluster
from matplotlib import pyplot

np.random.seed(123)
tests = np.reshape( np.random.uniform(0,100,60), (30,2) )
#tests[1:4]
#array([[ 22.68514536,  55.13147691],
#       [ 71.94689698,  42.31064601],
#       [ 98.07641984,  68.48297386]])

Скільки кластерів?

Це найважча річ щодо k-засобів, і існує маса методів. Давайте скористаємося ліктьовим методом

#plot variance for each value for 'k' between 1,10
initial = [cluster.vq.kmeans(tests,i) for i in range(1,10)]
pyplot.plot([var for (cent,var) in initial])
pyplot.show()

Локоть сюжет

Призначте свої спостереження класам та побудуйте їх

Я вважаю, що індекс 3 (тобто 4 кластери) настільки ж хороший, як і будь-який

cent, var = initial[3]
#use vq() to get as assignment for each obs.
assignment,cdist = cluster.vq.vq(tests,cent)
pyplot.scatter(tests[:,0], tests[:,1], c=assignment)
pyplot.show()

Діаграма розкиду

Просто опрацюйте, куди ви можете вставити все, що ви вже зробили, у цей робочий процес (і я сподіваюся, що ви кластери трохи приємніше, ніж випадкові!)


Ваша відповідь виглядає чудово. Чи можу я це ефективно використовувати для своїх даних. Неможливо отримати час, щоб спробувати.
користувач2721

@david w: Це одна з найкращих відповідей, які я бачив! Дуже дякую за публікацію окремого прикладу. Принаймні, я розумію суть вашої відповіді :) Ще раз дякую!
Легенда

@david w: Єдине питання, який у мене був, - метод ліктьового суглоба, показує зростаючі значення, а ваш сюжет показує зменшення. Це тому, що ви використовуєте значення спотворень безпосередньо з kmeans? Як я можу конвертувати це так, щоб він виглядав як ліктя Вікіпедії? І як останнє запитання, чи не знаєте ви, як це зробити для kmeans2 замість kmeans?
Легенда

Ви виявили індекс 3, просто дивлячись на перший сюжет?
Сигур

2

Можливо, спробуйте щось на кшталт Fastmap щоб побудувати свій набір знаків, використовуючи відносні відстані.

(досі) нічого розумного не написав Fastmap в python для створення рядків і може бути легко оновлений для обробки списків атрибутів, якщо ви записали власну метрику відстані.

Нижче наведена стандартна евклідова відстань, яку я використовую, яка приймає два списки атрибутів як параметри. Якщо ваші списки мають значення класу, не використовуйте його для обчислення відстані.

def distance(vecone, vectwo, d=0.0):
    for i in range(len(vecone)):
        if isnumeric(vecone[i]):
            d = d + (vecone[i] - vectwo[i])**2
        elif vecone[i] is not vectwo[i]:
            d += 1.0
    return math.sqrt(d)

def isnumeric(s):
    try:
        float(s)
        return True
    except ValueError:
        return False

0

Я не фахівець з пітонів, але надзвичайно корисно побудувати основні компоненти 1-го другого один проти одного на осях x, y.

Не впевнені, які пакунки ви використовуєте, але ось приклад посилання:

http://pyrorobotics.org/?page=PyroModuleAnalysis


Я не є експертом зі статистики. Не могли б ви пояснити більше про сюжетну ідею?
користувач2721

Основна ідея полягає в тому, що багато змінних співвідносяться між собою, і все можна звести лише до двох змінних, які не співвідносяться між собою і пояснюють "більшість" варіацій даних. Вам потрібно ознайомитися з аналізом основних компонентів і застосувати пакет, який дозволяє реалізувати його. en.wikipedia.org/wiki/Principal_component_analysis
Ральф Вінтерс
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.