# # - eval = FALSE ------------------------------------------------------------ # 库(LOLA) # # regionDB = loadRegionDB(“LOLACore / hg19”)# # regionSetA = readBed(“lola_vignette_data / setA_100.bed”)# regionSetB = readBed(“lola_vignette_data / setB_100.bed”)# regionSetC = readBed(“lola_vignette_data / setC_100.bed”)# activeDHS = readBed(“lola_vignette_data / activeDHS_universe.bed ") ## ---- eval = FALSE ------------------------------------------------------------ # userSets =GRangesList (regionSetA regionSetB regionSetC) # locResults = runLOLA (userSets、activeDHS regionDB,核= 1 ) ## ---- eval = FALSE ------------------------------------------------------------ # locResults [1:10 ,] ## ---- eval = FALSE ------------------------------------------------------------ # locResults[秩序(meanRnk,减少= FALSE),] [1:20,] # locResults[秩序(maxRnk,减少= FALSE),] [1:20,] # locResults[= =“sheffield_dnase”集合,][秩序(maxRnk,减少= FALSE),) # # locResults [userSet = = 2][秩序(maxRnk,减少= FALSE),] [1:10,] # locResults [userSet = = 3][秩序(maxRnk,减少= FALSE),] [1:10], # # locResults [userSet = = 1][= =“sheffield_dnase”集合,][15,]# locResults [userSet = = 2][= =“sheffield_dnase”集合,][15,]# locResults [userSet = = 3][= =“sheffield_dnase”集合,][1:15 ,] ## ---- eval = FALSE ------------------------------------------------------------ # writeCombinedEnrichment (locResultsoutFolder= "lolaResults"， inclesplits =TRUE)