读取表达矩阵创建SeuratObject

Bioinfoer

已于 2022-11-11 15:33:55 修改

阅读量1.5k

点赞数

分类专栏：生物信息学文章标签： r语言大数据

于 2022-11-03 11:11:57 首次发布

本文链接：https://blog.csdn.net/xiaozheng1213/article/details/127666377

版权

# R code
library(Seurat)
mat_5 <- read.table("GSM4037992_peripheral_donor_5_enriched_expression.tsv.gz")
row.names(mat_5) <- mat_5[,1]

#保存原始cluster labels
orig_cluster_labels5 <- mat_5$final_cluster_labels
mat_5 <- mat_5[,-c(1:3)]

#需要转置一下，行是细胞barcode
peripheral_donor_5 <- CreateSeuratObject(counts=t(mat_5),project="peripheral_donor_5")
peripheral_donor_5 <- AddMetaData(object = peripheral_donor_5, metadata = orig_cluster_labels5, col.name = "orig_cluster_labels")

mat_7 <- read.table("GSM4037994_peripheral_donor_7_enriched_expression.tsv.gz")
row.names(mat_7) <- mat_7[,1]
orig_cluster_labels7 <- mat_7$final_cluster_labels
mat_7 <- mat_7[,-c(1:3)]
peripheral_donor_7 <- CreateSeuratObject(counts=t(mat_7),project="peripheral_donor_7")
peripheral_donor_7 <- AddMetaData(object = pe