Как я могу показать названия образцов / генотипов для каждого кластера? - PullRequest
0 голосов
/ 10 октября 2018

Мой вопрос: как я могу показать названия образцов для каждого кластера?

пример:

Cluster 1: X20, X45, ..., X10 
Cluster2: X104, X50,..., X33

## создание кластеров hclust

## получить 10 групп из выборки

 out1 <- cutree(hc1, k = 50)

## нарисовать дендограмму с красными границами вокруг 5 кластеров

rect.hclust(hc1, k=50, border="red")
save(out1, file = "out1.rda") 
load("out1.rda")

## количество выборок в каждом кластере

> table(out1)

out1

 1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19  20  21  22  23 
 10  22 108 181 122  21  12  28   5  75  27  14  39   4   5   4  11   9  29   2   5   4  14 
 24  25  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46 
  4   2  22  23   3  10   1   6  11   8   2   5   9   6   4   5   1   1   1   2   9  14   1 
 47  48  49  50 
  4   1   1   1 

Дендрограмма кластера, показывающая 50 кластеров (щелкните здесь)

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...