3ca文件读取及注释信息添加

setwd("G:\\3CA\\Breast") 
library(Seurat)
# https://hbctraining.github.io/scRNA-seq/lessons/readMM_loadData.html
library(Matrix)
library(R.utils)
matrix <- readMM("Exp_data_UMIcounts.mtx")
barcode <- read.table("Cells.csv",sep = ",",header = TRUE,row.names = 1)
genes <- read.table("Genes.txt")
row.names(matrix) <- genes$V1
colnames(matrix) <- row.names(barcode)
scar.atlas <- CreateSeuratObject(matrix,project = "SCAR_Atlas_0020")
scar.atlas@meta.data$cell_type <- barcode[row.names(scar.atlas@meta.data),"cell_type"]
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容