R 分块处理列表并合并结果
例如,我有两个基因列表R 分块处理列表并合并结果,r,list,loops,R,List,Loops,例如,我有两个基因列表 gene_list_A $`STEARATE BIOSYNTHESIS I (ANIMALS)%HUMANCYC%PWY-5972` [1] "ELOVL1" "ACOT7" "ACSL1" "ACSL5" "ACSL4" "ELOVL6" "ACSL3" "ACOT2" "
gene_list_A
$`STEARATE BIOSYNTHESIS I (ANIMALS)%HUMANCYC%PWY-5972`
[1] "ELOVL1" "ACOT7" "ACSL1" "ACSL5" "ACSL4" "ELOVL6" "ACSL3" "ACOT2" "ACOT1" "ACSBG1"
[11] "ACSBG2" "SLC27A2" "ACOT4"
$`SUPERPATHWAY OF INOSITOL PHOSPHATE COMPOUNDS%HUMANCYC%PWY-6371`
[1] "PI4K2B" "MTMR14" "PTEN" "INPPL1" "PIK3CD" "PIK3C2G" "PIK3CB" "PIK3C2A" "PIK3CG" "PIK3C2B"
[11] "PLCZ1" "PPIP5K1" "PPIP5K2" "PIP5KL1" "PLCE1" "PIP4K2A" "PIP4K2B" "PIP4K2C" "SACM1L" "ITPK1"
[21] "IPMK" "OCRL" "ITPKB" "MINPP1" "ITPKC" "PLCB3" "PLCB4" "PIK3CA" "ITPKA" "PIK3C3"
[31] "PLCB1" "PLCB2" "PI4K2A" "TMEM55A" "IPPK" "TMEM55B" "MTMR3" "PIK3R4" "PIK3R3" "PIK3R2"
[41] "PIK3R1" "PIK3R6" "PIK3R5" "INPP5B" "INPP5A" "INPP5D" "INPP5J" "PLCG2" "PIP5K1A" "INPP5K"
[51] "PIP5K1B" "PLCG1" "PIP5K1C" "IP6K1" "CDIPT" "IP6K3" "IP6K2" "SYNJ2" "FIG4" "PIKFYVE"
[61] "SYNJ1" "PLCH1" "PLCH2" "PI4KB" "PLCD3" "PLCD4" "PLCD1"
名单很大
> length(gene_list_A)
[1] 10362
> length(gene_list_B)
[1] 5145
我通过以下函数处理这两个列表
gs.RNASeq <- 21196
gom.obj <- newGOM(gene_list_A,gene_list_B,gs.RNASeq)
Intersection <- getMatrix(gom.obj, name="pval")
gs.RNASeq
gs.RNASeq <- 21196
gom.obj <- newGOM(gene_list_A,gene_list_B,gs.RNASeq)
Intersection <- getMatrix(gom.obj, name="pval")
Intersection <- as.data.frame(matrix(NA,nrow = length(gene_list_A), ncol = length(gene_list_B)))
nSym <- length(gene_list_B)
chunkSize <- 10
for(i in 1:ceiling(nSym / chunkSize)){
Intersection[,] <- getMatrix(newGOM(gene_list_A,gene_list_B[((i-1)*chunkSize+1):min(nSym,(i*chunkSize))],gs.RNASeq),name= "pval")
}
dim(Intersection)