4.13 实战十三 | 10X | 小鼠乳腺上皮细胞
刘小泽写于2020.7.21
library(scRNAseq)
sce.mam <- BachMammaryData(samples="G_1")
sce.mam
# class: SingleCellExperiment
# dim: 27998 2915
# metadata(0):
# assays(1): counts
# rownames: NULL
# rowData names(2): Ensembl Symbol
# colnames: NULL
# colData names(3): Barcode Sample Condition
# reducedDimNames(0):
# altExpNames(0):
# 样本信息
sapply(names(colData(sce.mam)), function(x) head(colData(sce.mam)[,x]))
# Barcode Sample Condition
# [1,] "AAACCTGAGGATGCGT-1" "G_1" "Gestation"
# [2,] "AAACCTGGTAGTAGTA-1" "G_1" "Gestation"
# [3,] "AAACCTGTCAGCATGT-1" "G_1" "Gestation"
# [4,] "AAACCTGTCGTCCGTT-1" "G_1" "Gestation"
# [5,] "AAACGGGCACGAAATA-1" "G_1" "Gestation"
# [6,] "AAACGGGCAGACGCTC-1" "G_1" "Gestation"
依然是整合行名 + 添加染色体信息
library(scater)
rownames(sce.mam) <- uniquifyFeatureNames(
rowData(sce.mam)$Ensembl, rowData(sce.mam)$Symbol)