"站长,Mapping之后counts怎么合并成一个表?"
虽然旧文 但是非常有料~~~
在这之前先看下面的教程
总结 从零到壹:10元转录组分析小结~干货~ 然后,重点看批量处理数据的技巧~从零到壹:10元转录组分析 从零到壹:10元转录组分析~硬盘不够用咋办 从零到壹:10元~Mapping神器STAR的安装及用 从零到壹:从SRA下载到分析~纯干货 10元转录组分析:这次真的是干货了~灰常干
得到ReadsPerGene数据后
得到每个基因的Counts数之后,你需要将这些不同文件中的提取出来,以制备DEseq2所需要的原始文件,组数少的情况下很好吧,看好第几列、第几行,用R语言按照下面的命令就可以x<-Counts[-(1:4),2] #去掉的1到4行,选取第2列然后用cbind把所有x并在一起就OK了。
但是数量巨大怎么办
比如以下这样的300+样本
"少废话,来干货~"将R语言工作环境设置为这些文件所在文件夹注意这些文件夹中不能有其他文件如果你的样本是链特异性(Reverse)测序“啥是链特异性,需要解释的留言”用下面这组命令library(stringr)library(dplyr)fall <- dir()data.out<-df.use[1:4,]for (fnow in fall) {str <- str_sub(fnow, start = 4, end = 10)num <- as.numeric(str)df.read <- read.table(fnow)df.use <- data.frame(v1 = df.readV1,v4=df.readV4)V1,v4=df.rea)str_c<-str_c('SRR',str)colnames(df.use)<-c('V1',str_c)data.out <- full_join(data.out, df.use,by="V1")}data.out1<-data.out[-(1:4),-2] #这个是对data.out修整write.csv(data.out1, file = 'F:/out.csv')data.out1 就是DEseq2包中需要用的文件
相关文章
- Android怎么查看手机中的本地数据库
- java怎么输入字符_java怎么输入一个字符
- 内存频率有哪些?怎么看内存频率
- D2C 设计稿转代码是怎么实现的?自己做一个可行吗?
- Echarts漏斗图形的,怎么保障当数据是0的时候,漏斗图形一直存在。
- 怎么快速 把一个 监听事件 变成 一个 内部类 在编程一个新的 java文件
- 怎么消除苹果手机中点击输入框会出现一个半透明的灰色背景详解手机开发
- 学习Oracle如何编写函数(oracle函数怎么写)
- Linux自学指南:成为一个超级用户!(怎么自学linux)
- 自学Linux:一个全新的技术之旅(怎么自学linux)
- 怎么去转换任何系统调用为一个事件:对 eBPF 内核探针的介绍
- 让你快速上手:Linux系统的使用指南(linux 系统怎么用)
- 判断程序是否连接Redis的技术方法(怎么看是否连接redis)
- 测试Redis连通性5步简易指南(怎么测试redis连通性)
- 微擎用Redis让网络无忧(微擎怎么开启redis)
- 从BI到MySQL 链接一个改变着你的世界(bi怎么链接mysql)
- 如何使用ADS连接MySQL(ads怎么连接mysql)
- Redis储存优势点赞技术开发之道(点赞怎么使用redis)
- 如何在Oracle中实现循环(Oracle中怎么循环)
- 查看Oracle交集表一个新手指南(oracle交集表怎么看)
- Redis阻塞队列操作指南(redis阻塞队列怎么用)
- 如何在MySQL中正确地编写单个字段(mysql一个字段怎么写)
- 数Oracle上月数据获取的简单方法(oracle上月怎么取)