前面两篇文章介绍到pheatmap返回行和列的hclust结果:
1、按照前两篇文章的基础数据获得的热图结果,对行(基因的聚类结果)聚类树单独进行绘制,可以借助plot()函数:
默认画法
> plot(list$tree_row)
而有一个有趣的函数rect.hclust()可以将hclust结果绘图,并圈选出不同的cluster,譬如设定3个cluser:
> rect.hclust(list$tree_row,k=3)
得到如下结果:
2、而将hclust结果转换成另类树形图as.dendrogram()可以选择type画出其他形状的树形图:
这个树形图一共包含两个type,可以查询?as.dendrogram() 或者?plot.dendrogram得知
> plot(as.dendrogram(list$tree_row),type = "triangle")