如何定量、提取表达某基因的细胞（及比例）

 一，提取表达某基因的所有细胞

expr <- sce_endo@assays$RNA
gene_expression <- expr %>% .['PLAT',] %>% as.data.frame() %>% t()
gene_expression <- as.data.frame(gene_expression)
colnames(gene_expression) <- 'PLAT'
gene_expression$cell <- rownames(gene_expression)
gene_expression_sel <- gene_expression[which(gene_expression$PLAT > 0),]
sce_select <- sce[,rownames(gene_expression_sel)]
dim(gene_expression_sel) / dim(gene_expression)

二，计算每个cluster中表达某基因的细胞比例

for (i in c(0:15)){
  rate = length(WhichCells(object = sce_endo,
                           expression = VCAM1 > 0,
                           idents = i
  )) /table(sce_endo$seurat_clusters)[[i+1]]
  print(rate)
}