微生物测序分析LEfSe
LEfse分析定义 LEfse分析即LDA Effect Size分析,可以实现多个分组之间的比较,还进行分组比较的内部进行亚组比较分析,从而找到组间在丰度上有显著差异的物种(即biomaker);
主要是通过非参数因子Kruskal-Wallis秩和检验来实现的。
运行LEfSe软件主要分三大步骤:第一步:需要把普通的物种、基因等等的丰度信息的表格转化成LEfSe识别的格式。这一步会生成.in结尾的文件
第二步:这一步也是最关键的一步,统计显著差异的biomarker、统计子组组间差异、统计effect sizes(LDA score),会生成.res格式的文件。如下图所示Step1:两组或两组以上的样本中采用的非参数因子Kruskal-Wallis秩和检验检测出biomarker。
Step2:基于上步的显著差异物种基因,进行两两组之间的Wilcoxon秩和检验,检测出组间差异。
Step3:线性判别分析(LDA)对biomarker进行评估差异显著的物种的影响力(即LDA score),最终获得biomarker。第三步:基于第二大步的数据,绘制各种图片。如下图所示
LDA值分布柱状图:
展示了LDA score大于设定值有差异的物种,即具有统计学差异的biomaker。展现不同组中丰度有显著差异的物种,柱状图的长度代表显著差异物种的影响大小;
进化分支图:
由内至外辐射的圆圈代表了由门至属(或种)的分类级别。在不同分类级别上的每一个小圆圈代表该水平下的一个分类,小圆圈直径大小与相对丰度大小呈正比。
着色原则:无显著差异的物种统一着色为黄色,差异物种 Biomarker跟随组进行着色,红色节点表示在红色组别中起到重要作用的微生物类群,绿色节点表示在绿色组别中起到重要作用的微生物类群,其它圈颜色意义类同。图中英文字母表示的物种名称在右侧图例中进行展示。
biomaker在不同组各样本中的丰度比较图:
将biomaker丰度最高的样本的丰度设定为1,其他样品中该 biomarker 的丰度为相对于丰度最高样品的相对值。
微生物测序分析LEfSe相关推荐
- linux下微生物软件,微生物多样性专题 | 扩增子测序分析实战(三)软件安装
上一篇我也想装个linux系统!讲到linux系统的安装,实则就是为了本期准备的,因为测序数据分析涉及的大部分软件都是基于linux的,目前还没有linux系统或者打算只在windows下用R去做一些 ...
- 121个人电脑搭建微生物组分析平台(Win/Mac)
121个人电脑搭建微生物组分析平台(Win/Mac) 本节作者:刘永鑫 中国科学院遗传与发育生物学研究所 版本1.0.3,更新日期:2020年8月27日 本项目永久地址:https://github. ...
- iMeta:青岛大学苏晓泉组开发跨平台可交互的微生物组分析套件PMS(全文翻译,PPT,视频)...
Parallel-Meta Suite:跨平台可交互的微生物组快速分析套件 Parallel-Meta Suite: Interactive and rapid microbiome data ana ...
- MPB:扬州大学王梦芝组-反刍动物瘤胃原虫18S rRNA测序分析技术
为进一步提高<微生物组实验手册>稿件质量,本项目新增大众评审环节.文章在通过同行评审后,采用公众号推送方式分享全文,任何人均可在线提交修改意见.公众号格式显示略有问题,建议电脑端点击文末阅 ...
- QIIME 2教程. 04人体微生物组分析Moving Pictures(2020.11)
文章目录 QIIME 2用户文档. 4人体各部位微生物组 本节视频视频教程 启动QIIME2运行环境 样本元数据 下载和导入数据 拆分样品 序列质控和生成特征表 方法1. DADA2 方法2. Deb ...
- QIIME 2用户文档. 4人体各部位微生物组分析实战Moving Pictures(2018.11)
文章目录 前情提要 QIIME 2用户文档. 4人体各部位微生物组 启动QIIME2运行环境 样本元数据 下载和导入数据 拆分样品 序列质控和生成特征表 方法1. DADA2 方法2. Deblur ...
- 斯坦福大学统计系教授带你玩转微生物组分析
Susan Holmes 主要作品简介 1. 微生物组可重复分析R包 2. 微生物组数据统计方法讨论 3. 非聚类OTU算法dada2 课程推荐 可重复计算 孕妇微生物组的时空变异 孕妇阴道微生物组分 ...
- 二代三代转录组测序分析实战班
本文原创"生信宝典"公众号,作者陈同. 转录组大家都很熟悉了,我们之前也有几篇介绍: 转录组分析的正确姿势 39个转录组分析工具,120种组合评估(转录组分析工具哪家强-导读版) ...
- 表观调控高通量测序分析培训开课啦
在广大粉丝的期待下,<生信宝典>联合<宏基因组>在2018年4月14日在北京鼓楼推出<ChIP系列高通量测序分析专题培训>,为大家提供一条走进生信大门的捷径.为同行 ...
- 单细胞测序分析之小技巧之for循环批量处理数据和出图
"harmony"整合不同平台的单细胞数据之旅生物信息学习的正确姿势 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChI ...
最新文章
- Test2 unit4
- matlab 折线图_MATLAB作图实例:13:折线图和条形图的双轴图
- 《方舟生存进化》单机存档损坏解决方法
- matlab怎么复制相同的列,怎样将数组中某一列数据相同的行提取出来?
- 2019 年备受争议的 Facebook 能否走出去年的阴影?| 畅言
- 【论文笔记】From Facial Parts Responses to Face Detection: A Deep Learning Approach
- 盘点下半年EI国际会议,值得收藏!
- 3个方法解决百度网盘限速 (2019-08-04 更新)
- 微信开放平台----微信扫码登录
- IT 接口对接:足迹第十二步接口对接的定义(接口对接分三种:中间库方式的接口对接,Rest格式URL对接和HTTP格式URL对接;)
- uva10105 - Polynomial Coefficients(多项式系数)
- 【论文阅读】DeepIM: Deep Iterative Matching for 6D Pose Estimation
- Android 自定义View之咖啡杯动画
- 安装Glusterfs
- 屏幕录制软件camtasia studio 8序列号激活
- 【Flutter 异步编程 - 捌】 | 计算耗时? Isolate 来帮忙
- scrapy爬取伯乐在线文章
- 创建 Laravel 项目
- 26期20180702 shell特殊符号 cut sort_wc_uniq trr_tr_spli
- Python zipfile 简单使用(暴力破解)