# #——选项,包括= FALSE,回声= FALSE --------------------------------------- 库(BiocStyle) knitr: opts_chunk美元集(= FALSE警告、错误= FALSE消息= FALSE ) ## ----------------------------------------------------------------------------- 库(SingleCellExperiment)计数< -矩阵(rpois(100年,λ= 10),ncol = 10, nrow = 10) sce < - SingleCellExperiment(计数)altExp(,“飙升”)< - SingleCellExperiment(矩阵(rpois(20,λ= 5),ncol = 10, nrow = 2)) altExp (,“蛋白质”)< - SingleCellExperiment(矩阵(rpois(50,λ= 100),ncol = 10, nrow = 5)) altExp(,“CRISPR”)< - SingleCellExperiment(矩阵(rbinom (p = 0.1, 80 1), ncol = 10, nrow = 8))预计 ## ----------------------------------------------------------------------------- totalCount < -函数(x, i = 1,乘数= 1,subset.row = NULL){垫< -试验(x,我)如果(! is.null (subset.row)){垫< -垫(子集。行,,= FALSE)} colSums(垫)*乘数 } ## ----------------------------------------------------------------------------- 总数< - applySCE (sce、有趣= totalCount)总数 ## ----------------------------------------------------------------------------- 总数。手册< -列表(totalCount (sce),峰值= totalCount (altExp(,“飙升”)),蛋白质= totalCount (altExp(,“蛋白质”)),CRISPR = totalCount (altExp(,“CRISPR”)))stopifnot(相同(totals.manual总数 )) ## ----------------------------------------------------------------------------- totals10。手册< -列表(totalCount (sce乘数= 10),峰值= totalCount (altExp(,“飙升”),乘数= 10),蛋白质= totalCount (altExp(,“蛋白质”),乘数= 10),CRISPR = totalCount (altExp(,“CRISPR”),乘数= 10 ) ) ## ----------------------------------------------------------------------------- totals10。apply <- applySCE(sce, FUN=totalCount, multiplier=10) stopifnot(相同的(totals10。运用,totals10.manual )) ## ----------------------------------------------------------------------------- totals10。lapply <- lapply(c(List(sce)， altExps(sce))， FUN=totalCount, multiplier=10) stopifnot(相同的(totals10.)运用,totals10.lapply )) ## ----------------------------------------------------------------------------- 总数。custom <- applySCE(sce, FUN=totalCount, multiplier=10, ALT.ARGS=list(Spike=list(sub - set.row=2)， Protein=list(sub - set.row=3:5)))自定义 ## ----------------------------------------------------------------------------- 头。sce <- applySCE(sce, FUN=head, n=5)头。南加州爱迪生公司 ## ----------------------------------------------------------------------------- altExp head.sce altExp(头。sce，“蛋白质”)altExp(头。预计,“CRISPR ") ## ----------------------------------------------------------------------------- head.sce.list < - applySCE (sce、有趣=头,n = 5,简化= FALSE) head.sce.list ## ----------------------------------------------------------------------------- 手册。head <- head(sce, n=5) altExp(手动。head， "Spike") <- head(altExp(sce， "Spike")， n=5) altExp(手动。head, "Protein") <- head(altExp(sce, "Protein"), n=5) altExp(manual.head, "CRISPR") <- head(altExp(sce, "CRISPR"), n=5) manual.head ## ----------------------------------------------------------------------------- sessionInfo()