# #——回声= FALSE -------------------------------------------------------------- 美元knitr: opts_chunk集(错误= FALSE,警告= FALSE,消息= FALSE)图书馆(BiocStyle) set.seed (0 ) ## ----------------------------------------------------------------------------- 库(S4Vectors) df < DataFrame (x = rnorm (10000), y = rnorm(10000),基因= paste0(“GENE_”,示例(100、10000、取代= TRUE),细胞= paste0(“CELL_样本(10000,取代df = TRUE))) ## ----------------------------------------------------------------------------- 库(BumpyMatrix)垫< - splitAsBumpyMatrix (df (c(“x”、“y”)],行= df基因美元,美元列= df细胞)垫垫[1][[1 ]] ## ----------------------------------------------------------------------------- 选择< - df (1:10 0,) smat < splitAsBumpyMatrix(选择[c(“x”、“y”)],行=选择基因美元,美元列=选择细胞,稀疏smat = TRUE) ## ----------------------------------------------------------------------------- 暗(垫)dimnames(垫)rbind(垫,垫)cbind t(垫(垫子,垫子) ) ## ----------------------------------------------------------------------------- 垫(c(“GENE_2”、“GENE_20”)]垫(“GENE_10[1:5]垫 ",] ## ----------------------------------------------------------------------------- 出去了。X <- mat[，，" X "]退出。x。x (1 ] ## ----------------------------------------------------------------------------- pos < -。X > 0 pos[，1] shift <- 10 * out。X + 1 shift[，1]出来。Y <- mat[，，" Y "]大于- out。X < out。Y greater[，1] diff <- out。Y - out。x diff [1 ] ## ----------------------------------------------------------------------------- 我< -垫(,' x '] > 0 &垫(,' y '] > 0我[1]子< -垫[我]子子[1 ] ## ----------------------------------------------------------------------------- 垫(,“x”)垫(“x”,放弃= FALSE ] ## ----------------------------------------------------------------------------- 垫(1,1,' x ']垫(1,1,“x”。dropk = FALSE)垫(1,1,“x”,放弃= FALSE)垫(1,1,“x”。dropk = TRUE, = FALSE ] ## ----------------------------------------------------------------------------- 复制< -垫复制(,“x”)< -副本(,' x '][1 * 20副本 ] ## ----------------------------------------------------------------------------- 意思是(out.x)[1:5, 1:5) #矩阵var (out.x)(1:5, 1:5) #矩阵 ## ----------------------------------------------------------------------------- 1:5,分位数(out.x)[1:5]范围(out.x)(1:5、1:5)# #----------------------------------------------------------------------------- pmax(。x, out.y ) ## ----------------------------------------------------------------------------- sessionInfo ()