1. par(mar = c(4, 8, 2, 1))#展示在细胞中表达量最高的gene
    2. C=sce.all.filt@assays$RNA@counts
    3. dim(C)
    4. C=Matrix::t(Matrix::t(C)/Matrix::colSums(C)) * 100
    5. # 这里的C 这个矩阵，有一点大，可以考虑随抽样
    6. C=C[,sample(1:ncol(C),1000)]#对C取1000个sample即cell
    7. most_expressed <- order(apply(C, 1, median), decreasing = T)[50:1]#对C的每一行进行取中位数
    8. pdf("TOP50_most_expressed_gene.pdf",width=14)
    9. boxplot(as.matrix(Matrix::t(C[most_expressed, ])),
    10.         cex = 0.1, las = 1, 
    11.         xlab = "% total count per cell", 
    12.         col = (scales::hue_pal())(50)[50:1], 
    13.         horizontal = TRUE)
    14. dev.off()