0

我有一个马氏距离的邻近矩阵（相异）。

矩阵（样本）：

> dput(MD[1:5,1:5])

structure(c(0, 10.277, 8.552, 8.592, 9.059, 10.277, 0, 10.917, 
9.489, 8.176, 8.552, 10.917, 0, 8.491, 8.104, 8.592, 9.489, 8.491, 
0, 9.375, 9.059, 8.176, 8.104, 9.375, 0), .Dim = c(5L, 5L), .Dimnames = list(
    c("2", "4", "5", "6", "9"), c("X2", "X4", "X5", "X6", "X9"
    )))

矩阵有 1900 人，行名是一个 ID。我需要对这些人进行聚类，并在该人的 ID 旁边获得一些聚类。

我知道如何使用 k-means 进行聚类，但是当您已经有一个相异矩阵时，我不知道如何进行聚类。

4

2 回答 2

0

您可以使用层次聚类，从马氏距离矩阵开始：

MD
      X2     X4     X5    X6    X9
#2  0.000 10.277  8.552 8.592 9.059
#4 10.277  0.000 10.917 9.489 8.176
#5  8.552 10.917  0.000 8.491 8.104
#6  8.592  9.489  8.491 0.000 9.375
#9  9.059  8.176  8.104 9.375 0.000

hc <- hclust(as.dist(MD))

clusters <- cutree(hc, k = 3) # obtain 3 clusters
clusters
#2 4 5 6 9 
#1 2 3 1 3 

plot(hc)
rect.hclust(hc, k = 3, border = "red")

在此处输入图像描述

于 2017-02-13T17:53:55.400 回答
0

如果我理解正确，请使用dendextend：

fit<-hclust(MD)

fit %>% as.dendrogram %>% 
  set("branches_k_color", k = 3, value = c("purple", "orange","red")) %>% 
  plot

簇 # 由颜色指定......从 . 传递过来cutree。

于 2017-02-13T18:39:52.630 回答