输出聚类相似度矩阵

时间:2015-04-11 23:28:57

标签: r heatmap hierarchical-clustering

我已经生成了一个皮尔逊相似度矩阵,并使用pheatmap绘制了结果(使用hclust进行聚类,方法="完成")。我想输出有序矩阵,但在R中,默认似乎只是按字母顺序排列所有内容。

这是我的代码:

df <- cor(t(genes), method = "pearson")  
pheatmap(df, clustering_method = "complete")

头(基因)

       pre           early          mid            late       end
AAC1   2.0059007     3.64679740     3.0092533      2.4936171  2.2693034
AAC3  -1.6843969    -1.62572636    -0.7654462     -1.5827481 -1.6059080
AAD10  2.6012529     2.05759631     1.3665322      1.4590833  0.3778324
AAD14  0.5047704     0.76021375     0.1825944      0.6111774  0.1174208
AAD15  7.6017557     8.52315453     7.2605744      6.9029452  5.9028824
AAD16  1.2018193    -0.03285354     0.2229450     -0.1337033  0.2198542

这是当前输出(df)的样子:

     A     B     C    D  
A    1   0.5  0.25  0.1  
B  0.1     1   0.1  0.5    
C  0.5   0.2     1  0.2  
D    0   0.1   0.7    1

如何按hclust的顺序输出相似度矩阵?

我看过了,但我还没能找到任何可以完成我需要的东西。在此先感谢您的帮助!

(也很抱歉,我不知道如何正确格式化所有内容)

编辑:也许一些视觉效果会有所帮助。我的聚集pheatmap输出如下所示:ordered heatmap

我可以看到表现相似的基因组,但因为有很多基因不可能/没用来阅读标签。我想找出哪些基因聚集在一起,但我无法输出有序矩阵。

当我绘制没有聚类的数据时,它看起来像这样:unclustered heatmap

因此,我可以获得的输出/数据对于进一步分析几乎没用。

0 个答案:

没有答案