相关性分析热点图_高分文章中物种与代谢物相关性热图是怎么画的?
图1 物种代谢物热图(2015,Cell Host& Microbe,IF=15.753 )[1]
图2 物种代谢物热图(2018,NatureMedicine,IF=30.641)[2]那么,该如何画出此类高分文章中的相关性热图呢?这里,以16s与代谢组的数据为例,向大家分享如何使用R语言进行两个组学数据的相关性计算、绘制相关性热图。1.加载R包
#psych包用于计算相关性、p值等信息library(psych)#pheatmap包用于绘制相关性热图library(pheatmap)#reshape2包用于输出数据的整合处理library(reshape2)
2.读入数据
#读取微生物丰度信息表#表头需带有分类水平、物种名称等关键信息#第一列为样本名称信息phy "phy.xls", sep = "\t", header = T,row.names=1)
图3 微生物丰度信息表格
#读取代谢物信息表met "met.xls", sep = "\t", header = T,row.names=1)
图4 代谢物丰度信息表格3.计算相关性、p值
#计算相关性矩阵(可选:”pearson”、”spearman”、”kendall”相关系数)、p值矩阵cor "pearson",adjust="none")#提取相关性、p值cmt pmt head(cmt)head(pmt)
4.数据保存
#输出相关系数表格,第一行为代谢物信息,第一列为物种信息cmt.outwrite.table(cmt.out,file="cor.txt",sep="\t",row.names=F)
图5 相关性系数表格
#输出p值表格,第一行为代谢物信息,第一列为物种信息pmt.outwrite.table(pmt.out,file="pvalue.txt",sep="\t",row.names=F)
图6 p值表格
#以关系对的形式输出表格#第一列为物种名,第二列为代谢物名,第三、第四列对应显示相关系数与p值df "cor")df$pvalue as.vector(pmt)head(df)write.table(df,file="cor-p.txt",sep="\t")
图7 关系对信息5.绘制显著性标记
0.01的以“**”标注,p值
6.绘制相关性热图
#自定义颜色范围
图8 R语言绘制的物种+代谢物相关性热图
#图片保存,代码中输入”filename=”,或在R语言软件中点击“文件-另存为” 进行保存pheatmap(cmt,scale = "none",cluster_row = T, cluster_col = T, border=NA, display_numbers = pmt, fontsize_number = 12, number_color ="white", cellwidth = 20, cellheight =20,color=mycol,filename="heatmap.pdf")
参考文献[1]Kostic AD, Gevers D, Siljander H, et al. The dynamics ofthe human infant gut microbiome in development and in progression toward type 1diabetes. Cell Host Microbe. 2015;17(2):260–273.doi:10.1016/j.chom.2015.01.001[2]Hoyles, Lesleyet al. “Molecular phenomics and metagenomics of hepatic steatosis innon-diabetic obese women.” Nature medicine vol. 24,7 (2018):1070-1080. doi:10.1038/s41591-018-0061-3
实用科研工具推荐
详实生信软件教程分享
前沿创新组学文章解读
独家生信视频教程发布
相关性分析热点图_高分文章中物种与代谢物相关性热图是怎么画的?相关推荐
- java 线图_在Java中绘制一个简单的线图
您的代码和建议的问题: >再次,您需要更改组件的preferredSize(这里是Graph JPanel),而不是大小 >不要设置JFrame的边界. >在添加组件之后并在调用se ...
- r语言聚类分析_技术贴 | R语言pheatmap聚类分析和热图
点击蓝字↑↑↑"微生态",轻松关注不迷路 本文由阿童木根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 pheatmap默认会对输入矩阵数据的行和列同时进行聚 ...
- 服务器画热图显示无效的变量,使用pheatmap包绘制热图
加载所需R包 library(pheatmap) 设置工作路径 setwd("/Users/Davey/Desktop/VennDiagram/") # 清除当前环境中的变量 rm ...
- 34.35.热图(heatmap)、创建带注释的热图、使用辅助函数的代码样式、图像显示、图像插值、将图像数据导入Numpy数组、将numpy数组绘制为图像
34.热图(heatmap) 34.1.创建带注释的热图 34.2.使用辅助函数的代码样式 35.图像显示 35.1.图像插值 35.2.将图像数据导入Numpy数组 35.3.将numpy数组绘制为 ...
- origin做相关性分析图_高分文章的相关性热图,十分钟教你学会无代码复现
相关性热图可以便捷的显示多个变量之间的相关性,因此在组学数据分析中得到广泛应用,例如基因的共表达分析.样本重复检验.微生物群落的共发生网络分析等. 在昨天推文下就有人问热图怎么画?画热图最常规的办法自 ...
- 相关性分析热点图_车联网数据分析(二):驾驶行为分析
出门在外,遇上交通拥堵.事故,那可真是头大.许多交通拥堵和交通事故与驾驶员的驾驶行为有关.驾驶行为通常也是驾驶习惯,习惯成自然.在开车的过程中,如果有不良的驾驶行为,驾驶员自己很难意识到,但是后果可能 ...
- python 加载动图_在浏览器中使用TensorFlow.js和Python构建机器学习模型(附代码)...
大数据文摘授权转载自数据派THU 作者:MOHD SANAD ZAKI RIZVI 本文主要介绍了: TensorFlow.js (deeplearn.js)使我们能够在浏览器中构建机器学习和深度学习 ...
- 相关性分析p值_一行代码掌握皮尔逊相关分析,洞察变量关系
变量类型与推荐的假设检验方法 可以看到,当我们探索两个连续变量之间的关系时,相关分析是一个很好的选择.那么,相关分析的原理是什么?如何在Python中实现相关分析呢? 一.Pearson相关系数 针对 ...
- spss相关性分析看结果_利用spss做Pearson相关性分析步骤详解
有蛮多的学生私信老徐问如何利用spss做相关性分析,其实相关性分析应该是spss分析中较为基础的一个功能应用,很多学生可能是因为跨专业或者对统计软件了解较少,在没有经过系统学习的前提下,感觉云里雾里. ...
最新文章
- 如何在nuget上传自己的包+搭建自己公司的NuGet服务器(新方法)
- 特斯拉前员工承认备份AutoPilot代码,但入职小鹏汽车前已删除
- 苹果CMS V10 播放记录_苹果cms新手入门安装配置教程
- firefox自动化测试的常用插件
- 内核ko模块strip使用
- ceph 分布式存储安装
- gitglone 指定分支_安装说明 · harryxu/gitube Wiki · GitHub
- shiro服务器取不到用户信息,SpringMvc+Shiro有时候shiro取不到用户名
- Python数据科学学习进阶
- java 链接kafka单机版_kafka单机环境搭建及其基本使用
- Kepware OPC与服务器通讯
- [区块链]对话区块链人工智能顶级实力公司ObEN
- e的ax次方怎么求导
- 集群断电重启The connection to the server apiserver.k8s:6443 was refused - did you specify the right
- 笔试——分频电路设计
- flask_mail通过qq发送邮箱
- ubuntu安装chrome浏览器64位
- clear both
- seleniumbase学习总结4 - 运行测试用例
- sql 除法中如何防止被除数为0,以及SUM和舍入为指定的小数位数函数