Цветные метки ggdendro plot

Question

Я пытаюсь применить dendextend, чтобы раскрасить ярлыки ggdendro.Я получаю эту ошибку:

Error in UseMethod("as.dendrogram") : 
  no applicable method for 'as.dendrogram' applied to an object of class "dendro"

Почему конвертация не работает?Есть идеи, как это решить?

Код:

library(IncDTW)
library(ggplot2)
library(ggdendro) 
library(dendextend)

A <- matrix(1:50, nrow = 50, ncol = 1)
B <- matrix(1:75, nrow = 75, ncol = 1)
C <- matrix(25:49, nrow = 25, ncol = 1)
D <- matrix(1:50, nrow = 50, ncol = 1)
treeList <- list(A,B,C,D)

names(treeList)[1] <- "A"
names(treeList)[2] <- "B"
names(treeList)[3] <- "C"
names(treeList)[4] <- "D"

result <- dtw_dismat(treeList, dist_method = "norm2", return_matrix = F)
distMatrixResult <- result$dismat
hc <- hclust(distMatrixResult, method = "average")
hc <- dendro_data(hc)
dict <- setNames(names(treeList), 1:length(names(treeList)))
hc$labels$label <- sapply(hc$labels$label, function(x) dict[[as.character(x)]])
ggdendrogram(hc, rotate = TRUE, theme_dendro = FALSE)
dend <- as.dendrogram(hc)
d5=color_branches(dend)
plot(d5)
d5g=color_branches(dend,5,groupLabels=TRUE)
plot(d5g)
d5gr=color_branches(dend,5,groupLabels=as.roman)
plot(d5gr)

Answer 1

В текущей версии IncDTW вам нужно добавить имена после использования dtw_dismat (), в будущей версии я это исправлю.А пока попробуйте:

result <- dtw_dismat(treeList, dist_method = "norm2", return_matrix = TRUE)
distMatrixResult <- result$dismat
dnames <- LETTERS[1:4]
dimnames(distMatrixResult) <- list(dnames , dnames)
distMatrixResult <- as.dist(distMatrixResult)

также взгляните на этот пост: ggdendrogram label

Answer 2

Кажется, вы путали дендрограмму с данными дендрограммы.Это будет работать:

dend <- as.dendrogram(hc)

dict <- setNames(names(treeList), 1:length(names(treeList)))
hc$labels$label <- sapply(hc$labels$label, function(x) dict[[as.character(x)]])
ggdendrogram(dend, rotate = TRUE, theme_dendro = FALSE)

d5=color_branches(dend)
plot(d5)
d5g=color_branches(dend,4,groupLabels=TRUE)
plot(d5g)
d5gr=color_branches(dend,4,groupLabels=as.roman)
plot(d5gr)

Результат:

Редактировать:

Поскольку @ maxatSOflow указал , мы можем использовать usedist::dist_setNames(), чтобы установить имена для объекта "dist".Итак, чтобы, наконец, пометить дендрограмму как положено, вы можете сделать следующее:

library(IncDTW)
result <- dtw_dismat(treeList, dist_method="norm2", return_matrix=FALSE)
distMatrixResult <- result$dismat

# install.packages("usedist")
library(usedist)
distMatrixResult <- dist_setNames(distMatrixResult, LETTERS[1:4])

hc <- hclust(distMatrixResult, method="average")

library(ggdendro) 
ggdendrogram(hc, rotate=TRUE, theme_dendro=FALSE)
dend <- as.dendrogram(hc)


library(dendextend)
d5=color_branches(dend)
plot(d5)
d5g=color_branches(dend,4,groupLabels=TRUE)
plot(d5g)
d5gr=color_branches(dend,4,groupLabels=as.roman)
plot(d5gr)

Дает:

Данные:

A <- matrix(1:50, nrow=50, ncol=1)
B <- matrix(1:75, nrow=75, ncol=1)
C <- matrix(25:49, nrow=25, ncol=1)
D <- matrix(1:50, nrow=50, ncol=1)

treeList <- setNames(list(A, B, C, D), LETTERS[1:4])