DoubletFinder去双胞根据不同矩阵质量耗时长短差异很大,所以很多样本的建议多个任务并行跑,不然要等很久很久很久(亲测)。有时候用DoubletFinder也不能完全去掉所有双胞,但是双胞一般会自动聚成一个亚群,可以手动删除这种亚群。
对单个样本取双胞。返回的是一个Seurat对象。 obj <- rm_doublet(name,input) 可保存为RDS格式。 saveRDS(obj,paste0(name,".rds")) 建议写一个R脚本传参批量投递脚本得到RDS（这部分就不写了，如果有需求可留言），之后也可批量读取RDS，具体可参考下面的函数：以上方式适合多样本（...