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"]
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容