2013-12-10 27 views
0

我产生的热图基于一个4500行×400山口矩阵m:提取热图的特定成员中的R

hm<-heatmap(m, na.rm=TRUE, col=bluered(256)) 

那么,如何看到特定集群的成员?我读了一下,知道我必须使用'cutree()',但我仍然不确定如何正确实现它。

谢谢。

回答

1
hm<-heatmap.2(m, na.rm=TRUE) 
clstrs <- cutree(as.hclust(hm$rowDendrogram), 1:dim(m)[1]) 

并根据您需要多少个群集,选择正确的列。因此,如果你想要7个簇:

clstrs[, 7] 

会告诉你每个元素属于哪个簇。

+0

我明白了。谢谢你的答案。这是否仅在heatmap.2()中起作用,或者我是否也可以将它应用于heatmap()? –

+0

我认为它只适用于heatmap.2;如果你不想安装任何新的软件包,你可以试试 hcls < - hclust(dist(m)); cutree(hcls,k = 6:7) –