r语言dataellipse_几行R语言代码搞定菌群与环境因子或临床指标相关性的可视化...
相关性分析是生物信息学中常用的分析方法,可以用来分析菌群与菌群的关联,菌群与因子的关联等等。本文使用R语言内置函数cor()计算变量之间的相关系数,并用corrplot包进行可视化。(本文测试数据为R语言内置数据集mtcars)。
计算相关性矩阵
运用R语言内置函数cor()来计算相关系数:cor(x,method=c("pearson","kendall","spearman")),其中x是数据框矩阵,而pearson,kendall和spearman是计算相关性的三种方法,函数默认为pearson。
可视化相关系数矩阵
本文使用corrplot包进行相关性矩阵的可视化,corrplot(https://github.com/taiyun/corrplot)是发表在github上的开源R包。
具体过程如下:
install.packages("corrplot") #安装corrplot包
data(mtcars) #加载数据集
mydata
head(mydata, 7) #查看数据前7行
cordata
round(cordata, 2) #保留两位小数,得到相关性矩阵
##corrplot()函数进行相关矩阵的可视化
该函数通过颜色深浅可视化显著相关程度。参数主要有:
corrplot(corr,method=c("circle","square","ellipse","number","shade","color","pie"),type=c("full","lower","upper"),order=c("original","AOE","FPC","hclust","alphabet"),tl.cex=1,tl.col="red",order=c("original","AOE","FPC","hclust","alphabet"),hclust.method=c("complete","ward","ward.D","ward.D2","single","average","mcquitty","median","centroid"),tl.cex = 1,tl.col="red",...)
接下来简单介绍几个参数用法及展示其可视化:
1.矩阵类型,通过type参数,“upper”,“lower”, “full”,分别显示上三角,下三角还是全部。
corrplot(cordata, type="lower") #下三角矩阵
corrplot(cordata, type="full",method="circle") #全矩阵
2.method参数设置图形可视化的形状,method=c("circle","square","ellipse","number","shade","color","pie"),共有7种形状。
corrplot(cordata,type="full",method="ellipse") #设置图形为椭圆
corrplot(cordata,type="full",method="pie") #设置图形为饼图
3.通过函数corrplot.mixed绘制图形数值混合矩阵。
corrplot.mixed(cordata) #绘制图形数字混合矩阵
corrplot.mixed(cordata,lower="ellipse",upper="circle",tl.col="black") #下三角为圆形,上三角为椭圆
4.order参数将结果按照不同的聚类方式排序,addrect参数添加分组,具体如下所示。
corrplot(cordata,order="hclust",addrect=2) #以hclust聚类,分成2组
corrplot(cordata,order="hclust",addrect=2,col=heat.colors(100)) #修改颜色
小锐课堂,收获颇丰
更多好货,敬请期待
我们下期再见!!!
r语言dataellipse_几行R语言代码搞定菌群与环境因子或临床指标相关性的可视化...相关推荐
- 几行代码搞定Flash应用的多语言实时切换问题
在<WPF/Silverlight的数据绑定设计的真糟糕>文中批评了WPF/SL的数据绑定机制,拿的是Flex的数据绑定与其对比.本文介绍应用Flex的数据绑定来解决一个实际问题--多语言 ...
- python做事件研究法_35行代码搞定事件研究法(上)
作者简介: 祝小宇,个人公众号:大猫的R语言课堂 这期大猫课堂将会教大家如何用35行R代码写出最有效率的事件研究法. 注意,本代码主要使用data.table完成,关于data.table包的相应知识 ...
- 35行代码搞定事件研究法(下)
作者简介: 祝小宇,个人公众号:大猫的R语言课堂 前文推送: 35行代码搞定事件研究法(上) Hello亲爱的小伙伴们,上期已经讲到如何对单一事件日计算超额收益,本期将会教大家如何针对多个股票多个事件 ...
- 35行代码搞定事件研究法(上)
作者简介: 祝小宇,个人公众号:大猫的R语言课堂 这期大猫课堂将会教大家如何用35行R代码写出最有效率的事件研究法. 注意,本代码主要使用data.table完成,关于data.table包的相应知识 ...
- 万万想不到 10行代码搞定一个决策树
01决策树模拟实验 文章目录 01决策树模拟实验 要求 决策树简单介绍 搭建环境 产生数据集 划分训练集和测试集 生成决策树 Cross-Validation法 可视化决策树 10行代码搞定决策树 要 ...
- table中加表单元素每行怎么验证_Validform 一行代码搞定整站的表单验证 - 文章
Validform 一行代码搞定整站的表单验证,为什么能如此方便?插件的核心思想就是把所有的验证条件及验证提示信息绑定到每个表单元素,让验证代码在执行时只是核对表单下各元素的值是否跟绑定的验证条件相符 ...
- resnet50代码_13、SOTA论文实践-学习ResNet(80行代码搞定残差backbone网络)
0.论文 Camera Distance-aware Top-down Approach for 3D Multi-person Pose Estimation from a Single RGB I ...
- 零基础如何用 15 行 Python 代码搞定网易云热门歌单?
健身.吃饭.敲代码:等车.撸猫.下午茶--若能佐以合适的音乐当"配餐",总是惬意非常.本文就将带你爬一爬网易云的那些热门歌单! 作者 | 上海小胖 责编 | 仲培艺 心情好或心情坏 ...
- 国外stripe支付,超简单几行代码搞定
国外stripe支付,超简单几行代码搞定 海外的项目 需要stripe支付 很简单 几行代码 先加入依赖: compile 'com.stripe:stripe-android:6.1.2' 总共两种 ...
最新文章
- 2020最新直播源地址下载txt_TXT追书免费小说app安卓版下载-TXT追书免费小说最新版下载v5.0.0...
- 网站 html 中英文切换 - API 总结篇
- 著名Youtuber “NELK “发布有关比特币购买游艇的视频
- A2W和W2A :很好的多字节和宽字节字符串的转换宏
- 假如你有20万用来创业,你最想进入什么行业?
- Java 集合系列04之 fail-fast总结(通过ArrayList来说明fail-fast的原理、解决办法)
- DELPHI学习---简单类型
- RabbitMQ学习——生产者与消费者入门例子
- [原]ActiveReport6 for net使用(一)
- android禁止下拉状态栏
- 计算机专业英语词汇1500词
- usb打印机命令_自己动手不求人 日常打印机故障全知晓
- npm安装依赖包 版本冲突怎么办 --legacy-peer-deps的正确使用方法
- Vray和PS的冬景效果图制作教程蓝海创意云渲染
- 大厂Offer拿到手软啊!隔壁都馋哭了
- 面试部分梳理 - 计算机网络
- PHP app登录qq,APP配置QQ登录的详细教程
- jvm.option是什么,它是如何加载的
- 2018前端校招总结
- 怎样提高英语思维能力?
热门文章
- Nsight Compute使用出现does not have permission to profile错误如何解决
- Si8621AB-B-IS兼容品2Pai荣湃π122M30 3.0kV10Mbps 双通道数字隔离器
- PCI Express学习篇---Power Management(一)
- 模拟CMOS 第2章 MOS器件物理基础
- 【图像分割】扩散张量成像(DTI)脑白质纤维束成像【含Matlab源码 2174期】
- sql-labs 闯关 21~25
- android 控件上下排列_android 实现控件左右或上下抖动教程
- 写一篇高考作文 | 教育研究的创新之法:移用、化用和独创
- 基于HTK中文普通话语音识别DEMO,LilyMandarinSpeechRecognition V1.0
- 论文基于RF的3D骨架(RF-Based 3D Skeletons)--关于4D卷积分解的理解