edgeR中怎么实现两组间差异分析操作
发表于:2024-09-22 作者:千家信息网编辑
千家信息网最后更新 2024年09月22日,edgeR中怎么实现两组间差异分析操作,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。1. 读取文件需要读取基因在所有样本
千家信息网最后更新 2024年09月22日edgeR中怎么实现两组间差异分析操作
edgeR中怎么实现两组间差异分析操作,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。
1. 读取文件
需要读取基因在所有样本中的表达量文件,示例如下
gene_id ctrl-1 ctrl-2 ctrl-3 case-1 case-2 case-3geneA 14 0 11 4 0 12geneB 125 401 442 175 59 200
每一行为一个基因,每一列代表一个样本。读取数据的代码如下
# 读取表达量的表格counts <- read.table( "gene.counts.tsv", header=T, sep="\t", row.names=1, comment.char="", check.names=F)# 设置样本分组groups <- factor(c(1,1,1,2,2,2))# 构建edgeR中的对象y <- DGEList(counts=count,group=group)
2. 过滤count数很低的基因
根据CPM
表达量对基因进行过滤,代码如下
keep <- rowSums(cpm(y)>1) >= 2y <- y[keep, , keep.lib.sizes=FALSE]
3. 归一化
默认采用TMM
归一化算法,计算每个样本的 sizefactor, 代码如下
y <- calcNormFactors(y)
4. 进行差异分析
代码如下
design <- model.matrix(~group)y <- estimateDisp(y,design)et <- exactTest(y)
5. 提取结果
将差异分析的结果保存到文件中,代码如下
res <- et$tablewrite.table(res, "edgeR.xls", header = T, col.names = NA, sep = "\t" )
看完上述内容是否对您有帮助呢?如果还想对相关知识有进一步的了解或阅读更多相关文章,请关注行业资讯频道,感谢您对的支持。
代码
基因
样本
差异
分析
文件
结果
帮助
清楚
一行
代表
内容
对此
对象
数据
文章
新手
更多
知识
示例
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
中国网络安全知识竞赛2019
崇明区市场软件开发服务介绍
宝山区正规软件开发品质保障
工程项目软件开发计费
金鑫互联网金融科技
数据库开发工程师面试常见问题
山东微山网络安全案
服务器erlang
学生会网络技术部工作思路
福建中职网络安全比赛答案
计算机网络技术期末网络技术
网络技术和信息处理专业
小米帐号服务器出错
软件开发项目 分成
excel工作表 数据库
比较安全的服务器
网络安全书法作品内容
中国十大软件开发的公司
kegg数据库hsp
读三级数据库心得体会
php数据库导出命令
在线直播服务器
redhat服务器安装
售后服务规范监控服务器
网络安全整体解决方案公司
静安区参考网络技术服务质量保障
战术小队服务器不显示
网络安全手抄报四年级 上册
校园网络安全论文1500
数据库正整数用什么类型