单细胞多文库数据批读取-R

单细胞一个样本通常构建多个测序文库,所以拿到的下机文库数据通常需要合并起来分析,在R中通过dir()函数进行文件检索,配合lapply()进行数据读取,再merge()在一起会快很多

setwd(".")
infiles <- dir(".",pattern = "*_raw.rds", full.names = TRUE,recursive = T)
scRNAlist <- lapply(infiles, readRDS)#读取数据为list
scRNA <- merge(scRNAlist[[1]],scRNAlist[2:length(scRNAlist)])
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容