Hmisc包

推荐一个Hmisc包

下载好后使用#res=rcorr(as.matrix(rt))#这个函数可将rt列表中的任意两个基因展示相关系数与P值

cormatrix=function(cor,p){

  ut=upper.tri(cor)

  data.farme(row=rownames((cor)[row(cor)[urt],]

                          colum=rownames(cor)[col(cor)[ut]],

                          cor=(cor)[ut],

                          p=p[ut]))

}

cormatrix(res$r,res$p)

将数据以表格形式呈现

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容