2024-02-23
2024-01-19
2023-03-30
##细胞注释##使用VlnPlot画marker小提琴图VlnPlot(scobj, features = c("CD8B","CD3D","CD8A")) # CD8 TVlnPlot(scobj, features = c("CD3D","CD3E","CD4","IL7R"))# CD4 TVlnPlot(scobj, features = c("MS4A1", "CD79A","CD79B")) # BVlnPlot(scobj, features = c("PPBP", "PF4")) # PlateletVlnPlot(scobj, features = c("S100A9","CD14")) # CD14+ MonocyteVlnPlot(scobj, features = c("CDKN1C","FCGR3A")) # CD16+ MonocyteVlnPlot(scobj, features = c("CD74","CD1C")) # DCVlnPlot(scobj, features = c("CD3D","NCAM1","NKG7")) # NKVlnPlot(scobj, features = c("CYTL1","SPINK2")) # hspcVlnPlot(scobj, features = c("MKI67")) # 增殖细胞VlnPlot(scobj, features = c("HBM")) # 红细胞FeaturePlot(scobj,pt.size = 0.5,reduction = "tsne",c("CD8B","CD4","NCAM1","MS4A1","CD14","FCGR3A"),order = TRUE,ncol=3)
DimPlot(scobj, reduction = "tsne", group.by ="celltype",pt.size = 1,label.size = 5,label = T,raster = F)
library(ggplot2)DoHeatmap(subset(scobj, downsample = 200),features = marker_genes,size = 3) +theme(axis.text.y = element_text(size = 10))DoHeatmap(subset(scobj, downsample = 200),features = marker_genes,size = 3,group.by = "celltype",assay = "RNA",group.colors = c("#C77CFF","#7CAE00","#00BFC4","#F8766D","#AB82FF","#90EE90","#00CD00","#008B8B","#FFA500"))+scale_fill_gradientn(colors = c("navy","white","firebrick3")) +theme(axis.text.y = element_text(size = 8))
