Блог

Главная - Вопросы по программированию - Возможно ли изменить порядок столбцов в heatmaply?

Возможно ли изменить порядок столбцов в heatmaply?

#r #cluster-analysis #heatmap #dendextend #heatmaply

#r #кластерный анализ #тепловая карта #дендекстенд #тепловая карта

Вопрос:

Я создал дендрограмму с помощью heatmaply с помощью объекта dendrogram из пакета dendextend. Мне было интересно, есть ли способ изменить порядок столбцов: vs, am, carb, wt, drat, gear и т.д. Например, я хочу отсортировать их в алфавитном порядке. Возможно ли это?

 library(heatmaply)
library(tidyverse)
library(RColorBrewer)
library(dendextend)

custom_colors <- c("#B6E2D3", "#887BB0", "#FFD53D", "#EF7C8E", "#A1F7A1", "#EFE7BC", "#C26DBC", "#14C2DD",
                   "#5CD85A", "#AB6B51", "#39918C", "#FFA384", "#57DDF3", "#D0B49F", "#EEB5EB", "#74BDCB")

dend_example <- mtcars %>% dist %>% hclust(method = "ward.D2") %>% 
  as.dendrogram %>% color_branches(k=5) %>% color_labels

heatmaply_mtcars<- heatmaply(mtcars, hclustfun = hclust, hclust_method = "ward.D2",
                             Rowv = dend_example, 
                             main = "mtcars",
                             show_dendrogram = c(TRUE, FALSE),
                             showticklabels = c(TRUE, TRUE),
                             scale_fill_gradient_fun = ggplot2::scale_fill_gradient2(high = "#025492"),
                             file = "heatmaply_mtcars.html",
                             column_text_angle = 30,
                             colorbar_thickness = 50)

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

Ответ №1:

Я включил дендрограмму столбца в heatmaply с аргументом colv. Дендрограмма столбца была повернута с помощью функции rotate из пакета dendextend:

 dend_example <- mtcars %>% dist %>% hclust(method = "ward.D2") %>% 
  as.dendrogram %>% color_branches(k=5) %>% color_labels

column_dend <- t(mtcars) %>% dist %>% hclust(method = "ward.D2") %>% 
  as.dendrogram
column_dend  <- column_dend %>% rotate(c(vector_with_labels))

heatmaply_mtcars<- heatmaply(mtcars, hclustfun = hclust, hclust_method = "ward.D2",
                             Rowv = dend_example, 
                             Colv = column_dend,
                             main = "mtcars",
                             show_dendrogram = c(TRUE, FALSE),
                             showticklabels = c(TRUE, TRUE),
                             scale_fill_gradient_fun = ggplot2::scale_fill_gradient2(high = "#025492"),
                             file = "heatmaply_mtcars.html",
                             column_text_angle = 30,
                             colorbar_thickness = 50)
Метки: