单细胞数据分析-Seurat 提取细胞亚群（特定几个细胞群）

单细胞数据分析
已有数据库
Seurat
提取细胞亚群（特定几个细胞群）
Dotplot

R脚本代码如下：

hv1 <- readRDS(file = "./singlecell/bairong/fastgenomics-Hv1h5ad_2000vf-202304231117.rds")
hv1
sub_class = c("Basal resting","Suprabasal","Multiciliated (nasal)","Multiciliated (nasal)","Multiciliated (non-nasal)","Club (non-nasal)","Club (nasal)","Goblet (nasal)","Goblet (bronchial)","Goblet (subsegmental)","AT1","AT2","AT2 proliferating")
# hv1_finest_sub
hv1_fs = hv1[,hv1@meta.data$ann_finest_level %in% sub_class]

saveRDS(hv1_fs, file = "./singlecell/bairong/fastg-hv1-subcells_202308081058.rds",compress = F)#快



# remotes::install_github("mojaveazure/seurat-disk", force = TRUE)
library(SeuratDisk)
library(patchwork)
library(dplyr)
library(Seurat)




hv1_fs <- readRDS(file = "./singlecell/bairong/fastg-hv1-subcells_202308081058.rds")

library(ggplot2)


Proteases = c("ST8SIA1")
Proteases = c("TMPRSS2")
Proteases = c("ACE2")
Proteases = c("FURIN")
Proteases = c("ST8SIA1","TMPRSS2")
Proteases = c("ST8SIA1","TMPRSS2","ACE2")
Proteases = c("ST8SIA1","TMPRSS2","ACE2","FURIN")


DotPlot(hv1_fs, features = Proteases, group.by="ann_finest_level" ,        
        cols = c("yellow","red"))+  RotatedAxis()