资源见:
https://github.com/reactome/ReactomeGSA-tutorials/blob/master/notebooks/Jerby_Arnon_Seurat.ipynb
Get B cells
Clusters representing B cells are identified based on their expression of CD19, MS4A1 (CD20), CD79A, and SDC1 (CD138).
# 加载数据
jerby <- readRDS("data/jerby_processed.rds")
# 绘图
options(repr.plot.width = 12, repr.plot.height = 12)
VlnPlot(jerby, features = c("MS4A1", "SDC1", "CD19", "CD79A"), ncol = 1)
image.png
选取亚群
# 选取亚群
b_cell_clusters <- c(6, 10, 22, 23, 26, 29, 34)
b_cells <- subset(jerby, idents = b_cell_clusters)
网友评论