Тепловая карта дендрограммы на основе корреляции в R
У меня есть тепловая карта, которая выглядит так в R:
col<- colorRampPalette(c("red","white", "blue"))(10)
library("gplots")
heatmap.2(qq,scale="none",col=col,trace="none",density.info="none",dendrogram="column")
Но затем я сделал отдельный кластерный анализ на основе корреляции, который получился следующим образом:
library(Hmisc)
plot(varclus(qq,similarity="spearman"))
Как я могу изменить свою тепловую карту так, чтобы кластеризация была идентична кластерному анализу, который я сделал с корреляцией? Мне нужно как-то изменить heatmap.2
функция (или, возможно, использовать другую функцию), чтобы основываться на корреляции Пирсона. Есть идеи?
1 ответ
Решение
Пытаться
col<- colorRampPalette(c("red","white", "blue"))(10)
library("gplots")
library(Hmisc)
v <- varclus(qq,similarity="spearman")
devtools::install_github('talgalili/dendextend')
library(dendextend)
dend <- as.dendrogram(v) # comes from dendextend. The same as as.dendrogram(v$hclust)
heatmap.2(qq,scale="none",col=col,trace="none",density.info="none",dendrogram="column", Colv = dend)
(поскольку qq не существует, я не могу воспроизвести изображение)