2016-11-07 115 views
0

我有一個熱圖,看起來像這樣在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") 

enter image description here

但後來我做了基於相關的單獨聚類分析,說出來了這樣的:

library(Hmisc) 
plot(varclus(qq,similarity="spearman")) 

enter image description here

哪有我修改了我的熱圖,以便聚類與我關聯的聚類分析相同?我需要以某種方式修改heatmap.2函數(或可能使用不同的函數)以基於皮爾森相關性。有任何想法嗎?

+1

確實添加'distfun =函數(x)的as.dist(1 - COR(T(X ),method ='sp')^ 2)'到你的''heatmap.2'' get it> – rawr

回答

0

嘗試

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不存在,我不能再現圖像)