edgeR的使用
来源:互联网 发布:centoa安装java 7u80 编辑:程序博客网 时间:2024/06/07 10:33
注:本文为作者个人理解,如有不同见解欢迎指正o(^∀^)o
edgeR包是进行RNA-seq数据分析非常常用的一个R包。该包需要输入每个基因关于每个样本的reads数的数据,每行对应一个基因,每一列对应一个样本。建议使用htseq-count进行统计,输出文件即可直接使用。如果需要算RPKM,需要自己统计基因长度信息。
第一步:构建 DGEList类变量
edgeR的大多数操作都是对 DGEList类型变量进行,所以第一步必须构建该类型变量,使用函数:DGEList
命令为:
y <- DGEList(counts=x)
这里的x就是每个行一个基因每列一个样本的reads数表达矩阵。
当然还需要一个分组信息的数据,需要自己建立,例如:
group <- c(1,1,2,2)y <- DGEList(counts=x, group=group)
第二步:基因筛选
对于在大多数样本中表达数量都很少的基因,需要进行过滤,这一步可以根据自己定义的标准过滤,edgeR推荐使用该包的CPM( count-per-million )值进行过滤,命令:
keep <- rowSums(cpm(y)>1) >= 2#至少在两个样本里cpm大于1y <- y[keep, , keep.lib.sizes=FALSE]
第三步:标准化
edgeR的标准化思想主要针对的是不同样本在建库时效应。这一点与RPKM不同,因为edgeR认为不同的基因对于所有样本的影响是相同的,所以不必考虑。因此为了消除这种建库时的效应,edgeR会更推荐你使用他的calcNormFactors函数,算出来的值叫做trimmed mean of M-values (TMM) ,命令为:
y <- calcNormFactors(y)y$samplesgroup lib.size norm.factorsSample1 1 10880519 1.17Sample2 1 9314747 0.86Sample3 1 11959792 1.32Sample4 2 7460595 0.91Sample5 2 6714958 0.83
第四步:聚类与热图
在表达量分析中经常会遇到根据基因表达量来对基因进行聚类分析的过程,edgeR建议使用logcpm这个指标来进行聚类和画热图,命令是:
logcpm <- cpm(y, prior.count=2, log=TRUE)
第五步:差异表达分析
edgeR的差异表达分析有两种方式,简单的一种是直接进行配对检验,命令是:
et <- exactTest(y, pair=c("A","B"))#检验A、B样本的差异表达基因topTags(et)
另一种算法比较复杂,其思想是:首先表达值需要服从某个分布,然后与这个分布不相符的基因才是差异表达基因。具体做法如下:
design <- model.matrix(~Subject+Treat)#设计一个实验设计矩阵fit <- glmFit(y, design)#构建一个基因表达值得广义线性模型lrt <- glmLRT(fit)#通过似然比检验得到最后的结果
- edgeR的使用
- edgeR(未完)
- GUID使用的使用
- ${}的使用
- 的使用
- @[],@()的使用
- ->的使用
- +、-、*、/、%的使用
- &&的使用
- %:的使用
- Makefile中使用$$的使用
- Makefile中使用$$的使用
- Red5的使用--使用日志
- prefernce的使用具体使用
- 【工具使用】git 的使用
- 关于使用AsyncTaskLoader的使用
- 使用CABasicAnimation的基本使用
- 使用RocketMQ的客户端使用
- tcl库、包、模块的创建与使用
- MySQL
- UITableView grouped之调整任意header、footer首尾高度
- Java常用类 集合框架
- 什么是实时操作系统和非实时操作系统?
- edgeR的使用
- Java语言基础
- SpringBoot和Swagger结合提高API开发效率
- 计算机小白突破自己的第一步
- esxi gui
- bzoj2318Spoj4060 game with probability Problem 期望DP
- 我们从那里来—子宫日记 Womb
- 异常和集合类型
- 升级windows10后 Oracle VM VirtualBox 无法启动导致 genymotion 启动失败