pheatmap返回的结果是什么

发布时间:2022-03-18 17:35:37 作者:iii
来源:亿速云 阅读:231

这篇文章主要介绍了pheatmap返回的结果是什么的相关知识,内容详细易懂,操作简单快捷,具有一定借鉴价值,相信大家阅读完这篇pheatmap返回的结果是什么文章都会有所收获,下面我们一起来看看吧。

pheatmap是一个热图绘制的R包,全称pretty heatmap。利用的绘图函数是pheatmap(),对应的数据则是一个数值矩阵,譬如基于如下的矩阵10X6:

> mat
              CK-WT-1      CK-WT-2      CK-WT-3   CK-tdr1-1   CK-tdr1-2   CK-tdr1-3
AT1G01010    3.741490    7.3618000    5.8173400   5.7113100   7.9705400  10.3762000
AT1G01030    5.235280    2.7707000    2.6685900   3.2263200   1.3210500   1.9672600
AT1G01040    2.821317    1.5633947    1.9316282   3.1948090   2.6008540   2.3012776
AT1G01050  118.660000  140.1430000  123.3830000  97.2229000  95.2539000  91.8525000
AT1G01060    1.873769   16.9090246    0.9559375   0.4774184   0.5273923   0.4333881
AT1G01070    1.710346    0.7802436    2.7996091   4.7297117   4.3637146   3.3732144
AT1G01080   30.261105   33.6571056   31.4669073  30.9689028  28.0334025  27.6215018
AT1G01090   86.640500  116.4700000  111.7380000  82.8809000  85.7292000  79.6027000
AT1G01100 1448.627845 1634.0417300 1566.5625600 648.4734402 683.4113500 647.2749880
AT1G01120  173.996000   51.0019000   52.3322000  66.9486000  41.1148000  49.7169000

直接按照默认参数绘图获得的结果是返回一张热图:

library(pheatmap)
pheatmap(mat,scale = "row")

而实际是pheatmap()返回的结果涉及多个:

Value

Invisibly a pheatmap object that is a list with components


根据以上的解释,可以知道pheatmap返回的对象是列表,其中包含了多个聚类的结果,其中一个组成结果为行聚类结果,对应的是行基于hclust聚类返回的对象,而hclust返回的对象也是一个列表,包括了多组信息。

以上面的代码为例,将结果赋给一个新的向量,然后查看向量的结果:

list=pheatmap(mat,scale="row")
> summary(list)
         Length Class  Mode   
tree_row 7      hclust list   
tree_col 7      hclust list   
kmeans   1      -none- logical
gtable   6      gtable list

可以看到list$tree_row也是一个列表:

> summary(list$tree_row)
            Length Class  Mode     
merge       18     -none- numeric  
height       9     -none- numeric  
order       10     -none- numeric  
labels      10     -none- character
method       1     -none- character
call         3     -none- call     
dist.method  1     -none- character 
> list$tree_row$order
 [1]  2 10  5  7  8  4  9  1  3  6

而其中有一个重要的内容是order,他是进行行聚类重排后的排列顺序,基于这个信息可以输出排序后的矩阵,从而获得与热图显示顺序一致的文件结果:

> newOrder=mat[list$tree_row$order,]
> newOrder
              CK-WT-1      CK-WT-2      CK-WT-3   CK-tdr1-1   CK-tdr1-2   CK-tdr1-3
AT1G01030    5.235280    2.7707000    2.6685900   3.2263200   1.3210500   1.9672600
AT1G01120  173.996000   51.0019000   52.3322000  66.9486000  41.1148000  49.7169000
AT1G01060    1.873769   16.9090246    0.9559375   0.4774184   0.5273923   0.4333881
AT1G01080   30.261105   33.6571056   31.4669073  30.9689028  28.0334025  27.6215018
AT1G01090   86.640500  116.4700000  111.7380000  82.8809000  85.7292000  79.6027000
AT1G01050  118.660000  140.1430000  123.3830000  97.2229000  95.2539000  91.8525000
AT1G01100 1448.627845 1634.0417300 1566.5625600 648.4734402 683.4113500 647.2749880
AT1G01010    3.741490    7.3618000    5.8173400   5.7113100   7.9705400  10.3762000
AT1G01040    2.821317    1.5633947    1.9316282   3.1948090   2.6008540   2.3012776
AT1G01070    1.710346    0.7802436    2.7996091   4.7297117   4.3637146   3.3732144

关于“pheatmap返回的结果是什么”这篇文章的内容就介绍到这里,感谢各位的阅读!相信大家对“pheatmap返回的结果是什么”知识都有一定的了解,大家如果还想学习更多知识,欢迎关注亿速云行业资讯频道。

推荐阅读:
  1. 易语言中的返回命令是什么
  2. 浅谈Java中的this作为返回值时返回的是什么

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

pheatmap

上一篇:blastall的参数有哪些

下一篇:awk匹配怎么替换gsub

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》