愿武艺晴小朋友一定得每天都开心
## 按cellcycle对象画多层嵌套柱状图,按4个样本分开
umap <- fasting_memory@reductions$umap@cell.embeddings %>%
data.frame() %>%
cbind(cellcycle = fasting_memory$cellcycle) %>%
cbind(celltype = fasting_memory$Celltype) %>%
cbind(Sample = fasting_memory$Sample)
umap <- umap[umap$celltype %in% c("HSPCs","CLP","Pro-B","Pre-B","Immature-B","Mature-B"),]
umap <- umap[umap$cellcycle %in% c("Cycling"),]
library(dplyr)
sample_umap <- umap %>% group_by(Sample) %>% sample_n(size = 436) #等数量抽样
table(sample_umap$celltype)
colors <- c("fed"="#309EFF","memory_10d"="#63B8FF","memory_35d"="#87CEFF","memory_66d"="#B0E2FF")
sample_umap$celltype <- factor(sample_umap$celltyp