在R里面write.csv或者write.table时,quote默认是True,如果需要保存的文件中字符串两边不显示双引号的话,要设置quote=F
#本流程主要为了得到pbmc3k的细胞barcode信息,为了下一步过滤pbmc3k的BAM文件library(SeuratData)pbmc3k=pbmc3kpbmc3k_meta=pbmc3k@meta[,c("seurat_annotations")]pbmc3k_meta$cell_id=rownames(pbmc3k_meta)pbmc3k_meta$`grouping_level_1`=c(rep("Control",1350),rep("Case",1350))pbmc3k_meta=pbmc3k_meta[,c("cell_id", `grouping_level_1`, "seurat_annotations")colnames(pbmc3k_meta)=c("cell_id","grouping_level_1","grouping_level_2")pbmc3k_meta$cell_id=paste0('CB:Z:',pbmc3k_meta$cell_id,'-1')write.table(pbmc3k_meta,'pbmc3k_meta.tsv',row.names = F, col.names = F,quote = F)write.table(filter,'filter.txt',row.names = F, col.names = F,quote = F)
后面根据filter文件过滤BAM文件的教程见:Is there way to filter the BAM file produced by 10x pipelines with a list of barcodes?
