3ca文件读取及注释信息添加
作者:
皮尔洛_dys | 来源:发表于
2023-01-11 14:27 被阅读0次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"]
本文标题：3ca文件读取及注释信息添加
本文链接：https://www.haomeiwen.com/subject/opwbcdtx.html
网友评论