R语言如何做FDR校正(21)
1.什么是FDR校正?
FDR (false discovery rate),即校正后的P值,中文一般译作错误发现率。在转录组分析中,主要用在差异表达基因的分析中,控制最终分析结果中,假阳性结果的比例。
为什么要用FDR?
在转录组分析中,如何确定某个转录本在不同的样品中表达量是否有差异是分析的核心内容之一。一般来说,我们认为,不同样品中,表达量差异在两倍以上的转录本,是具有表达差异的转录本。为了判断两个样品之间的表达量差异究竟是由于各种误差导致的还是本质差异,我们需要根据所有基因在这两个样本中的表达量数据进行假设检验。常用的假设检验方法有t-检验、卡方检验等。很多刚接触转录组分析的人可能会有这样一个疑问,一个转录本是不是差异表达,做完假设检验看P-value不就可以了么?为什么会有FDR这样一个新的概念出现?这是因为转录组分析并不是针对一个或几个转录本进行分析,转录组分析的是一个样品中所转录表达的所有转录本。所以,一个样品当中有多少转录本,就需要对多少转录本进行假设检验。这会导致一个很严重的问题,当数据量很大时,在单次假设检验中较低的假阳性比例会累积到一个非常惊人的程度。
2.分析前的数据准备
一列P值
0.029842334
0.038417236
0.039422132
0.043541069
0.047857277
0.052881786
0.059129016
0.061978623
0.097525808
0.101549223
3. R语言怎么做FDR校正
p.adjust(c(0.01,0.05,0.01,0.05,0.03,0.02), # P值列表method ="BH" # FDR校正的方法
)
# [1] 0.030 0.050 0.030 0.050 0.045 0.040
4. BioLadder生信云平台在线做FDR校正
不想写代码?可以用BioLadder生信云平台在线绘制FDR校正。
网址:https://www.bioladder.cn/web/#/chart/58
4.1 导入数据
在文本框中输入您的P值,每个P值单独占一行。
4.2 调整参数
可以调整计算FDR的方法。
FDR的计算是根据假设检验的P-value进行校正而得到的。一般来说,FDR的计算采用Benjamini-Hochberg方法(简称BH法)
4.3 提交查看
上传完数据,并调好参数后,点击右上方“提交”按钮,图像会出现在左侧。如果没有上传自己的数据,会使用默认的demo数据绘图。
提交成功后,FDR结果表格会显示在左上方,同时会在左下方显示数据的分布情况的柱形图。
4.4 下载
点击右下角的下载按钮,可以下载结果表格。
5.BioLadder生物信息在线分析可视化云平台
网址:https://www.bioladder.cn/
BioLadder目前上线了生物学分析最常用的50多个模块,主要包括以下四类:
**数据可视化:**箱线图,南丁格尔玫瑰图,韦恩图,UpSet图,饼图,词云图,小提琴图,FDR校正,弦图,柱形图等
**组学数据分析:**序列的多重比对,表达数据的CV曲线图,FDR校正,T-SNE,热图,相关性热图等,趋势分析的mFuzz,差异分析的火山图,富集分析的气泡图,修饰位点上下游模体分析的seqLogo,Motif热图等
**功能分析:**GO弦图,相互作用网络图,富集分析
**数据预处理:**归一化,补值,FDR校正,长宽表互换等
R语言如何做FDR校正(21)相关推荐
- 随机森林回归预测r语言_使用随机森林(R语言)做回归
引言 随机森林( random forest) 是一种基于分类树( classification tree) 的算法,它可以用于分类和回归,本文在这里以广西地区1990-2014共25年的GDP数据作 ...
- r语言怎么做经验分布_R语言如何做马尔科夫转换模型markov switching model
原文链接: R语言如何做马尔科夫转换模型markov switching modeltecdat.cn 假设 有时间序列数据,如下所示.经验表明,目标变量y似乎与解释变量x有关.然而,乍一看,y的水 ...
- R语言—自动做统计学表格
R语言-自动做统计学表格 今天为大家介绍一个超级好用,能为大家节省很多时间的模型回归结果展示的包-epiDisply包,我们不用再一个一个手动填写结果.整理结果,R自动输出可以用于论文发表的结果格式. ...
- CSDN ARIMA R语言_SPSS做不了主成分分析,那就换R语言实现全过程吧
SPSS只能完成主成分分析的一部分环节,主成分得分等计算尚需结合其他工具(如Excel)来完成,这对SPSS用户来说,是极不方便的.我们一定要很清楚的知道,SPSS可以一条龙做因子分析,但不能通过菜单 ...
- 【R语言 数据挖掘】R语言如何做关联规则?
一.前言 提到数据挖掘,我们第一反应就是之前听到的啤酒和尿不湿的故事,该故事就是典型的数据挖掘中的关联规则.购物篮分析区别于传统的线性回归的主要区别为,关联分析针对离散数据:下面我们利用R语言的aru ...
- 跟着Nature Microbiology学作图:R语言ggplot2做散点图添加拟合曲线和p值
本地文件 s41564-021-00997-7.pdf 论文 Protective role of the Arabidopsis leaf microbiota against a bacteria ...
- r语言怎么做经验分布_训练宝宝语言能力应该怎么做
怎样正确引导宝宝的说话能力 宝宝的语言发育分为三个阶段,分别为:发音.理解和表达. 新生儿一出生就会哭了,两个月大的宝宝就能发出喉音,三个月大的宝宝就能发出类似"啊啊"的声音了. ...
- R语言如何做独立样本t检验
t.test(x绩效总分x绩效总分 ~ x绩效总分 x分组) t.test(要检验的变量~分组变量)
- 求助帖:用R语言lmer做多层线性模型HML,请问变量之间的关系如图所示应该怎么用语句呢?
材料的relevance是自变量 RT和Will是因变量 被试自身变量包括ppt_int1和ppt_int2 材料自身变量是item_int 随机误差是error
最新文章
- HbuilderX中的git的使用 git HbuilderXgit HbuilderX 使用git
- seaborn heatmap热力图用法
- 【JavaSE04】Java中循环语句for,while,do···while-思维导图
- abap 字符串处理
- Intellij IDEA 的使用
- lucene创建索引
- 用 Visual Studio 和 ASP.NET Core MVC 创建首个 Web API
- 从jHiccup开始
- 华南理工计算机接口技术随堂练习_研究生考试计算机408跟845有什么区别?
- java实现远程系统控制_java实现电脑远程控制完整源代码
- 深度学习2.0-39.RNN训练难题-梯度弥散与梯度爆炸
- 计算机能否代替老师英语作文,雅思大作文范文:电脑不可取代老师
- ubuntu设置截屏热键(区域截屏)
- 详解String的intern方法
- 重庆财经职业学院计算机一级考试题,2021年重庆财经职业学院单招语文考试模拟试题库...
- web页面实现剪切板图片黏贴功能
- OpenJdk各版本、各操作系统下载,亲测好用。
- VOT中的EAO评判指标
- 【开源】一款PyQT+Pyserial开发的串口调试工具
- android 百度地图 自定义地图标注,百度地图自定义标注