使用corrplot绘制行、列不同,且带有p值显著性标注的相关系数图

news2024/11/15 10:17:02

导读:

相关系数衡量两个变量之间的线性关系,通常以N*N的矩阵形式展示。例如样品vs样品,或者基因vs基因的相关性。本文介绍了使用corrplot R包绘制M*N的相关系数矩阵,例如M个基因表达与N个代谢物信号间的相关性,同时带上p值显著性标注。

Diabetes, Metabolic Syndrome and Obesity文章《Dysbiosis of the Salivary Microbiome is Associated with Hypertension and Correlated with Metabolic Syndrome Biomarkers》fig 9 展示了不同微生物与生化、表型特征间的相关系数。图中X轴为不同的微生物,Y轴为不同的血浆生化指标或人口统计特征。图中的点表示spearman相关系数,红色点表示负相关,绿色点表示正相关。红色点越大,表示负相关越大;绿色点越大,表示正相关越大。点上的星号(*)表示p值。***: p<=0.001, **: p<=0.01,*:p<=0.05。从图中可以看出,收缩压(Systolic BP)、舒张压(Diastolic BP)、葡萄糖(Glucose)等指标与众多微生物正相关,且显著。而乳制品摄入(Dairy_intake)、胰岛素(insulin)、吸烟(Smoking)等与微生物相关性不显著。

相关系数简介

相关系数(correlation coefficient)是一个介于-1和1之间的数,用于衡量两个变量之间的线性关系。如果两个变量之间的变化趋势是一致的,即一个变量增加时另一个也增加,那么它们之间的相关系数为正;如果一个变量增加而另一个减少,则相关系数为负。接近1的相关系数表示强正相关。接近-1的相关系数表示强负相关。接近0的相关系数表示两个变量之间没有或几乎没有线性关系。

常见的相关系数包括:

  • 皮尔逊相关系数(Pearson correlation coefficient):适用于测量两个连续变量之间的线性相关性
  • 斯皮尔曼等级相关系数(Spearman’s rank correlation coefficient):适用于测量两个变量的等级(或排序)之间的相关性,不要求数据是线性的

Spearman相关系数更适用于非线性关系和不满足正态分布的数据,适用于评估变量之间的顺序关系;而Pearson相关系数适用于线性相关性和正态分布的数据,用于评估两个连续变量之间的线性关系程度。

注意:相关系数只能告诉我们变量之间是否存在线性关系,以及这种关系的强度和方向,但它不能告诉我们这种关系的原因或机制。此外,相关系数对异常值比较敏感。

1打开作图URL

微生信-在线绘制带p值显著性标注的相关系数图,支持行列不同

2,示例数据

点击图片上方的示例数据,下载,并使用excel打开。

示例数据包括2个矩阵:

矩阵1:mRNA表达矩阵。行是样品,列是基因。交叉的cell为标准化表达值。

矩阵2:代谢物信号强度矩阵。行是样品(与表达矩阵完全一样,包括顺序),列是代谢物。交叉的cell为标准化信号强度。

注意:这里计算的是矩阵1的列元素和矩阵2的列元素间的相关性。它们的列数可以相同(N*N矩阵),也可以不同(N*M矩阵)。当然,两个矩阵也可以完全一样(即最常见的相关系数矩阵)。

3,输入检查

示例数据:点击输入框下面的“示例”按钮,将载入示例数据。

真实数据:数据放在excel中,调整好后,Ctrl+A选中数据,Ctrl+C拷贝,Ctrl+V粘贴数据到输入框中。

然后使用输入框下面的“输入检查”按钮先对输入数据进行检查。若检查不通过,请根据检查提示重复【修改-输入检查】步骤,直到检查通过(如下图所示),然后可以继续选择参数。

注:输入检查是新加功能,它会根据不同模块的输入要求,逐行逐列检查输入数据,并给出提示,确保数据符合模块要求。

4,选择参数

图片大小:包括图片宽度,图片高度。图片不易过小,太小的话,文字有可能被截断,此时可以使用acrobat illustrator或者inkscape软件编辑pdf/svg文件调整

相关系数计算方法:Pearson或者Spearman

图形:绘图时每个小矩形里边的图形形状。包括圆形,正方形,椭圆形,阴影(负的值加阴影),颜色(完全填充的颜色矩形),饼图(见文档最后)。

布局:矩阵形式,上三角或者下三角(见文档最后)

颜色:低、中、高三种颜色。若低=中,或者中=高时,仅两种颜色

Colorbar数值范围:留空则默认。填写的话则根据所填数值限定colorbar的范围,但是所填数值必须完全覆盖原始范围,否则还是默认。例如若实际范围是-0.4到0.5,那么可以填写-0.5,0.5,但是填写-0.4,0.4则还是会自动转成-0.4到0.5。

Colorbar的字体大小

行列名字:字体大小,颜色

显著性标注:默认是0.001,0.01,0.05,对应***:p<=0.001,**:p<=0.01,*:p<=0.05

显著性标注字符:默认是*,也可以使用其他符号,例如+

显著性标注字体大小,颜色

是否根据相关性排序:默认按照输入的顺序绘图,若要根据相关性进行聚类(根据相关性进行层次聚类),则可以选择排序,并可以选择添加黑色矩形边框,以达到突出显示的效果。必需是1-4。1就是整个矩阵添加一个边框。注意:排序选项仅支持N*N的矩阵(见文档最后)

是否显示对角线:针对N*N矩阵,可以选择不显示对角线都是1的相关系数。

字体:设置了期刊杂志中最常用的两种:Times New Roman和Arial。如需使用其他字体,可以使用acrobat illustrator编辑生成的pdf或者svg图片

5,提交出图

检查通过,并且参数选好后,点击“提交”按钮,约2s后,会在页面上显示带p值显著性标注的相关系数矩阵图。我们提供了pdf、svg两种矢量图,png、tiff两种标量图供大家下载使用。可以使用acrobat illustrator等软件编辑矢量图,进行组图,调整字体等操作,以满足论文要求。

同时提供了计算结果,相关系数矩阵,p值矩阵。

也可以绘制其他样式的图。例如:

微生信助力高分文章,用户200000+,谷歌学术3900

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2094217.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

国产芯片+国产操作系统打造办公系统

在《使用国产操作系统作为开发系统》一文中&#xff0c;我介绍了将开发系统从 Ubuntu 替换为 Deepin 系统的过程。经过一个多月的使用&#xff0c;Deepin 系统已然成为我的主力开发平台&#xff0c;其顺手程度让我对国产操作系统的信心大增。于是&#xff0c;我开始将目光瞄向公…

顶级开源许可证详解

目录 软件许可证类型&#xff1a;版权左派和宽容型 顶级开源许可证详解 GNU 通用公共许可证 (GPL) Apache 许可证 Microsoft 公共许可证 (Ms-PL) 伯克利软件发行版 (BSD) 通用开发和分发许可证 (CDDL) Eclipse 公共许可证 (EPL) MIT 许可证 了解你的开源许可证&#…

java编辑器——IntelliJ IDEA

java编辑器有两种选择——IntelliJ IDEA和VsCode。其中IntelliJ IDEA现在是企业用的比较多的&#xff0c;是专门为java设计的&#xff0c;而VsCode则是通过插件来实现Java编辑的。 1.IntelliJ IDEA 官网下载链接&#xff1a;https://www.jetbrains.com/idea/ 注意选择社区版…

AWS-亚马逊网络服务(基础服务)-AWS 定价计算器-概述与动手部署:

让我们来概述并亲身实践如何使用 AWS 定价计算器来计算 概述&#xff1a; AWS 定价计算器是 Amazon Web Services (AWS) 提供的基于 Web 的工具&#xff0c;可帮助用户估算其特定用例的 AWS 服务成本。欢迎来到雲闪世界。 它允许客户建模他们的基础设施并根据他们打算使用的…

【AI 绘画】更快?更省显存?支持 FLUX?使用绘世启动器安装 SD WebUI Forge

使用绘世启动器安装 SD WebUI Forge 下载绘世启动器 绘世启动器下载地址1&#xff1a;https://gitee.com/licyk/term-sd/releases/download/archive/hanamizuki.exe 绘世启动器下载地址2&#xff1a;https://www.bilibili.com/video/BV1ne4y1V7QU 新建一个文件夹取名sd-webui-…

中仕公考怎么样?公务员考试什么时候补录?

公务员考试补录的时间和方法通常因地区和职位的不同有所区别&#xff0c;一般来说&#xff0c;这一过程会在面试、体检和考核环节完成后启动。 如果在招录过程中出现职位空缺或者并未全部招满的情况&#xff0c;就会进行补录。用人单位会通过其官方或公告形式公布相关信息&…

指针5.回调函数与qsort

今天来学习回调函数与qsort 目录 1.回调函数实现模拟计算器代码的简化原代码运行结果简化代码运行结果 qsort函数排序整型数据代码运行结果 qsort排序结构数据代码 qsort函数的模拟实现代码运行结果 总结 1.回调函数 回调函数就是⼀个通过函数指针调用的函数。 如果你把函数的…

JavaEE第22节 TCP段(报文)结构剖析

目录&#xff08;关于字段有不理解的&#xff0c;哪里不会点哪里&#x1f618;&#xff09; 逻辑结构字段解析一、源端口&目的端口二、序号&确认序号三、头部长度四、保留位五、特殊标志位六、窗口大小七、校验和八、紧急指针九、可选选项十、数据 逻辑结构 如图&…

入门Java编程的知识点—>Http协议(day20)

了解http协议是什么掌握http请求信息、响应信息格式 项目目标&#xff1a; 实现本地客户端与服务器一问一答的请求与响应,了解http协议即可. 项目步骤&#xff1a; 服务器端代码编写 先在当前src文件下新建一个包: webserver,再该包下创建一个类Server&#xff0c;书写代码如…

带着耐心细心平常心和编程共舞

编程是什么&#xff1f;一个工具、一门技术还是一个爱好&#xff0c;不同的对待方法会带来不同的心态、产生不同的结果。编程需要扎实的基础、严密的思维和开阔的视角&#xff0c;新技术和框架日新月异&#xff0c;只有抱着科学、乐在其中的态度才能掌握高效的学习、实践方法。…

Matrix:重塑APM领域,以简驭繁的性能监控新纪元

在数字化转型的浪潮中&#xff0c;应用程序的性能监控&#xff08;APM&#xff09;已成为企业IT架构中不可或缺的一环。随着业务复杂度的提升和用户对体验要求的日益增高&#xff0c;如何高效、精准地监控并优化应用性能&#xff0c;成为了每个开发者和技术团队面临的重大挑战。…

机器学习(五) -- 监督学习(8) --神经网络1

系列文章目录及链接 上篇&#xff1a;机器学习&#xff08;五&#xff09; -- 无监督学习&#xff08;2&#xff09; --降维2 下篇&#xff1a; 前言 tips&#xff1a;标题前有“***”的内容为补充内容&#xff0c;是给好奇心重的宝宝看的&#xff0c;可自行跳过。文章内容被…

2.3导数与微分的基础与应用

1. 导数的基本概念 大家好&#xff0c;欢迎来到我们的数学大讲堂&#xff01;今天我们要聊聊一个有点酷又有点恐怖的东西——导数。别担心&#xff0c;不是让你在黑板上画曲线的那种&#xff0c;而是关于“变化率”的一种数学表达。 那么&#xff0c;什么是导数呢&#xff1f…

利用实用规模量子计算模拟宇宙中最极端的环境

华盛顿大学和劳伦斯伯克利国家实验室最近的研究展示了可扩展的技术&#xff0c;有朝一日可以实现最高能量下的基础物理实验模拟。 目录 核物理和高能物理的实用规模模拟 我们的模拟方法 Qiskit 如何使我们的实验成为可能 展望量子模拟技术的未来 粒子物理学的标准模型囊括了我们…

005.Python爬虫系列_浏览器开发者工具(详解)

我 的 个 人 主 页:👉👉 失心疯的个人主页 👈👈 入 门 教 程 推 荐 :👉👉 Python零基础入门教程合集 👈👈 虚 拟 环 境 搭 建 :👉👉 Python项目虚拟环境(超详细讲解) 👈👈 PyQt5 系 列 教 程:👉👉 Python GUI(PyQt5)文章合集 👈👈 Oracle数…

web渗透:XXE漏洞

XXE漏洞的原理 XXE&#xff08;XML External Entity Injection&#xff09;漏洞是由于应用程序在解析XML输入时&#xff0c;没有正确处理外部实体&#xff08;External Entities&#xff09;的引用而产生的安全漏洞。XML外部实体可以引用本地或远程的内容&#xff0c;攻击者可…

大学生房屋租赁推荐系统

一、项目概述 Hi&#xff0c;大家好&#xff0c;今天分享的项目是《基于推荐算法的大学生租房交流平台》。 对于初入社会的大学生来说&#xff0c;租房是一个复杂且繁琐的过程。租房交流平台可以提供房源浏览、信息分享等&#xff0c;帮助大学生快速找到合适的房源&#xff0…

随笔2优化算法

目录 2.1 线性规划&#xff08;Linear Programming&#xff09; 2.2 非线性规划&#xff08;Nonlinear Programming&#xff09; 2.3 整数规划&#xff08;Integer Programming&#xff09; 优化算法 是数学建模中非常重要的工具&#xff0c;用于在给定的约束条件下找到某个目…

HCIP笔记13-交换(2)

一台交换机上的一个接口可以映射多个不同的MAC地址&#xff0c;但一台PC上一个MAC只能对应一个接口 三层架构--冗余--二层的桥接环路 导致的问题&#xff1a; 1.广播风暴 2.对同一数据帧的反复拷贝 3.MAC地址表翻滚 这三个问题将不断占用硬件资源&#xff0c;逐渐导致设备…

大数据基础:实时数仓发展趋势

文章目录 实时数仓发展趋势 一、​​​​​​​实时数仓现状 二、​​​​​​​​​​​​​​批流一体 实时数仓发展趋势 一、​​​​​​​​​​​​​​实时数仓现状 当前基于Hive的离线数据仓库已经非常成熟&#xff0c;随着实时计算引擎的不断发展以及业务对于实时…